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THERMAL PROPERTIES AND GROWTH OF CRYSTALS 


The following six contributions were read at the Society’s meeting at 
Southampton on 19th December 1951. 


The Growth from the Melt of Single Crystals of 
Low Melting Point Metals 


By A. J. GOSS anp S. WEINTROUB 


University of Southampton 


MS. received 10th March 1952 ; read before the Physical Society at Southampton 
on 19th December 1951 


ABSTRACT. The results of a series of experiments on the preparation of tin, lead, zinc, 
and bismuth single crystals in the form of rods are given. ‘The effects on the crystallization 
of various thermal factors, in particular the temperature gradient and the rate of growth, are 
studied and the results for the four metals are compared. 


SEN ERO DWC ON 
sl HE study of the physical properties of metals is greatly simplified by the 
use of specimens in the form of single crystals. Various techniques have 
been proposed and developed during the past forty years for the growth 
of such crystals but a close examination of the methods reveals that they are 
still largely empirical. Holden (1949), to whom we are indebted for the most 
recent comprehensive review of single crystal techniques, makes it clear that 
' the information available as to the conditions governing the successful growth 
' of metal single crystals is but meagre and very often contradictory. Several 
_ investigators, however, consider that the purity of the metal, the temperature 
gradient G in the specimen, and the rate R of growth of the single crystals, are 
controlling factors, but no systematic investigation of the influence of these or 
other factors on the successful growth of metal single crystals has as yet been 
made. Experiments to provide information on this subject are being carried 
out in this laboratory and some of the results so far obtained for pure tin*, lead, 

zinc and bismuth are described and discussed in this paper. 
The horizontal travelling furnace method of Andrade and Roscoe (1937) 
was used and specimens of the four low-melting point metals mentioned above 
* The results for tin have been briefly reported in Nature, Lond., 1951, 167, 349. 
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were crystallized, by slow solidification from the melt, into single crystal rods | 
approximately 5mm in diameter and some 20-30cm in length. About | 


100 specimens of each metal were investigated. The temperature gradient in 


the solidifying metal was varied over as wide a range as possible, and various _ 


rates of growth, from 0-06 to 30mm/min, were used. ‘The variation in the 
temperature of the metal was recorded by means of a thermocouple embedded 
in a dummy specimen and the average temperature gradient at the solid—liquid 
boundary deduced. This is the gradient which is quoted in the tables. The 
rods, after crystallization, were carefully examined, the crystalline structure 
determined and the crystallographic orientations of the single crystals measured. 
The results, together with other relevant details, are listed and discussed in the 
following sections. Details of the experimental arrangements, theoretical 


considerations and the influence of varying amounts of added impurities are 


omitted. It is hoped to describe these later. 


§2. THE INFLUENCE OF G AND R ON THE FORM OF THE CRYSTALS 
It is well known that single crystals are not perfect but frequently contain 


perceptible imperfections of various kinds which are referred to in the literature 


by such names as lineage, mosaics, striations, etc. The rods, after crystallization, 


were first examined visually and the imperfections seen carefully noted. The — 
rods were then chemically etched and a sketch made of the crystalline structure 


made visible by the etching. When 10cm or more of a rod was found to consist 
of one predominant single-crystal structure, that rod was classified as single. 


The orientation of the cylindrical axis of the rod, relative to the crystallographic. 
axes of the crystals formed in the rod, was determined, by means of a specially, 


designed optical goniometer, by the method used by Chalmers (1935). 


For the four metals studied, it was found that, on an average, though 


exceptions did occur, the crystalline forms of the rods of any one metal which 
were produced with he same values of G and R followed the same general 
pattern. ‘This ‘average’, over 10-20 or more specimens, is the subject of the 
rows in the tables 1-4. In the columns are listed the number of specimens 
examined in the particular group, the percentage of single crystals obtained, 
the percentage with imperfections and the values of G and R used. 
(i) Tin 

In table 1 the results for tin are summarized. The metal was 99-997% pure 
and was obtained from Messrs. Johnson, Matthey & Co. Ltd. In most cases 
the first few centimetres of the rod to solidify consisted of small crystals and 
the remainder of the rod usually formed a single crystal. Varying G had little 


or no effect on the crystalline form but varying R altered the results considerably. 
The change of shape of the temperature-time trace at solidification was also 


altered by the variation in R (see §4(i)). 

The polycrystals within the first centimetre of the rods were 1 mm? or less 
in surface area and larger crystallites appeared as the distance from the end of 
the rod increased. ‘The inter-crystalline boundaries were generally jagged, 


especially near the beginning of the rod. Experiments performed with a 


thermocouple embedded within this portion of the specimens indicated that 


supercooling occurred in this case. For the slowest speeds, the first section of 


the rod consisted very often of a single crystal but of different orientation from 
that of the remainder. 


ee 
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For rates of growth of 0-5 mm/min and 20mm/min, apart from the first 
| few centimetres, the remainder consisted of a single crystal for which the 
(9 orientation was constant throughout the whole length of the rod. Three main 
_ types of imperfections were observed in the single crystal portions: (a) small 
inserts of crystals of different orientations; (b) two or more crystals within the 
_ length and (c) lineage, i.e. thin strips of crystallite, less than 1 mm in width. The 

crystals in (6) often formed a bicrystal with the boundary parallel to the rod. 


Table 1. 'Tin—Tetragonal Structure 99-997 %, Purity 


: (1) (2) (3) (4) (5) 
20 16 100 0 11—13 
, 10 20 70 20 8—10 
5 24 42 50 9—19 
gs 24 67 13 10—28 
0-5 24 96 4 10—35 


(1) rate of growth R (mm/min); (2) number of specimens; 
(3) percentage of single crystals; (4) percentage of crystals 
with lineage; (5) range of gradient G (°c/cm). 

- Rates of growth of 5, 1-5 and 0-5 mm/min produced the (0) type of imperfection. 
The strips of crystallite in (c) grew along the length of the rod and usually had 
orientations differing from one another by one or two degrees. In the case of 
tin, for specimens containing lineage, the angle between the rod axis and the 
¢-axis was always 70—90° and the orientations of the individual strips were 

* usually related to one another by small angles of rotation of the c-axis about 

the rod axis. This type of imperfection was greatest for the 5 mm/min rate of 

growth. Lineage formations in tin have been studied, quite recently, by 

Teghtsoonian and Chalmers (1951). 

Another type of macro-mosaic which occurred, and which was observed 
both before and after etching, consisted of fine lines about 0-1mm apart, but 
no variation in orientation was detected on examination of the surfaces containing 
_ these lines. 

On crystals grown at the rate of 20mm/min various types of rapid 
solidification phenomena, such as dendritic markings on the surface or occasional 
shrinkage cavities, were seen near the end of the specimens. 

(ii) Lead 

Tadanac lead of 99-998° purity, kindly presented by Messrs. Goodlass 
Wall and Lead Industries, Ltd., and 99-99% pure lead obtained from Messrs. 
Johnson, Matthey & Co. Ltd., were studied. No significant difference could 
be detected between the two in regard to crystallization. ‘The results 
summarized in table 2 refer only to 'Tadanac lead. 

As for tin, the effect of varying the temperature gradient had negligible 
effect, but it is just possible that higher gradients might influence the 
crystallization. In more than half of the 143 specimens examined one type of 
imperfection, lineage, was observed. ‘The strips constituting the lineage were 
about 1mm wide but they were not uniformly parallel to the length of the 
specimen. ‘The orientations of the various strips were related to one another 
by rotations of the lattice around the specimen axis irrespective of the actual 
orientation of the crystal. In columns 4 and 5 of table 2 the percentages of 
specimens with lineage, where the variation in orientation is greater than 1° or 5° 
respectively, are given. 
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It can be seen that the greatest proportion of crystals with lineage occurred 
for rates of growth of 1-5 and 5mm/min and for these speeds the number of 


crystals with a variation in orientation of more than 5° was also a maximum. 
Only for 31 mm/min were rods crystallized free from lineage, although, even 
for these, variations in orientation of less than 4° could be detected in some cases. 


Table 2. Tadanac Lead—Face-Centred Cubic Structure 99-998 % Purity 
(1) (2) (3) (4) (5) (6) 


31 13 100 0 0 11 
16°5 36 33 67 8 . 6—21 
593 il) 20 UP 53 23—50 
15 33 18 79 64 5—65 
0°5 6 50 50 50 17 
Dep 40 62 38 17 5—50 
(1) rate of growth R (mm/min); (2) number of specimens; 


(3) percentage of single crystals; (4) percentage of crystals with lineage 
>1°; (5) percentage of crystals with lineage >5°; (6) range of 
gradient G (°c/cm). 


With rates of growth up to 16mm/min almost every specimen consisted of 
20 cm or more of one predominant orientation, but in other specimens inclusions 


of crystals of other orientations appeared along the length of the rod. For 


31 mm/min there was a tendency to form a series of single crystals of 5—15cm 
in length with the crystal boundaries across the rod. For 16 or 31mm/min 
rapid solidification patterns were formed similar to those for tin. 

An average was taken of the number of seed crystals formed in the first 
millimetre or so of the specimens and, for rates of less than 5 mm/min, it 
would appear that more than two crystallites could have acted as seeds and, 
for higher rates, less than two. 

To observe the effect of oxidation several crystals were left for 24 hours in 
air before etching and coloured oxide films formed on some of the crystal faces. 
Observations on 15 specimens treated in this way showed no oxidation on the 
(100) planes. This is in disagreement with Gwathmey et al. (1945) who found 
these planes to be active, though at a higher temperature 150—250° c. 


(iti) Zinc 

Crystals were grown from 99-997°% pure zinc kindly presented by the 
National Smelting Co. Ltd. and also from zinc of 99-99% and 99-95% purity. 
Table 3 contains results of the experiments on the purest metal. Results for the 
less pure metal were similar but with a definite inclination for the formation of 
slightly less successful single crystals. For example, for the 5 mm/min rate, 
the crystals in some cases consisted entirely of strips and, for 1-5 mm/min, 
bicrystals were frequent. 

The size and form of the crystals changed rapidly with the rate of growth but 
changing the temperature gradient had little effect. ‘The general picture is 
shown in table 3, but in particular the following are worthy of note: 

(a) 16:5mm/min. ‘The beginning of each rod always consisted of polycrystals 
and then strips of crystal about 1mm wide which extended along the length of 
the specimen. ‘The orientations of these strips were not simply related to one 
another as for lineage in tin and lead. 
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(b) 5-3 mm/min. 'The rods in this case differed from those in (a) in that the 
last 10cm or so consisted either of a bicrystal with the boundary parallel to the 
specimen, or of a single crystal. 

(c) 1:5mm/min. Large strips of crystals approximately 1mm wide occurred 
at the beginning of the rod and then for over 80% of the specimens the remainder 
of the rod was a single crystal, the others being bicrystals. 

(d) 0-5 mm/min. 'The beginning of the rod consisted of one or two crystals 
and for 90% of the specimens the remainder of the rod was monocrystalline. 
Inclusions within this length seldom occurred but the crystal often lost its 
_ original uniform cross section with the production of a neck or flat surface, or 
_ by filling the mould completely over a short distance. 


Table 3. Zinc—Hexagonal Structure 99-997 % Purity 


(1) (2) (3) (4) (5) (6) 
16:5 20 0 Strips 70—90 8—30 
ses! 27. ZL. Bicrystals 55—90 2—40 
1:5 26 85 Bicrystals 12—89 3—45 
& Single Average—76 
0:5 30 90 Single 12—90 8—30 


Crystals Average—60 
(1) rate of growth R (mm/min); (2) number of specimens; 
(3) percentage of single crystals; (4) crystal form; (5) %°; (6) range of 
gradient G (°c/cm). 


For examination of the crystals thermal etching was sufficient in most cases. 
Thermal etching produced a surface with ridges parallel to the basal plane. 
Crystals with small angles between the hexagonal axis and the rod axis, and which 
were subjected to prolonged thermal etching, exhibited flat faces with ridges. 
Hexagonal prism planes of the first and second type were developed and the 
angle of 30° between them could be measured quite accurately. 

(iv) Bismuth 

The purest bismuth available in sufficient quantity was of 99-97% purity 
and was supplied by Messrs. Mining and Chemical Products Ltd. The crystals 
produced were less perfect than those of the other three metals. In fact, no 
particular rate of growth or temperature gradient resulted in a good yield, although 
rates ranging from 0-06 to 21 mm/min and gradients in the liquid metal from 
8 to 32°c/cm were tried. 

As little effect was expected from a change of the temperature gradient only 
the highest and lowest gradients available were used and these produced no change. 
The considerably wide range of rates covered is indicated in table 4. 

Change of the rate of crystallization had a marked effect on the crystal form. 
This is considered under the two following headings: 

(a) 1:5-0-2mm/min. These rates resulted in rods with, usually, one or two 
main crystals and with small inserts of varying crystal orientation. Of the 
33 specimens examined only four had lineage with orientations varying by less 
than 1°; the others had lineage similar to that in tin and lead. This lineage was 
generally less continuous than that in tin and of a more jagged appearance but 
it was invariably along the specimen. 

(b) 5-3-10-6 mm/min. Rods were produced consisting of one or two crystals, 
but broken up transversely, sometimes by short lengths of crystals of different 


a 
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orientations but frequently by ‘block’ structure. These ‘blocks’ consisted of : 
units, a few millimetres or less in length, across the rod and with near orientations 
so that the same crystal orientation extended approximately along the whole 
length of the rod. Seventeen of the crystals grown at the above rates showed _ 
block type imperfection, the others being single crystals. No lineage was visible 
in these specimens. Except for one example at 1-5 mm/min the block structure _ 
was confined to crystals grown at rates of 5-3-10-6mm/min. 


Table 4. Bismuth—Rhombohedral Structure 99-97% purity 


(1) (2) (3) (4) (5) (6) (7) 
10-6 112 3 Short Lengths 71-7 31 15 
of Crystal 
or ‘ Blocks’ 
a8) 12 + * Blocks ’ 71-6 2, 20—30 
or ‘ Inserts ’ 
5 12 2 Lineage, 71°8 3 15—40 
; some Inserts 
0-5 — 12 2 Lineage 78°7 0 20—50 
and Inserts 
0:2 9 0 Lineage 77-4 0 40 


and Inserts 


(1) rate of growth R (mm/min); (2) number of specimens; (3) number of single 
crystals; (4) crystal form; (5) average %°; (6) maximum supercooling (°c); Zz 
(7) range of gradient G (°c/em). 


§3:°DHE “CRYSTAL ORIEN LADIONS 

A specially designed optical goniometer was used to measure the orientations 

of the various crystals grown. ‘The orientation could be determined to within 1°. 

The variation in orientation along the rod could be observed directly by moving 

the rod vertically without rotation. Values of the orientation %, i.e. the angle 

between the axis of the rod and the significant crystallographic axis of the metal 
discussed, are listed in the appropriate table or mentioned below. 

(i) Tin 
The angle between the rod axis and both the [001] and [010] axes were 
measured. The results are shown in the stereographic projections (figs.-1 and 2). 
to} [Ho] 
- 5S MM/MIN : 


x 1.5 MM/MIN 
© O05 MM/MIN 


* 20 MM/MIN 
xX 10 MM/MIN 


x 


oo} ToI9) (Sol [o1g) 

Liiyee, ole Stereographic projection of the Fig. 2, Stereographic projection of the rod 
rod axes of tin crystals grown at fast axes of tin crystals grown at medium 
rates. and slow rates. 


For the crystals grown at rates of 10 and 20mm/min the distribution of the 
directions of the rod axis relative to the crystallographic axes is approximately 
uniform but with some preference for the higher values of %. For the lower 
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speeds, however, the angle % (relative to the [001] axis) is always between 
68° and 90° and there is some preference also for the [110] axis to be near to the 
rod axis. ‘The above values refer to the single crystal portions of the rods. The 
crystallites in the first few centimetres of the rods had values of 5 usually smaller, 
but occasionally larger. 
: (ii) Lead 

_ The readings taken on the optical goniometer for the various single crystals 
of lead were plotted on a Wulff net and the angles between the nearest [100], 
[110] and [111] axes and the rod axes were determined. A stereographic plot 
of the directions of the rod axes did not show any particular preferred orientation. 
The distribution of the orientations was examined statistically, there being 
values available for more than one hundred specimens. A description of the 

method used is given elsewhere, and the term ‘uniform distribution’ found 
in this way is defined more explicitly by Goss (1951). The examination showed 
that the distributions for the crystallographic axes were uniform. This applied 
to both the pure and less pure lead. | 
(iii) Zinc 
The main imperfection has been indicated previously as strips of crystal, 

about 1mm wide, along the length of the specimen. These strips are not of 
near orientations as for the lineage specimens of the other metals but consist of 
_ crystals with the hexagonal axis approximately at right angles to the rod axis 
and in all directions around the rod axis. 

y% for the c-axis was between 75° and 90° for all the crystals grown at 
5-3 mm/min except one for which % was 55°. For the fastest rate all the crystals 
consisted of strips with high values of 4. ‘The tendency for high values of % was 
still very marked even for the rate of 1-5mm/min where 78% of the crystals 
and the bicrystals had ’s above 75° (table 3). 

For the slowest rate, 0-5 mm/min, a large scatter of values of 4 occurred, but 
_ there was a slight preference for both the high and low orientations. Low 
orientation crystals have perfect cleavage planes but the crystals were not 
usually examined by cleaving. Lineage was.observed in a few cases but any 
lineage at the higher rates was masked by the production of strips instead of 
single crystals. F 

(iv) Bismuth 

The bismuth crystals, both the single and the more usual imperfect crystals, 
grown at all speeds, favoured one particular orientation. 5 for the [111] axis, 
the axis of trigonal symmetry, was 60°-90° for 85% of the specimens. In 
column 5 of table 4 the average values of y% obtained are given. For a uniform 
distribution + should have an average of 60°, but the average for the rates of 
0-5 and 0-2 mm/min is 78°, and for the higher rates 72°. An extremely small, 
but perhaps significant, effect due to the temperature gradient appears to occur, 
for the average value of % for the crystals grown under a low gradient was at 
least 5° less than the corresponding value for the crystals grown under a high 
gradient in each of the rate groups. 

Within the first centimetre or so of the rod a crystal with a low value of 
sometimes occurred, but the other crystallographic axes showed no preference 
for any particular direction. 
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§4. THE TEMPERATURE MEASUREMENT IN THE METAL gf 

Chromel-alumel thermocouples were used for the measurement of the 
temperature, a silver-soldered couple being immersed in the liquid metal of a 
dummy specimen identical with the crystal specimens. The e.m.f. was | 
recorded on a Kent potentiometric recorder the scale of which was 250 mm wide 


measuring to 6my. This maximum corresponds to about 150°c. Figure 3 is 
a reproduction of a typical trace. It represents the variation with time of the 


temperature of the mid-point of a rod and, for approximately uniform conditions | 


over the centre section of the rod, it gives, in addition, the temperature gradient 
at any instant along the rod. The temperature gradients were calculated from 
the traces. The traces were confirmed by independent measurements using a 
direct reading potentiometer. Further, in order to assess the direction of heat 
exchange between the specimens and their surroundings measurements were 
made with differential thermocouples and pairs of thermocouples. ‘The results 
are of considerable interest for the discussion of the heat flow during 
crystallization. 

The value for the temperature gradient calculated from the curves and 
quoted in the tables is the average of the gradient in the liquid metal immediately 
above the melting point and of the gradient in the solid just below the melting 
point. This value is not an indication of the shape of the cooling curve. 

(i) Tin 

Change of the rate of growth produced a distinctive change in the 
corresponding traces. Three types of trace were obtained: (i) For the rate of 
10 mm/min or more, the temperature at solidification was steady at 232°C, in 
some cases for a period of minutes, and then fell sharply. The solid gradient 
was greater than the liquid gradient (fig. 3). (ii) For the rate of 5 mm/min the 
cooling curve was smooth with no appreciable kink. (iii) For low rates there 
was a small kink at the melting point and the solid gradient was less than the 
liquid gradient in contrast to (i) (fig. 3). 
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Fig. 3. Tin: temperature-time traces. Fig. 4. Bismuth: temperature-time 


traces. 


The form of the cooling curve is thus determined by the rate of growth of 


the crystal, i.e. the speed of travel of the furnace, and the change in the conditions — 


under which the crystals were grown was reflected in the form and orientation of 
the crystals, as given in table 1. 

The equivalence of the rate of growth and the speed of travel of the furnace 
was checked by using the Chalmers (1949) ripple method. 
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(ii) Lead 

The temperature records were similar to those for tin. For the slow rates 
of 0-2, 0-5 and 1-5 mm/min there was often a decrease in the temperature gradient 
as the metal cooled through the temperature of 327°c (the melting point), the 
effect being more marked at the higher gradients. Unlike tin there was no small 
kink in the trace at low rates. For rates of 16:5 and 31 mm/min the temperature 
Was approximately constant at or below the solidification point, with the 
gradient after solidification greater than that immediately before solidification. 
This behaviour is similar to that of tin. 

The form of the cooling curve corresponded directly with the movement of 
the furnace and confirmation of the results was obtained by using thermocouples 
connected to a direct reading potentiometer and by using ripples. 

(iii) Zinc 

The chromel—alumel thermocouples were attacked by the zinc vapour, but 
by using a fresh thermocouple for each run the temperature in the solidifying 
metal was recorded. At the same time comparison readings were taken on a 
similar thermocouple shielded from the vapour but still within the furnace tube. 
For the actual crystallization runs, records were taken of the shielded 
thermocouple only and the temperature of the solidifying metal was deduced. 

For the highest rate of 16-5 mm/min there was an increase in the temperature 
gradient after solidification. A slight increase in the gradient occasionally 
occurred at lower rates, but no decrease in gradient was ever observed for low 
rates of crystallization. ‘This is in contrast to the results obtained with tin 
and lead. 

At low rates and gradients the temperatures of the shielded thermocouple 
and of the metal corresponded closely, but at high rates and gradients the 
temperature gradient in the furnace tube was considerably greater than that in 
the metal. In an extreme case there was as much as 100° difference in temperature 
_ between the metal and its surroundings. 

The gradient values in table 3 are those deduced for the metal. 

(iv) Bismuth 

Although bismuth slowly attacked the thermocouple wires, measurement of 
the temperature in the solidifying metal was still possible and records were 
obtained in the same way as for tin and lead. 

Supercooling was observed in almost all the records for rates of 5-3 mm/min 
and higher, but at slow rates only one instance of supercooling was detected. 
However, in several cases, a change in gradient for these slow rates took place 
at 10°-20°c below the solidification temperature. 

Bismuth has a much smaller thermal conductivity in the solid than in the 
liquid (Kyo4i4/Kiiquia= 1/2). Hence the slope of the temperature-time trace 
increases sharply as the temperature drops through the solidification point. 
This occurred regularly except for the lowest rates of 0-2 and 0-06mm/min. 
No change in gradient was detected in those cases. ‘These effects are shown in 


ne. +. 
§5. CONCLUSION 


The metals tin and lead were similar in their crystallization. Good crystals 
were produced for both fast and slow rates of growth but, for medium rates, 
the crystals were less satisfactory and had the greatest amount of lineage. 
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For bismuth, all rates produced single crystals though many had slight — 


imperfections. Fast rates tended to break up the crystals into blocks, low rates 


to produce lineage. For this metal, at rates of 5 or more mm/min, supercooling | 


This 9 


occurred, and this could be seen clearly on the temperature-time traces. 
supercooling was peculiar to bismuth. 

Zinc behaved quite differently from the other three metals. 
rates were single crystals produced. At faster rates strips of crystal of similar 
orientations appeared. 

In no case was the magnitude of the temperature gradient an influence on the 
growth, but the form of the temperature-time trace was dependent (in zinc only 
to a small extent) on the rate of growth. 

The rate has, in addition, an effect on the orientation. 
tin and zinc and to a slight extent also for bismuth, but the influence varies from 
metal to metal. The effects on the various metals together with a summary of 
the other results are given in table 5. 


Table 5. Summary 
Rate 
(mm/min) Tin Lead Zinc Bismuth 
Crystal Single crystal Single crystal No single Single crystals 
form crystals with 
imperfection 
5 Orientation Random Random 7 tending ws tending to 90° 
to 90° 
Temp. grad.* Zeroatm.p.: Zeroatm.p.: Increasein Supercooling : 
increase in increase in solid increase in solid 
liquid liquid 
Crystal Lineage Single crystal Bicrystal Single crystals 
form with lineage with 
d imperfection 
5 Orientation 70°<%<90° Random 55°<b<90° x tending to 90° 
Temp. grad.* Littlechange Littlechange Nochange  Supercooling : 
observed Increase in solid 
Crystal Single crystal Singlecrystal Single Single crystals 
form crystal with lineage 
— Orientation 70°<%<90° Random Random is tending to 90° 
‘Temp. grad.* Decrease in Decrease in No decrease Increase in solid 
solid solid in solid 
observed 


* At solidification. 


Though the matter is not discussed in this paper it can be shown that the 


information given by the temperature-time curves for tin, lead and bismuth 


may be related theoretically to the heat flow from the liquid to the solid by taking 
into account the release of latent heat, and a value for the ratio of the thermal 
conductivities of the solid and liquid metals may then be deduced. 

The importance of the influence of the purity of the metal on the successful 
growth of single crystals should not be overlooked, and the results of a further 
series of experiments in which the effects of added impurities to tin have been 


a ee 


Only at slow © | 


This is the case for — 


~ 
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_ investigated will be reported elsewhere. 'The effect observed with the less pure 
‘ zinc mentioned earlier appears to fall into line with these subsequent observations. 
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Mr. K. E. Puttick (British Iron and Steel Research Association). A little while ago 
when growing single crystals and bicrystals of tin I made some observations, which though 
they were made during the course of another investigation and not followed up in any way, I 


_ think are worth mentioning here. 


The single crystals were grown in Pyrex tubes and in moulds of graphite and sindanyo ; 
the latter were split cylinders about 1-5 in diameter with two holes for the metal, each about 
0-8 cm diameter and 14cm long. The crystals grown in Pyrex were of comparable size. 
‘The bicrystals were grown by Chalmers’ method (1940) in boats of graphite or sindanyo. 
‘The speed of furnace travel was about 10 cm/h and the average temperature gradient, 
though not measured, was of the order of 10 deg/ cm. 

It was found that there were differences between crystals grown in these different 


materials. First, the lineage or macro-mosaic structure was found to be very marked in 


those crystals grown in glass, with respect both to the relatively wide range of angles in the 
crystallites and to their small size. In crystals grown in graphite and sindanyo, this 
substructure, though usually observable, appeared to have wider bands and smaller dispersion 
of orientation. Figure 1 shows a Laue photograph of a crystal grown in graphite which 


_ shows no evidence of macro-mosaic: in fig. 2, grown in sindanyo there is only slight splitting 


of spots. Figure 3 and 4 of crystals grown in glass, show that the crystallites were small 
relative to the diameter of the x-ray beam and there is a range of angles of order 10°. In 
fig. 3 the disorientation of the structure seems to correspond to rotation about a certain 
direction. "This was frequently found to be the case, the crystallites forming a fibrous 
structure with the specimen direction as the common axis. Figure 4, however, shows a 
substructure which cannot be so described. 

The orientation of different crystals was observed in all cases to have a most probable 
value for the c-axis of about 90° relative to the specimen axis, but the dispersion of orientations 
about this value was much greater for specimens grown in graphite than for those grown in 
glass. The orientation of the a-axis showed no tendency to a particular value. (These 
crystals showed no sign of polycrystals at the end of the rod as did some of those reported by 


- Goss and Weintroub.) 


In growing bicrystals from two single crystal seeds, symmetrically disposed with respect 
to the axis of the specimen, it was found that with graphite moulds the boundary, while 
being under these conditions roughly parallel to the specimen axis, was irregular. In 
sindanyo moulds, however, the boundary was microscopically straight. Straight boundaries 
confined to the median plane of the specimen were only consistently obtained when the 


_ ¢-axes of the seed crystals were perpendicular to the specimen axis. 


It seems likely, then, that the thermal properties of the mould have a considerable 
influence on the growth of crystals, presumably because they affect the temperature gradient 
and rate of growth. Teghtsoonian and Chalmers (1951) have confirmed that the rate of 
growth does affect the size of the macro-mosaic units. It would be particularly interesting to 
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know whether within the limits imposed by other considerations, the rate of cooling at 
nucleation does exercise a significant effect on the orientation of crystals, relative to the 
specimen axis. One can see that such an effect should exist, if we assume that the 
probability of a particular orientation is a function of a characteristic time of formation ofa 
stable nucleus. Altering the rate of cooling implies varying the period of time spent in the __ 
neighbourhood of the freezing point, and so affecting the probability of formation of a — 
nucleus of any particular orientation. sz 

Suppose, for example, that the orientation « of a crystallographic axis relative to the | 
specimen axis for a number of crystals has a normal probability distribution about a mean — 
Xm, With a standard deviation o. The probability P(«) of a particular orientation is then 
inversely proportional to the probable time of formation of a nucleus of that orientation, and 
we suppose that nucleation occurs when this is comparable with the time spent in the 
temperature range (0, 0+66) 

ie. Pia) {1/(dt/d0)}80 or o 12a) -!?exp[— (a—dm)*/ 207] = b(d0/dt) 

where } is a constant. Putting ~=a,, we see that the standard deviation of orientation is. 
inversely proportional to the rate of cooling at nucleation. The question we have to 
answer, then, is: why is there a minimum nucleation time characteristic of a certain 
orientation ? “One possible reason is that the surface energy is anisotropic, so that the 
freezing point and rate of growth are functions of orientation (Chalmers 1949). 

This picture, however, offers no explanation of the curious fact that the conditions which ~ 
apparently reduce the range of orientations between different crystals increase the range of 
the substructure orientation. 


ee 
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Asymmetry of Texture in Cubic Crystals Grown from the Melt 
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ABSTRACT. Careful electro-polishing of single crystals pheres of copper (face-centred 
cubic system) grown from the melt has revealed that they appear to consist of a number of 
parallel crystallites of square cross section with their lengths parallel to the [001] axis which 
lies nearest in direction to the ‘ direction of growth’ in the melt. 


INGLE crystal copper spheres about 15 mm in diameter on shanks (6-28 mm 
diameter) were carefully turned from cylindrical ingots which had been 
formed by solidification of the melt from the bottom upwards. They were 

then anodically polished to remove the superficial layer of polycrystalline material 
left by turning. This was done very carefully in a bath of 50% v/v H,PO, and 
at current densities of about 1 amp/cm*, holding the shank slantwise, rotating the 
crystal in the bath and dipping it in and out continually so that scoring due to the 
movement of bubbles upwards or the downward movement of dense electrolyte 
near the anode (‘tear’ formation) was minimized. This treatment aimed at 
producing as uniform a surface as possible. After about 15 minutes, when all 
polycrystalline material had been removed, dimples or pits were observed on 
the two cube faces which had been uppermost and lowermost respectively during 
growth of the crystal from the melt. The dimples lay along lines parallel to the 
cube face sides of the generating cube, forming a sort of cross-hatched pattern. 
On the other four cube faces, which lay nearer the vertical during growth of the 
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crystal, the dimples lay along lines parallel to that cube edge direction which 
was nearest the vertical. 

This is illustrated in figs. 1 and 2. The line AB in fig. 1 (a) is the direction 
of steepest temperature gradient and the line CD that of the nearest cube edge 
direction to it. Figure 1(d) indicates the lines along which the dimples would 
lie on a cubic crystal cut from the ingot. Figure 2 shows the dimples on a 
spherical crystal. The line CD is the same in all three illustrations. 

The dimples were about 0-1 mm diameter with about 0-3 mm between rows. 
The rows were not perfectly straight on the uppermost and lowermost faces 
_ but waved slightly rather like the markings in a finger print. Cracks could be 
observed on the sides of the crystal on the {011} (cube edge) regions. They were 


® about 0-45 mm apart and appeared to form the boundaries of crystallites. In 


general, the dimples tended to avoid the cracks rather than to coincide with 
them. The only conclusion that can be drawn from these observations is that 
the crystals grow along the cube axis nearest to the direction of steepest 
temperature gradient, rather than along the temperature gradient itself. It will 
be interesting to note whether similar non-equivalence along main crystal axes 
is observed in the mechanical, magnetic and electrical properties of metal crystals. 


Fig. 1. Fig. 2. 


It is generally supposed that crystal growth occurs at dislocations and it is 
therefore to be expected that the direction of growth would be a singularity 
since these could hardly be randomly disposed. ‘This may explain why the 
centres of faces in crystals growing from solution are frequently the points at 
which fresh layers commence (Bunn and Emmett 1949). 

A very tentative explanation of the pattern is that the crystals grow as a 
series of closely aligned crystallites, each exposing cube facets and extending 
by deposition of material on the uppermost of these. ‘This raises interesting 
speculations as to the initiation of the crystallites at their lower ends which must 
lie on the periphery of the ingot (see fig. 1(a)). On the other hand, dendritic 
growth may give rise to the effect. 

The author is indebted to Mr. A. J. Goss for drawing his attention to the 
fact that a similar structure, but on a rather finer scale, had been observed in the 
grains of polycrystalline copper-rich alloys by Northcott and Thomas (1939), 
using repeated mechanical polishing and etching. ‘These authors assumed that 
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the effect was due to dendritic growth but in that case it would have been expected 
that considerable intermeshing of the dendrites would have taken place at the 
boundaries between crystals. Most of their photographs show little indication 
of this. 

Similar patterns have been observed by the author on the upper menisci of 
nearly all copper single crystals formed by upward growth from the melt. This 


is evidence that the effect is due to facet formation during growth, contrary to | 


HN a Pett ides Mes oe arent = Cael 
if Lb i ~~. 


the result found by Anderson for zinc (Anderson 1932). Since our crystals — 


were formed by very slow growth, the crystal—liquid interface may be expected 
to display the equilibrium form for the crystal, whereas Anderson allowed more 


rapid growth, and it has been shown by the author (unpublished) that this would — 


tend to suppress facet formation. 

The copper used was spectroscopic standard material* but, although grown 
in purified hydrogen, may have acquired a trace of oxygen, sufficient to give rise 
to an adsorption complex on the slowly growing crystal and hence to aid facet 
formation. Occasional occlusion of oxygen during this process may have given 


rise to the regions which become pitted on polishing. Alternatively, the pits may ~ 


have centred on dislocations round which the growth of each crystallite took 
place. 


The author wishes to thank Miss S. Treadgold and Mr. C. Engel for the 
illustration and photograph. 
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ABSTRACT. Crystals of thallium-bromo-iodide, or KRS-5, are transparent to radiation 
of wavelengths between 0:5 yw and about 40 w. This note briefly describes the growing 
of KRS-5 single crystals and their optical working. The crystals are grown in vertical 


cylindrical crucibles by a progressive freezing of the melt from bottom to top. The physical, — 


chemical and mechanical properties of the material are summarized. 


per cent form a mixed crystal, known as KRS-5, by solid solution (Koops 


[erent bromide and thallium iodide in the proportion 44:56 mole 


1948) which is transparent to radiations in the region 0-5u—-40. Since 


KRS-5 can be formed into optical components by grinding and polishing and 
is chemically stable, it is a very valuable addition to the range of optical materials. 


Although KRS-5 is cubic in structure and therefore isotropic it is for various © 


reasons desirable to have it in the form of single crystals; for instance, inter- 

crystalline cracks are liable to form along the boundaries when the material is 

being ground or polished. ‘The most suitable method of growing the single 

crystals appears to be that devised by Bridgman and applied to KRS-5 by 
* Sunplied by Messrs. Johnson, Matthey, Ltd., London. 
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Smakula and others (Elliott and Mayes 1945, Tuttle and Egli 1946). 
In this method the highly purified bromide and iodide, mixed in the proper 
_ proportions, is contained in a cylindrical crucible of platinum or Pyrex, of which 
_ the closed end is conical. The steady and progressive freezing through the 
_ melt which is necessary for the growth of a single crystal is encouraged by 

lowering the crucible, pointed end first, through two coaxial cylindrical 

furnaces. ‘The temperature of the upper furnace is about 460°c, that is, some 
_ 45°c above the melting point of the mixture ; the lower furnace is at roughly 330° c. 

The furnace temperature must of course be closely controlled to prevent a too 
_ rapid wandering of the region of freezing. 

A steep temperature gradient in the region of freezing is an essential feature 
of the method. It has the obvious purpose of ‘anchoring’ the region of freezing 
and it is also believed to assist the rejection of impurities which might disturb 
uniform crystal growth or the optical properties of the crystal. In order to 
achieve a very steep temperature gradient we have tried lowering the crucible 
into a bath of liquid metal in the lower furnace. This also has the advantage that 
after solidification the crystal enters a region of nearly uniform temperature. 
However it seems that the temperature gradient near the freezing boundary may 
be too high and that the crystal structure, which is easily deformed at high 
temperature may be dislocated. At present, during modifications to the 
apparatus, we have returned to the use of empty air space in the lower furnace 
and the results are very promising. An annular baffle between the two furnaces 
helps to maintain a reasonably high temperature gradient. With this 
arrangement the rate of growth for a 2 in. diameter crucible is about 1 mm/hour. 

When solidification is complete the crystal must be annealed and since it is 
plastic even at room temperature and is subject to delayed elastic effects the 
whole cooling process must be carried out very slowly. 

KRS-5 has a density of 7-4g cm™ and a linear coefficient of expansion about 
seven times greater than that of glass. It has that peculiar combination of 
toughness and plasticity commonly described as ‘horny’. It cannot be cleaved 

_ but it can be sawn, turned on a lathe, abraded with emery and optically polished 
with the usual polishing powders, though the technique of optical working 
differs in some ways from that for glass (Cooke 1947). In particular, flow lines 
formed below the surface during grinding appear as raised lines some time after 
the surface has been polished. To avoid this the blank is annealed again after 
grinding (Smakula and Klein 1950). 

The refractive index of KRS-5 in the infra-red is 2-45—2-25; the dispersion 
falls steeply from 0-1 at 1 to 0-002 at 10u, and above that wavelength it remains 
low (Plyler 1947, Tilton et al. 1949). 

The solubility of KRS-5 in water is comparatively low (0-05 g/100 cm®) and 

polished surfaces do not deteriorate in a humid atmosphere. 
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ABSTRACT. ~ Following the development of methods for growing single crystals of 
materials having good optical properties, the search has been intensified for new substances. 
One such is caesium bromide, and this has been grown by the Kyropoulos—Pohl method 
of drawing from the melt; a modification of the standard process is described. 

In the search for new materials, it is useful to have some way of predicting the useful 
long-wavelength transmission limit, and a rule is stated which forms a useful guide. 


Si INTRODUCTION 


RYSTALS such as the alkali halides which are used as optical materials are 
( for the most part grown by one of two methods. In one, which owes 

much to the work of Bridgman and of Stockbarger, the powdered material 
is moved slowly across a small finite temperature discontinuity maintained by 
two adjacent furnaces held at differcnt temperatures. 

In another process, developed by Kyropoulos and extended by Pohl and others, 
the crystal is drawn from the melt. 

Some crystals are best grown by one of these two methods. In the writer’s 
laboratory the second method has been used extensively for growing large crystals 
(9 in. diameter) of rocksalt and potassium bromide, and smaller ones of sodium 
iodide with added thallium, for scintillation counters. Occasionally small ones 
of sodium bromide, potassium iodide, silver chloride and rubidium bromide have 
been grown. In the growing of these materials a modification of the standard 
method has been developed, and this is described. 

In the latter part of the paper is stated a rule which has been found useful for 
predicting whether the infra-red transmission of a material is likely to be favourable 
and therefore whether it is worth while to attempt to grow it. 


§2. GROWING FROM THE MELT 


A seed crystal is selected and mounted in a chuck, preferably of nickel or of 
nickel-plated material, which is water-cooled. ‘The powdered material is placed 
in a crucible (often of unglazed porcelain) standing in a furnace. ‘This furnace — 
has a heating element at the bottom as well as the sides, to prevent the formation 
of a stalagmite. ‘The seed crystal in its chuck is suspended over the top of the 
powdered material with the water running, and the heat is switched on. The 
material slowly melts, and for a 9 in. crystal this part of the process takes a week— 
hurrying leads to cracked crucibles and spoiled furnaces. 
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The seed is then lowered into the melt which is raised about 20° above the 
melting point, and the surface is melted off to form a new surface. The temper- 
ature is then. lowered to a few degrees above the melting point and growth 
commences. <A ‘bun-shaped’ mass forms on the seed, and extends sideways until 
it would seem to be in danger of reaching the wall of the crucible. At this stage, 
orthodox procedure is to raise the seed and its attached bun, and form a new bun on 
the underside. If this is done the finished crystal has a corrugated surface. 
It has been found in the writer’s laboratory that by paying special attention to the 
heating arrangements, especially avoiding any asymmetrical cooling, and by 
proper selection of the size of the vessel, it is not necessary to lift the crystal, but 
one can rely on the liquid level dropping. ‘This comes about because the solid is 
denser than the liquid. A smooth-faced crystal results, instead of the corrugated 
face normally obtained. This growing process also takes a week. Finally, 
another week is spent in the crystal cooling down. 
While from first to last three weeks is necessary for a 9 in. crystal, a 3in. one 
need only take about two days. 
A special advantage of drawing from the melt is that one may observe the crystal 
from time to time to see if it is growing as a single or a poly-crystal, and, in the latter 
ase, one may melt back and start again without great loss of time. 


§3. CAESIUM BROMIDE 


Powdered caesium bromide of high purity was obtained from Dr. Norman 
Wright, of the Dow Chemical Company, and was grown into a crystal by the 
method just described without any particular difficulty. A seed was made by 
melting some and allowing it to cool, and then picking out a sufficiently large piece 
from the cooled mass. 

From the crystal was cut a sixty degree prism of good quality, and this is being 
used to measure the refractive indices of the material. 

This material is likely to be important as a rival to thallium bromo-iodide as an 
infra-red material; it is likely to be usable to a wavelength in excess of 40 microns. 
Recent measurement by Plyler gives a transmittance of 76° at 38 for a 

~7mm-thick window. Unlike the thallium compound, it is colourless and it has 
mechanical properties which, from the optical worker’s point of view, are little 
worse than those possessed by potassium bromide. It is hygroscopic but not 
extremely so. The crystal was grown by Mr. J. Skinner. 


§4. INFRA-RED TRANSMISSION LIMITS 


The search for useful infra-red transmitting materials continues, and it is 
convenient to have some guide as to what one may expect beforehand. In com- 
piling a list of wavelength limits of useful transmissions—by which is meant here 
the wavelength beyond which a spectroscopist would not seek normally to use a 
material—the author noticed a rule which so far appears to be obeyed by these 
materials. The upper wavelength limit of useful transmission is about one-third 
of the residual ray (reststrahlen) wavelength. ‘This limit corresponds well with 
the wavelength at which the transmission of a two centimetre thick sample has 
dropped to about 10%. 
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The extent to which the rule is followed may be gauged from the following ~ 
table : 
Transmission Limits and Residual Rays 


Material Useful wavelength limit (2) Residual rays 
Fused silica — 3-6 9 
LiF 6 17 
CaF, 9-5 32 
NaCl 7 5D 
KCl a 63 
AgCl PIRES 82 
KBr Dy 82 
KI 32 94 
TIBrI 38 120 
[CsBr >38 < 134] 


The residual ray wavelength fits better than the infra-red absorption maximum, 
as may be seen, for example, by considering lithium fluoride. The residual ray 
wavelength for caesium bromide is so far not known, but the absorption maximum 
is 134, and the residual ray maximum will be less than this, by an amount due 
to the Forsterling correction. One may conclude that caesium bromide should 
have a useful limit in the above sense of just over 40 y, and the likelihood of this 
being so seems to be indicated by the measurement of Plyler referred to above. 


Crystal Growth Ideas and Crystallographers* 


The writer was brought up to include, in the latter term of the title, anyone who worked 
on the enormous variety of problems connected with the properties of crystals: optics, 
hardness, glide and slip, chemical failure, solubility of discrete crystal planes, crystal structures 
as evaluated by x-rays and similar methods, and the many problems and difficulties asso- 
ciated with their growth. It has not escaped the attention of the writer that a ‘ crystallo- 
grapher ’, at least in England, need be conversant with only one of these categories and can 
be quite ignorant of anything else happening around him, indeed, need not be very conver- 
sant with the elements of descriptive morphology known for over a hundred years. How 
else would it be necessary to refute the current implications of a theory, possibly valid 
without these latest frills, wherein a phenomenon occurring at times on a face which scarcely 
grows, is supposed to supply the ‘ motif’ of growth for all the crystal. The side facets, on 
which most growth occurs, are dismissed as ‘ kinked’ or ‘ incomplete’ faces, though 
{1010} of either CdI, or SiC is quite as complete as {0001}, it being forgotten that a unit-cell, 
e.g. of a hexagonal crystal, must be clothed in {1010} as well as {0001}. What we require, 
then, is a theory which makes them grow more rapidly, and solid phase to separate via them, 
as against those on which spirals appear, and on which, from the very shape of the completed 
crystal, anyone versed in elementary morphology would immediately recognize that little 
or no addition to the crystal had taken place during growth. 

Briefly, in addition to the beautiful observations of Verma (1951) and Amelincx (1951) 
on SiC and further ones on Cdl, by Forty (1951), on mica by Amelinex (1952), on colloidal 
gold by Forty (1952) and AIB, from Al (Horn et al. 1952), the writer has observed all the 
phenomena on a scale up to nearly a millimetre between loops, both closed and straight- 
edged, all practically visible to the naked eye but not commented upon by others. He 
appeals to anyone to examine for himself a lump of SiC crystals—not specially selected 
ones—and verify that the incidence of reasonably good crystals is (in our present sense) 
far above one in ten: that really good specimens are far from common and that the great 
bulk of SiC crystals are dendritic, often curved-surface masses, and it would be difficult, with 
such unreliable material as this, to connect up with any theory of crystal growth. The writer 
has all along considered that some vortex of vapour, possibly containing the material for 
growth, was present at the conclusion of growth. This formed a smear of SiC at the end 


* Read at Physical Society Meeting at Southampton on 19th December 1951. 
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and imposed some kind of spiral on it. Often (fig. 1), the spiral-covered dendrite is visible. 
The crystal is punctured in depth by holes of various shapes and these are covered with a 
thin smear of SiC, the latter being detected from more solid SiC underneath through total 
reflection at the surface SiC-air, rendering the crystal paler in reflected light at these 
numerous places. The deposition of the spiral is, in this view, an event occurring near the 
end of crystallization. 

The picture of Cdl, (fig. 2) shows excessively thin sheets, 500 A thick by } mm across. 
and several thousands of times wider than they are thick. Negligible amounts enter the 
crystal by growth on such {0001} surfaces. Even in an average Cdl, crystal the sides are 
practically responsible for over 99° of the solid phase. While the drop evaporated (3 hours) 
the shade of Newton’s colours (straw, Ist, and green, 2nd order, latter bigger) persisted, 
but the triangular seed observable on the straw-coloured one fringed itself in, all this time, 
with red, Ist order, and became itself blue, 2nd, the crystal growth, all this time, on the bulk 
of {0001} was nil. In practically all cases given so far, we know nothing of how the crystal 


Figure 1. A straight-forward spiral on SiC, Figure 2. Excessively thin plates of CdI, 
visible with a hand-lens. showing Newton’s lowest orders of colour, 
formed under a cover-slip. ; 


grows; alternatively, we artificially restrict it in some way (coverslip with CdI,, innumerable 
cavities and vapour streams in SiC). Mica grows along its optimum direction which is 
along a sheet from side to side. How can a growth spiral come about under such circum-. 
stances ? ‘There is something suspicious in all other cases quoted so that bona-fide crystal 
growers drop the idea like hot cakes as inapplicable to their work. In the writer’s opinion 
the early views of Frank (1949) have something to be said for them. ‘The surface spirals 
and other similar features are also well attested as a phenomenon, so far unexplained, but 
to make the former theory depend upon the latter observations does not appear to the writer 
to make sense. 
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ABSTRACT. A brief description is given of Kirkpatrick and Baez’ x-ray reflecting micro- 
scope. The types of aberration possible in optical systems without axial symmetry are 
discussed and the effect of spherical aberration on the diffraction pattern of an object point 
is evaluated. The resolving power of a circular mirror is calculated, and a method of 
correcting field obliquity is described. A method is given for the design of a figured mirror 
to give smaller spherical aberration with zero field obliquity. 

Design equations are given for a system employing two circular mirrors in each plane, by 
which spherical aberration, obliquity and curvature of field can be corrected. The computed 
performance of the system is described. 


§1. INTRODUCTION 

T has been pointed out (Jentzsch 1929, Kirkpatrick and Baez 1948, 
| Kirkpatrick 1950) that the phenomenon of total external reflection of 

x-rays by solid surfaces makes possible the formation of real images of 
objects illuminated by x-rays. In particular, such images may be formed with 
a magnification greater than unity, thus suggesting the possibility of an x-ray 
microscope. 

The scheme proposed by Kirkpatrick and Baez is illustrated in fig. 1. A small 
object illuminated by x-rays is placed near a concave mirror so disposed that the 
central ray of the x-ray beam strikes the mirror at a glancing angle 7. If the angle 
be sufficiently small, total reflection of the x-rays takes place at the surface of the 
mirror, which exerts a focusing action on the beam. The image is formed at a 
distance from the mirror which is greater than that of the object, and under these 
conditions is larger than the object. However, the image is strongly astigmatic, 
and to overcome this defect a second reflector is placed in the x-ray beam with 
its plane perpendicular to that of the first, and so arranged that the astigmatism 
of the first mirror is cancelled by that of the second. The magnification of the 
image is then different for two mutually perpendicular directions. Kirkpatrick 
and Baez outline a method for overcoming this defect, but in any case it is not 
serious and can be corrected later in the optical enlarger. 

Apart from any geometrical aberrations which may arise, the attainable 
resolution is given by the maximum aperture angle and the x-ray wavelength. 


For large magnifications the emergent rays will be nearly parallel, and under 


these conditions it is easily seen that the maximum semi-aperture angle is equal 
to the glancing angle z. ‘The resolving limit is therefore equal to A/2z where 2X is 
the x-ray wavelength. 

‘The maximum glancing angle for which total reflection will take place from 
a given substance is, in first approximation, proportional to the wavelength, 
so that the resolving limit tends to be independent of the wavelength for a given 
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mirror material. Kirkpatrick and Baez give i=7 x 10°A (where A is in cm) for 
the glancing angle from the platinum metals which, they state, give the largest 
angles for a given wavelength of any known substance. This value gives a 
minimum resolving limit of 704, but this cannot be approached closely in 
practice because of the geometrical aberrations of the system. Prince (1950) 
has calculated the resolution on the basis of the dimensions of the geometrical 
confusion patch for a system using mirrors of circular longitudinal cross section, 
and estimates the resolving limit to be in the region of 3 400A. 

The effects of surface conditions and accuracy of the mirrors have been 
discussed by Jentzsch (1929), Ehrenberg (1947, 1949 a,b) and Lucht and 
Harker (1951). . 

The object of this paper is to examine the types of aberrations which can 
arise in a system of mirrors used at very small glancing angles, and to indicate 
means by which such aberrations may be corrected or reduced. 


§2. ABERRATIONS OF NON-AXIAL SYSTEMS 
The most noticeable feature of optical systems such as that illustrated in 
fig. 1 is that there is no axis of symmetry as in the case of more conventional 
systems. This suggests that the possible types of aberrations may be different 
from those occurring in more familiar cases. 


Pbject 


Cylindrical Mirrors 


Fig. 1. Schematic arrangement of the Kirkpatrick—Baez Fig. 2. Spherical aberration. 
microscope. 


Considering one plane at a time, the position of emergence of a ray from the 
exit pupil can be described by a single coordinate x. Rays with the same value 
of x but separated in a direction perpendicular to the plane of incidence on the 
mirror of the central ray will suffer the same aberration if the mirror be of 
cylindrical shape, and, if it be spherical, the aberration will differ only by an 
insignificant amount for the values of angular aperture considered here. Also, 
as the deflections produced by the two mirrors are at right angles to each other, 
interactions between them are negligible for small aperture and fields. ‘The 
distance of the object point from some arbitrarily chosen origin, measured at 
right angles to the imaging beam (which is assumed small in angular aperture) 
can be designated by a ‘field coordinate’ h. Then the angular aberration of the 
ray can be conveniently expressed as a series of products of powers of x and h: 


Angular aberration = @49% + dygX" + dggX? . 2. . + Ag h + Ay,xh + dg,x*h 
+ Ag R+ .... + Aggh® +ay.xh? + dopx*h? + AgoXh Ato... ceceee (1) 


In a system with axial symmetry some of the terms in (1) would be prohibited 
by symmetry considerations, but here all the terms may have finite values. ‘The 
first few terms may be interpreted as follows: 

a,,x, change of focal plane; a, x2, primary spherical aberration; dggx°, 
secondary spherical aberration; a ,h, change of magnification; a,,xh, obliquity 
of field; a,,x?h, coma; dy.h?, distortion; a,,xh?, curvature of field. 
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.The term ‘spherical aberration’ is used for dy9x?, dggx°, as these quantities © 


depend only on the aperture. The primary spherical aberration varies with 


the aperture in a manner similar to that of the tangential component of coma in — 


axially symmetrical systems, but does not vary with the field. 


§3. SPHERICAL ABERRATION 


Consider a mirror of circular cross section and radius of curvature R on ~ 


which is incident an infinitely narrow pencil of rays diverging from an object 
point O (fig. 2), distant / from the point of contact P of the principal ray with the 
mirror. The pencil is focused at a point I, distant l’ from P. ‘Then, if the 
glancing angle made by the principal ray with the mirror be 7, the distances / and I’ 
are connected by the relationship 1//+1/l'=2/iR. A ray emanating from O 
and striking the mirror at a small but finite distance r from P will intersect the 
principal ray at a point I’, distant a from I. Kirkpatrick and Baez give an 
expression for the distance a which may be put into the rather more convenient 


form aa — SHIR) _ Srl GR=S) (2) 
= Ut 


or, in terms of the aperture angle « between the aberrant ray and the principal ray, 


_ 3x0'(L-iR) 3a? GR-L) 
a=— i(21—iRP => 5 plas isnexoani ts (3) 


This gives the primary spherical aberration only, but is a fair approximation for 
small apertures. In what follows, computed results are obtained by ray tracing 
and not from this formula. It will be noted that the spherical aberration is zero 
when /=I'=7R, 1.e. when object and image both lie on the Rowland circle of 
the mirror (the circle of radius R/2, passing through the centre of curvature of 
the mirror and tangent to it at P). This is, then, a particularly favourable case for 
image formation, but as the magnification is unity it is of little interest in 
microscopy. 

The size of the geometrical aberration patch is equal to «a, but in the practical 
case the distribution of intensity is modified by diffraction effects, and it is 
necessary to take these into account. 

If the longitudinal aberration is equal to ka the deviation from circularity 
of the wave front for any angle « is equal to ka3/3. 

The intensity at a point in the image plane through I and normal to PI is 


given by 1 Tre De 
I=z5 | __ C085" (Ra /3 + Sa) ds | jot octal eee (4) 


where S is the distance of the point considered from I, and the unit of intensity 
is the value at the focus if the system were free from aberration. For k=0 
(4) gives the well-known result J=sin? (27S«/A)/(27S«/A)?. The integral has 
been evaluated numerically for values of ka?/3 equal to A/4, A/2 and 3A/4. The 
results are shown in fig. 3. 

The most obvious feature of the curves as the aberration increases is a 
progressive shift of the maximum towards the side to which the rays are deviated, 
as is to be expected. ‘There is also a fall of the maximum value of the intensity, 
a progressively smaller fraction of the total energy being concentrated in the 
central lobe and a corresponding increase of intensity in the side lobes on the 
side to which the central lobe is displaced. ‘The lobes on the other side are 


eS 
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diminished in height, but a rather unexpected phenomenon is the development 


of a long plateau or ‘toe’ to the intensity curve on this side. 
Since the width of the central lobe does not vary greatly for moderate amounts 
of aberration, it is evident that the resolution is not profoundly affected. The 


_ contrast will diminish, however, by reason of the large side lobes, and this effect 


may well be of importance if the instrument is to be used for the measurement 


of differential x-ray absorption, as, for example, in the location of specific elements 


(Engstrom 1947, 1949). In this case it is evidently desirable to keep the spherical 


_ aberration down to the lowest possible value. 


Sa/ a 


Fig. 3. Intensity distribution in focal plane for varying amounts of spherical aberration. 


For a mirror of circular cross section, and of focal length f, used with a 
glancing angle of 7 for the principal ray, the wave front distortion p for a 
semi-aperture angle « and infinite magnification is given by p=fa'/47. If a 
quarter of a wavelength of spherical aberration be allowed, the resolving limit d 
is given by d=4(fr2/i)18. nee) 
To obtain the variation with wavelength, the aperture angle is obtained by 
substituting 7 =7 x 10°A in (5), when the resolving limit is given by d=0-0056(Af)#. 
For i=0-01, A=1-4A and f=1cm, d is equal to 1350. It is evident from (5) 
that this value changes only slowly with the wavelength or the scale of the system. 


§4. OBLIQUITY OF FIELD 

Kirkpatrick and Baez point out that a circular mirror used in this manner 
gives a focal surface which is strongly inclined to the normal to the principal ray. 
However, because of the large spherical aberration, the inclination of the focal 
surface for ray pencils of small aperture is a function of the position of the aperture 
stop. For the particular case of infinite image distance, as shown in fig. 4, the 
angle PB between the focal surface and the principal ray is given for. small 
apertures by 

tanB=1i1R—2D){(iR+D), = = = — wseveee (6) 

where D is the distance, measured in the space of the short conjugate, between 
the point of incidence on the mirror and the aperture stop. 

In the case considered by Kirkpatrick and Baez, the aperture stop is the 
mirror itself, ie. D=0; the focal surface thus makes an angle 7 with the principal 
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ray. This is evidently a very unfavourable condition, for the object would have 
to lie in a plane almost coinciding with the line of sight. However, it can be seen — 
from (6) that B becomes 90° if D=—iR. This corresponds to a virtual stop on — 


the short-conjugate side, but the same effect is produced by a real stop on the | 


long-conjugate side at a distance 7R/3 from the mirror, which is the image formed 
by the mirror of the required virtual position. The system is then as shown 
in fig. 5. : | 

The focal surface will, of course, only be flat in first approximation, and 
computation shows that the sharp field is limited by field curvature. Figure 6 
gives the results of computation on a system for which R=200, 7=0-01, giving 
a focal length of 1 cm. 


i 
7. Focal \ Focal 
7 Surface \ Surface 
< [/Aperture Stop \ 
\ 

\ 

Wee \ 
> ik/3 \ 
Aperture Stop Sie 


Fig. 4. Obliquity of focal surface. Fig. 5. Correction of field obliquity by 
, aperture stop. 
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Fig. 6. Field curvature of circular mirror Fig. 7. Variation of magnification 
with aperture stop. with field. 


As the glancing angle varies over the field, the system must be used at a 
wavelength somewhat greater than the 1-4A suitable for z=0-01. ‘Taking 
A=1-6A, which is suitable for the field shown, the semi-aperture angle « for a 
quarter wavelength of aberration is 0-545 x 10-3, giving a resolution d=1 500A. 
The physical depth of focus A/2«? is then 0-0275cm. It is seen that the system 
is suitable for a field of 37-5 microns. ‘The mirror is about 0-26cm in length, - 
and the maximum length contributing to any one image point is 0-062 cm. 

This may be compared with the case of a circular mirror with the aperture 
stop at the mirror surface. ‘Taking 7=0-01, the physical depth of focus is reached 
at a total field of 5-5 microns, i.e, the sharp field is very small in extent. 
Reasonably extended fields can therefore only be obtained by employing much 
smaller aperture angles. 

The aperture angle varies somewhat over the field for a fixed size of stop, 
so the depth of focus also varies ; however, its variation is such as to give a somewhat 
larger field than the 37-5microns quoted above. It follows that the 
magnification also varies over the field, but this is not a serious defect. Its 
variation is shown in fig. 7, where the magnification is given in terms of its value 
at the point where the glancing angle has the value 0-01. 
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§5. FIGURED MIRRORS 


It is evident that spherical aberration can be removed completely if a parabolic 
mirror be used in the case of infinite image distance, but the field obliquity is 
very large. ‘The radius vector from the focus of the parabola must meet the curve 
at a very small glancing angle, i.e. the latus rectum of the parabola must be very 
small compared with the radius vector of the portion of the curve used. Under 
these conditions the parabola can be represented with sufficient accuracy over 
the portion used by the polar equation r=4r,i?/(27+«)®, where 7 is the glancing 
angle of the principal ray from the focus and « is the semi-aperture angle taken 
as a running coordinate. Incident rays parallel with the axis of the parabola 
will then come to a focus at the origin. A ray making an angle 6 with the axis 
and incident at the point for which «=0 will be reflected so as to pass the focus 
at a perpendicular distance h=Ory. A ray parallel with the first passes the focus 
at a perpendicular distance given by h+Ah=46rgi?/(2i+«)? =6r)(1—«/i) very 
nearly. Hence Ah= — Orox/t. 

This ray and that from the point « = make an angle « with each other; 
therefore the distance from their point of intersection to the plane through the 
focus and normal to the axis is Ah/x=—6r,/i. Since R=6ro, the intersection 
lies on a plane passing through the focus and making an angle 8 =2 with the axis. 
This is the focal surface for the small part of the parabola employed, and the 
obliquity of the focal surface is therefore equal to that for a circular mirror with 
the aperture stop at the mirror surface. As there is no spherical aberration in 
the case of the parabola, it is evident that the field obliquity cannot be rectified 
by a suitable placing of the aperture stop. 

This suggests that by departing somewhat from the parabolic form it might 
be possible to design a combination of mirror and stop which, while possessing 
smaller spherical aberration than the circular mirror, could be given an acceptable 
value of field obliquity. 

The method by which this can be achieved is as follows (fig. 8). The surface 
AB is represented with reference to the point F by the polar equation 


WY = 2ir-t2— f(r), C7) 


where ‘I’ is the angle between the radius vector and a fixed direction, and f(r) 
is small compared with VY’. The line FP is taken asa principal ray and FP assumed 
equal to unity. The ‘unperturbed’ curve ‘i’ =27r-"? then represents a parabola. 
The glancing angle z’ between the radius vector and the curve is given by 


, ~ i ; 
i’ >tanz Steg sat (7): DEAS Ar (8) 


The angle 0 made with the axis by a ray from F after reflection is 
6 = —21' = —f(r) —2rf'(r). 


The aperture stop is taken to be at S, distant / from F, and the focal surface is 
taken to be a plane CD through F and normal to FP. 

The principal ray for an object point F’, distant k from F, is SF’, cutting the 
curve in P’. Then F’P=FP’=r’, say. The angle made with the axis by this 
ray after reflection is 


6’ = —f(r')—2r'f'(r') + FPF’ = —f(r')—2r' f(r’) thr’. vases (9) 


,| 
i | | 
586 F. Dyson a 
q 
Another ray, F’P” from F’, meeting the curve at F’, where F’P" =r’ + Ar’, is, | 
after reflection, inclined to the axis at an angle | 
h 

z. 
0 + AW’ = —f(r' + Ar’) —2(r' + Ar’) f(r’ +49) + sag ‘ 7 


Then the angle between the rays F’P’ and F ee after reflection, is 


0’ = —f(r' + Ar’) —2(r' + Ar’) f(r’ + Ar’) + ——\s -| — f(r’) —2r'f'(r') + =| 
~ Ar’ | 3rr) + 2rP (r+ 7 A a ee (10) 


The degree of approximation is such that (10) takes account of variations of A@’ 
which are proportional to Ar’ but neglects those due to the square and higher — 
powers of Av’, i.e. (10) gives the extent to which the focal surface conjugate to 
a plane at infinity deviates from CD but does not take account of spherical 
aberration. -From (10) the condition for CD to be the true focal surface for 
parallel rays is 


Tr — 


sf’) EP) Phirt=0, ee (11) 


omitting the primes. 

The solution of (11) is complicated by the fact that the relationship between 
h and r involves f(r). The resulting differential equation can be solved, but 
leads to a power series which is only very slowly convergent and is of little practical 


~ 


|] Displacement of 
P PoP roots tere 


15 Field 
(microns) 


Fig. 8. Figured mirror. Fig. 9. Field distortion of figured mirror. 


value. However, an approximate solution may be found by numerical means. 
The steps are as follows: assuming the unperturbed curve ‘’ =2zr—1?, calculate 
h for a suitable range of values of r and an assumed value of J; calculate h/r? 
for each r and approximate by a polynomial of form A/7?= A,r"; substitution 
in (11) and solution of the resulting differential equation gives f(r) directly; 
substitute f(7) in (7) and check the design by computation. This has been done 
for /=1 andz=0-01. ‘The resulting equation for the curve is 
Y’ =0-027-42 — 0-0012457(47 — 1/1573). 
The computed form of field is shown in fig. 9, ‘The spherical aberration is 
nearly constant over the field and is such that a semi-aperture angle of 0-77 x 10-3 
can be used with a wavelength of 1:44, giving a resolving limit of 900A, and a 
corresponding depth of focus of 0-0118 cm; thus the field is adequately flat | 
over at least the extent shown. ‘The computation was carried out using a 
two-term approximation to h/r?. It is evident that further refinement of the 
design would give a better approximation to flatness of field. The total length 


of mirror employed for a focal length of 1 cm and a field of the extent shown 
in fig. 9 is approximately 0-7 cm. 
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It is found that the chief limitation on the design is that in order to reduce 
the spherical aberration the aperture stop must be nearer the focal surface, and 
_ hence a greater length of the mirror surface must be used for a given field. This 
is undesirable, for not only must a greater area of mirror be finished accurately, 
but also the ratio in scale of the two mirrors in mutually perpendicular planes 
must be greater, with a corresponding increase in the spherical aberration of 
the large mirror. 


§6. SYSTEMS EMPLOYING TWO MIRRORS IN EACH PLANE 
Harker (private communication) has pointed out that the field obliquity can be 
corrected by using two mirrors in each plane, with the concavities facing in 
opposite directions. An analysis of such compound systems shows that if the 
concavities face in the same direction both field obliquity and primary spherical 
aberration can be corrected. The mirror arrangement is shown in fig. 10. 


Fig. 10. ‘Two-mirror system. 


The object point O is imaged by the mirror M, at I and the second mirror M, 
forms an image of I at infinity. ‘The radii and glancing angles of the principal 
ray are R,, R, and i,, 7, for the mirrors M, and M,. To simplify the algebra the 
aperture stop is taken to coincide with M,. Under these conditions the angle 
between the principal ray and the focal surface for M, can be shown to be 
‘tan B=R,i,?/2l'. For the second mirror eqn. (6) can be applied, setting D equal 
to L, the distance between the points of incidence of the principal ray. ‘Then, 
- for the overall field obliquity to be corrected 


Ryiy? _ 1,(2L—%Ry) _ 2a," 13 

pee eRe Btls ee 
where L has been replaced by J,’ +7,R,/2, the condition for infinite final image 
distance. Equating the spherical aberrations of the two mirrors, 

R,= 81,4 R, = L,')/ Ry’. 

Substitution of this value in (13) gives 7, =R,2,3/41,'(3l,/—2R,). ‘Thus, as 
an example, if 7,=0-01, R,=100 cm, /,’=0-75 cm, we find R,=112-5 cm, 
a= 0-133, D=4-5 cm, /=1-5 cm, 

The performance of this system has been ascertained by computation. ‘The 
residual angular spherical aberration is very small, amounting to about 
1-7 x 10-* for a semi-aperture angle of 0-001 for a point in the centre of the 
field, and is approximately proportional to the fourth power of the aperture. 
~ Little use, however, can be made of this very small value of spherical aberration 
because of the field aberrations. ‘The most important of these are coma and 
curvature of field. 
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The field curvature can be corrected with the help of the coma if the aperture 


stop be placed at a position different from that of the first mirror. To explain | 


this it should be noted that from eqn. (1) it is possible to write ane combination. 
of the two aberrations as 


Angular aberration =4@,,x*h+@,,xh®, =... se (14) 


If the aperture stop be moved from M, but kept centred on the principal ray 
from the object point at zero field the beam from this point is not affected, but 
the principal ray from any other point with field coordinate / is now moved 
laterally in the beam emanating from that point by an amount proportional to h. 
This can be taken into account in (14) by replacing x by X=x—kh, where k: 
depends on the new position of the stop. Equation (14) then becomes 


Angular aberration =a,,(« — kh)*h + a,.(« — kh)h? 
= Ay X7h + xh?( ay. — 2de,k) + h3(ay,k" — aypk). ...... (15) 
In (15) the term a,,x2h is the coma and is unchanged. The term xh?(a,)—2a,,k): 


represents the field curvature, and can be made zero by suitable choice of k. 
The remaining term represents distortion and is not of serious consequence. 


Angular 
Aberration 
2x10 Ae = 37-5 microns ‘ 
110-5 = 18-75 microns 
A =0 
-0-00! -0:0005 0 0:0005 0-00! 


Aperture Angle 


Fig. 11. Aberrations of two-mirror system. 


The most suitable position for the stop is between the mirrors and distant 
approximately 0-4 cm from the point of incidence of the principal ray on Mj. 

Under these conditions the angular aberration is plotted in fig. 11 for three 
values of h. 

The maximum value of the glancing angle is 0-0143, corresponding to a 
wavelength of 206A. With this value, and allowing spherical aberration A/4 
at a field of 37:5 microns, the useful semi-aperture angle is 0-89 x 10-° radians, 
corresponding to a resolving limit of 11504. However, most probably it will 
be sufficient to allow the spherical aberration to reach A/4 at the half-field of 
18-75 microns; this makes the semi-aperture angle 1-3 x 10- and the resolving 
limit 800 A. 

§7. GENERAL DISCUSSION OF RESULTS 

It can clearly be seen from the above that the resolving power of the x-ray 
microscope cannot easily be varied by more than a factor of about two from the 
value given by a simple pair of crossed circular mirrors. This was forecast by 
Prince, although his value for the resolving limit, which was based on geometric 
optics alone, appears to be a little pessimistic in view of the true distribution of 
light in the image as given by wave theory. 
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The chief advantage to be gained from the use of systems more complex than 
that of Kirkpatrick and Baez would then seem to lie in the possibility of correction 
of field obliquity. The single circular mirror with a correctly placed aperture 
stop may prove to be the most useful on account of its simplicity. It has the 
advantage that, as the number of reflecting surfaces is a minimum, scattering 
will also be a minimum; furthermore, it seems likely that circular (or spherical) 
surfaces can be more accurately made than figured ones. Very careful experimen- 
tation would be necessary to prove the reality of the difference between the 
1500A obtainable by this system and the lower values obtainable by more 
complex systems, though it may be that the additional factor of almost two to 
be gained is worth the effort. 


§8. CONCLUSIONS 

The aberrations characteristic of mirror systems used far from their axes of 
symmetry have been considered, and especially the effect of spherical aberration 
on the Airy pattern of a single object point has been evaluated. 

An investigation has been made into methods of correcting these aberrations. 
It has been indicated that the field obliquity of a single pair of crossed circular 
mirrors can be corrected by a suitably placed aperture stop. 

A method has been devised of designing a figured mirror which, when used 
with an aperture stop, eliminates field obliquity but has less spherical aberration 
than a single circular mirror. 

By using two circular mirrors in each plane it is possible to correct field 
obliquity, field curvature and spherical aberration, and the design equations 
and computed results are given for such systems. 

The more complex systems considered do not increase the resolving power 
by as much as a factor of two, however, and it seems possible that the most 
rewarding system to use will be the simple circular mirror with suitably placed 
aperture stop. 
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ABSTRACT. The focal properties and spherical aberration constants of a number of 
aperture electron lenses of varying diaphragm thickness and varying strength were deter- 
mined numerically using potential distributions obtained (a) from a rigorous solution due 
to Ollendorff, and (4) on a resistor network analogue. The focal lengths of ‘weak’ aperture 
lenses of finite diaphragm thickness agree closely with those given by Davisson and Calbick’s 
formula for the ‘infinitely’ thin aperture lens; ‘strong’ aperture lenses of all thicknesses 
show systematic deviations from this formula. The primary spherical aberration constants 
of diverging aperture lenses are small, and of such sign that they cannot be used for the 
correction of the spherical aberration of converging electron lenses. 


§1. INTRODUCTION 

N aperture or ‘pinhole’ lens exists where a diaphragm pierced by a ~ 

AN circular aperture separates two spaces with different electric field 
strengths #, and E,. The geometry of such a lens is shown in fig. 1. 

When the first field strength £,=(U)—U,)/d, between electrode A, and 

diaphragm D, of thickness T and aperture radius R, is greater than the second 


Fig. 1. 


E, =(U,—U,)/d,, between diaphgram D and electrode Ag, the lens is divergent 
and, conversely, when the first field strength is less than the second, the lens 
is convergent. ‘The aperture lens of the first kind is the only known practical 
electron lens which is divergent; this type of lens is found in most accelerating 
systems, such as those used in electron guns. 

Davisson and Calbick (1931, 1932) have shown that an aperture lens, with 
an indefinitely thin diaphragm and with small field strengths E, and E,, and a 
potential U)=U applied to the diaphragm, has the properties of a thin lens 
located in the plane of the diaphragm whose focal length is given by 


f=4U/(E,—;). ieee) 


* Now at Carnegie Laboratory of Physics, University College, Dundee. 
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If one of the fields is zero, e.g. E,=0 and E,=£, eqn. (1) can be written in the 
form //R=—485, where the dimensionless parameter 5=U/ER is a measure 
| of the strength of the lens. The most representative values of 5 used in practice 
lie between 10 and 20. The condition under which eqn. (1) was derived is 8 >1. 
Bedford (1934) discussed more fully the interpretation of eqn. (1) for finite 
values of 5. Gans (1937), using a numerical method, found small systematic 
deviations from the focal length predicted by formula (1) over the range 
e=3-5 to d=12. 
_ The present work was designed to evaluate, in a few cases, the spherical 
aberration constant of diverging aperture lenses about which nothing was known. 
It was also desired to obtain the focal length for values of § which are no 
longer large. 


$2. FIELD DISTRIBUTION WITHIN. THE LENS FIELD 


The focal length and spherical aberration constants were evaluated from 
electron trajectories traced through the lens fields. The potential distribution 
around an indefinitely thin diaphragm, T=0, perforated by a circular hole, 
has been calculated rigorously by Ollendorff (1932).* In the case of a diverging 


Fig. 2. Diverging lens T=9, axial potential values and first four derivatives (u=1, 6=3). 


lens where the space to the right of the diaphragm is field-free, E,=0, 
Ollendorff’s solution for the potential along the axis can be put into the form 


CATs AH Tame aris il (2) 
In eqn. (2) €=2/R, the origin being at the diaphragm, and 
e(=slemAsCtan fim, hee (3) 


Equation (2) is valid under the condition d=d,>R, i.e. 5>1. In practice 
the deviation of Ollendorff’s solution from the true potential distribution is 
negligible for 6 >3 and still small for 6=2. (The true value of 4(¢) is zero for 
f=—6.) 

Figure 2 shows curves of 4(¢) and its first four derivatives ¢’, 6”, ¢’”, and 
4°, where the primes denote differentiation with respect to ¢; in fig. 2 U=1 
and 6=3, 


* An identical solution for this problem was given by Fry (1932) and by Glaser and Henneberg 
(1935). 


 — 
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a | 
If E,=0 and E,=E>0, the aperture lens is convergent. Then Ollendorfi’s 
solution takes the form j 


$(f)=U+ UlE—g(2)I/6. renee (A) 


Figure 3 shows values of 4(£) and its first two derivatives ¢’ and ¢” for this case. ; 
For lenses of finite thickness, 7-40, the potential distributions were obtained — 
on a resistor network analogue* (Liebmann 1950). ‘To reduce the errors due to 
the finite mesh size the extrapolation method described by Liebmann was used. 
The equivalent model has a distance between diaphragm and each anode of 


-30 =20 =1-0 0 10 20 30 


Fig. 3. Converging lens T=0, axial potential values and first four derivatives (u=1, 5=3). 


about 160 mesh units and an aperture diameter of 40 mesh units. To obtain 
an indication of the remaining influence of the mesh size on the potential values 
measured on the resistance network analogue, potential values were compared, 
in the limiting case of the indefinitely thin diaphragm (T'=0), with Ollendorff’s 
exact solution, eqn. (2). ‘The extrapolated network values for ¢(¢) were found 
to agree with the calculated ¢(¢) values to within 0-0002U or better. The 
derivatives of the axial potential distribution 4(z), needed for the computation 
of the trajectories, were obtained from the measured potential values by 
differencing in the usual manner. 


53: DRATEC TORY TRACING 


Paraxial trajectories were traced with the help of Liebmann’s (1949) 
recurrence formulae 


Tati = OW, ty Oar,’ Pena, = Osr si: Oy, eee (5) 
where the QO are a power series in Az. 
For an initially parallel ray (7)’ =0) the focal length f in image space is given 
by 7)/r{. In a diverging lens the position 2, of the virtual image is obtained by 
producing the paraxial ray back to cut the axis. 


_* This was a new precision network with 160 meshes in the z direction and 50 meshes in the 
y direction and increased accuracy. 


Focal Properties of Aperture Electron Lenses 593 


The spherical aberration constant C, can be obtained from the transverse 
aberration Av, of the ‘extra-axial’ ray from the paraxial ray. ‘The recurrence 
formulae are then 


Arng =QAra+ QeArn, Aras’ =QsArq’ + QsAt,+ Qs, 00-0: (6) 
where Q, is an expression containing the derivatives of ¢ up to the fourth order, 
and the r, and r,,’ values of the paraxial ray. 


7 


§4. PARAXIAL PROPERTIES OF LENSES 
4.1. Diverging Lens, E,=E, E,=0 
(i) Thin lens, T=0. In the limiting case of the thin diaphragm Ollendorft’s 
potential values were used (eqn. 2). Trajectories were calculated for different 
values of the parameter 6, characteristic of the strength of the lens. A typical 


Table 1. R=0-1cm 


8 AIR a/R flom— ney _—dDish.(%) 
= Y — 6:78 —6:°75 —0-68 —0-80 20 
5 —18-51 —18-:57 —1-85 —2-00 8 
10 — 38°51 — 38-55 —3-85 —4-00 4 
20 —78°5 —78°5 —7°85 —8-00 D 
50 —198-4 —198-4 —19-8 — 20-00 1 
1 


100 —398 — 398 BS: —40-00 << 
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ny 
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Fig. 4. Paraxial trajectory showing 2} Fig. 5. Focal length and aberration 
(6=U/ER=2). constant as function of 6. 


trajectory is shown in fig. 4 ford6= U/ER=2.* ‘Table 1 gives f and 2, for various 
| values of 8 between 2 and 100, and fig. 5 gives a plot of fas a function of 8. 
Agreement with the approximate formula, eqn. (1) (dashed line), is shown 
to be very close at high values of 6, i.e for ‘weak lenses’; this was expected as 
» eqn. (1) was derived under the condition 6>1. For lower values of 6 (‘strong 


* Such a small value of 6 was chosen to illustrate the shape of trajectories, as the refractive 

- power of the lens is correspondingly great; this shows up better the bending of the trajectory in 

the lens field. Also, for values of 6>2, the virtual image point, coordinate 2;, would not be shown 
as — 2; would be much too large. 
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lenses’) systematic deviations from the Davisson and Calbick (DC) formula | 


are found, the lens being always stronger than predicted by this formula; thi 


deviation amounts to as much as 20% at 5=2. The results given in fig. 5 and | 
table 1 agree closely over the range 5=3-5 to 12 with those obtained by Gans ~ 


(1937). . H 

he principal planes coincide with the centre of the lens for large values. 
of 8. As the lens increases in strength they move slightly away from the centre_ 
by distances of the order of 0-04R. 

(ii) Lenses with finite thickness, TAO. To study the dependence of the 
focal length, for given field strengths Z, on the diaphragm thickness, lenses of 
diaphragm thickness T=0-2R, 0:4R, 0:8R and 1:6R were investigated. Paraxial 
_ lens constants were determined for the relatively thick lens T=0-8R for the 


Table 2. f/R for U=50kv, R=0-1cem 
Thickness  20kv/em 50 kv/cm 100 kv/cem 


0 98-4 38-6 18:5 
0-2R == 37-9 as 
O-4R — 38-1 5 
0-8R 98-4 38-6 18-6 
1:6R — 38-2 = 


range of 6=5 to 25, but only for 6=10 for the other lenses. ‘The results given” 
in table 2 show that the focal length is almost identical with that of the lens with 
an indefinitely thin diaphragm, T'=0, so that the curve shown in fig. 5 can be 
taken as representative of all divergent aperture lenses whatever the diaphragm 
thickness. 

4.2. Converging Lens, T=0. 


Using potential values obtained from eqn. (4), the focal length of a converging 
lens was calculated for the case 5=5 with a potential U=50kv applied to the 
diaphragm, the space to the left of the diaphragm also being at this potential. 
As the aperture lens is followed by a space of changing refractive index n, the 
focal length f depends on the position where f is measured, as f/m is constant 
in this space. If the position z for measuring f is chosen where a ray incident 
parallel to the optical axis crosses the axis, 2,=4:4cm, then f= +6-6cm. If the 
optical system were assumed to be a ‘thin’ pinhole lens followed by a parallel 
field refracting the rays parabolically, and Davisson and Calbick’s formula 
were applied to the pinhole lens, then one would find 3; =4-0 cm and f= +6-0 cm. 
Hence the deviation of the numerically determined focal length from the focal 
length calculated from Davisson and Calbick’s approximate formula, in the 
direction of greater refracting power, is only slightly greater (10%) in this case 
than for the divergent lens (8%). 


§5. SPHERICAL ABERRATION CONSTANT C, 
5.1. Diverging Lens of Diaphragm Thickness T =0 


A typical plot of the deviation Ar of the ‘extra-axial’ ray from the paraxial | 


ray, showing that the deviations reach a maximum near the centre of the lens, 
is given in fig. 6; the initial off-axis distance of this extra-axial ray incident 
parallel to the optical axis is 7)/R=0-1. The spherical aberration constant C, 
is defined by the relation Av; =C,«3, where Ar; is the transverse aberration and 
a the angle of the ray with the optical axis. Table 3 gives the values of C, for 


ioe, Aico alae sikh diene 
4 * 


| 
| 
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_ various values of 5 from 2 to 50, and in fig. 5 C, is plotted as a function of 6. For 
higher values of 5 C, follows a square law, Orsi: 162, but for small values of 6 
C, is smaller than this relation (broken line) would predict. 

In order to compare the aberrations of various types of lenses and to show 
the connection between focal lengths and aberrations it is convenient to plot 
the dimensionless aberration constant S=C,/f, dependent on the lens geometry, 
against the relative focal length f/D, where D=2R. If the results obtained 
experimentally by Liebmann (1949) for unipotential electron lenses, and by 
Bruck and Romani (1944) and Ramberg (1942) using numerical methods, are 


Zr 
30 2.0 -10 0 10 20 30 


a eeae 
x10") 
01 10 t/D 10 100 
° Pia phureae Lens (Ramberg) = Magnetic Lens (Liebmann 
Liebmann & Grad) 
° i » (Bruck &Romani) a Aperture Lens 
Fig. 6. Deviation of first-order ray from Fig. 7. Comparison of lenses. 
paraxial ray. 
Table 3. R=0O-1lcm 
a) CJR Cy (cm) S=C,/f f/D 
2 +20-7 +2-07 —3:05 — 3-39 
5 +177:1 +17:7 —9-57 —9-26 
10 +731 +73-1 —19-0 —19:3 
20 +3060 +306 — 39-0 — 39-2 
50 +19500 +1950 —98°-3 —99-2 


plotted on a double logarithmic scale, straight lines of slopes between 1-75 and 3 
are obtained (fig. 7). In the case of the divergent aperture lens the relation 
between 6 and f/D is nearly linear (slope =1 in fig. 7), and the value of S is much 
smaller than in the other lenses. 

In simple symmetrical divergent glass lenses the spherical aberration is 
always of such sign that the marginal focal length is smaller than the paraxial 
focal length. However, in the divergent electron lens investigated here the 
opposite has been found, the marginal focal length being greater than the 
paraxial focal length. In this and in the very much smaller absolute value of the 
spherical aberration constant C, this divergent electron lens differs from all other 
electron lenses studied so far. 

In view of this difference a check by Scherzer’s (1936) integral formula and 
by Goddard’s (1944) method was carried out for the case 6=10. ‘The spherical 
aberration constant obtained by Scherzer’s formula was identical in sign and 
value within a few per cent, whereas Goddard’s method gave a value of the same 
sign but about 25% lower in value.* 


* This appears due to the fact that in Goddard’s step-by-step method the past three points are 
involved in every step forward. This leads to a smoothing out of the effect of the sharp peaks in 
the higher derivatives of the field distribution, which are mainly responsible for the aberrations, 
unless a greatly refined interval size is used. 

2R-2 
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Because of the ‘sign’ of the spherical aberration constant of the divergent 
aperture lens, one cannot correct the spherical aberration of a converging lens 


by combining it with a diverging aperture lens. 


5.2. Lens of Diaphragm Thickness T =0-8R 


The spherical aberration constants for lenses of finite thickness were — 


determined only for the very thick lens with 7 =0-8R, as the computation of the 
spherical aberration constants is rather tedious. In table 4 values of C,/R are 


Table 4. C,/R for R=0-1cem, U=50kv 


E (kv/cm) L—( T=0°8R 
20 +5000 +2750 
50 st! +642 
100 se dlg7/ +166 


compared with those of the lens JT =0, and in fig. 5 C, is plotted as a function of 6 
for R=0-1cm. One sees that the spherical aberration constants of this very 
‘thick lens are of the same sign but are slightly smaller than those of the 
‘indefinitely thin’ lens. For very weak lenses the difference is of the order of 
20 to 50%, but at the greater refracting powers the difference becomes very 
small. It is therefore reasonable to assume that over the usual working range and 
for the diaphragm thicknesses used in practical lenses the curve given in fig. 5 
for the spherical aberration constant of the ‘indefinitely thin’ lens can be applied. 


§6. CONCLUSIONS 


Focal lengths of aperture lenses of finite diaphragm thickness vary only 
very slightly with thickness, so that in practice the Davisson and Calbick 
‘thin lens’ formula may be applied directly in the case of the diverging lens, 
and with a suitable modification, taking the following parallel field into account, 
in the case of the converging lens. 

The spherical aberration constant does not seem to vary appreciably with 
the diaphragm thickness, particularly at relatively high refracting powers; for 
the diverging aperture lens met in practice it is very small and of sign opposite 
to that found in simple diverging glass lenses. 
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ABSTRACT. Examination of the requirements for reducing the spherical aberration of an 
electron lens by modifying the field shows that the hypothetical lens field H(0, z)=y/z” 
should have small aberrations. 

It is shown that the properties of a certain Glass of this general field configuration can 
be described in terms of an exact analysis. In particular, the electron optical properties of 
the y/z* field are described in detail. 

The y/z* field is found to have small aberrations, but practical conditions limit appreciable 
improvement in resolving power. The characteristics of the y/z° field are also examined and 
are found to give little further advantage. 


§1. INTRODUCTION 

HE spherical aberration of electron lenses has been the subject of extensive 
| study, and it has been established that a positive aberration is an inherent 
property of lenses of the type in present use, the aberration being produced 
by an increase in the refractive power of the lens field with the radial distance 
from the optic axis (Scherzer 1936). It follows that the aberration cannot be 
corrected by combining existing lenses, and it remains to develop lenses which 

give a minimum of aberration. 

In principle the spherical aberration can be reduced to any value by suitably 
modifying the lens field, as is shown for example by the influence that asymmetry 
of the axial field has upon the aberration of the magnetic lens (Dosse 1941 a). 
Indeed, Glaser (1940 a) has derived the form of the magnetic lens field that gives 
zero aberration, but the refractive power of this particular field is too small to 
be of practical importance (Rebsch 1940). Accordingly, it is of interest to 
examine the requirements for the reduction of the aberration of a practical lens 
and to determine the characteristics of a suitable lens field. 

The approach to this problem is indicated by an examination of the general 
expression giving the spherical aberration in third-order approximation (Glaser 
1940 a), namely 

A= in|, [Lr,*+2Mr,?r,/2 + Nr,'*] dz, 
_ where ry’ is the slope of the trajectory at the object point z =X, and r,(z) is the 
specific paraxial trajectory that originates in the object point with unit inclination 
to the optic axis. L, M, and N are functions of the axial distributions of the 
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electric and magnetic components of the lens field, 4(0,2) and H(0,z) | 
respectively, and their derivatives, i.e. 


L=o~ [s (+ 3, it i) 19%: < HH" "| 
majo: eee 


Many of the variables entering into this expression are interdependent and — 
it is difficult to determine a practical approach for minimizing the aberration, ; | 
except to specify that a strong negative term in the integral is required. For a _ 
purely electrostatic lens this requirement is partially fulfilled if the fourth and 
second derivatives of the field are respectively positive and negative in the region 
between the object and the image. In a continuous lens field these two 
requirements are satisfied mutually only over a limited region about the maxima 
of the distribution. Thus, this condition, together with the simultaneous 
requirement for a rapid convergence of the electron trajectory, indicates that the __ 
optimum field is one that includes a sharp maximum, with the object situated in — 
that region of the field which is varying rapidly. Unfortunately, the field gradient — 
is limited in practice by auto-electronic emission from the electrodes, and it is 
also undesirable to locate the object in a region of high field strength. It appears 
that there is little improvement to be obtained in this manner. 4 

The corresponding condition for the aberration of a purely magnetic lens to 
be a minimum requires that the second derivative of the field be highly positive 
in the region between object and image. ‘This is also a necessary condition for 
the aberration of a composite electric and magnetic lens to be a minimum, the 
consideration given to a purely electrostatic lens also applying. These 
conclusions are in general agreement with those given by Rebsch (1938), who 
showed that the aberration may be reduced by immersing the object in a more 
rapidly varying field. It should be stressed that application of these relations 
is very limited by practical restrictions to the field intensities. Nevertheless, some 
small improvement is possible. These general conclusions will be considered 
in relation to the purely magnetic lens in order to ascertain the limit to practical 
improvement. 

The field of a typical magnetic lens is approximated to by the analytical field 
H(0, z) = H(0, 0)/{1 + (z/a)?}, and it is often convenient to discuss its performance 
in terms of this analytical field (Glaser 1941, Dosse 1941 b, Hesse 1950). The 
second derivative of this axial field is positive only for values (z/a)?>3. Thus, 
at the centre of the lens where maximum refraction occurs, the second derivative 
is negative so that a positive contribution to the basic spherical aberration results. 
This is indicated by the corresponding influence of the object position (Glaser 
1941), 

Consider now the hypothetical field having the general form H(0, 2) 1/3”. 
The second derivative of this axial field is always positive. Accordingly, this 
lens field should give a smaller aberration than does the 1/{1+(z/a)"} field; 
for this reason its characteristics have been examined in detail. In practice this 
field cannot be reproduced completely because of the discontinuity. However, 

a part of the field can be produced by a suitable choice of lens geometry. Thus, 
if the object is regarded as being self-luminous, as is often permissible, the 
characteristics of the lens may be determined by consideration of the analytical 


( 


Z Magnetic Lenses having the Axial Field H(0, z)=y/z" 599 


field alone if the object be placed in that part of the field approximating to the 
theoretical distribution, provided that the entrance aperture can be suitably 
adjusted. A description of the properties of the 1/z” magnetic lens field forms 
mthe content of the remainder of this paper. 


§2. ELECTRON MOTION IN THE 1/2” MAGNETIC LENS FIELD 


An analytical solution to the electron trajectory in a lens field has been 
obtained for only a few lens fields (Glaser 1941, 1949, Hutter 1945, Riidenberg 


- 1949). Therefore, it is interesting to find that a certain class of the general 


1/s” lens field is amenable to exact analysis in terms of elementary quantities. 
The equation of motion of a paraxial electron in the magnetic field 


_ H(0, z)=y/s", where y is a constant characteristic of the lens, may be expressed 


or + BPr/28" =O, loons: (1) 
where B=(ey?/8m¢z)"” is a constant which determines the power of the lens, 


and has the dimension of length alone. By making the substitution u=r'/r, 
the ray equation may be reduced to a differential equation of the first order in the 


_ form of the Riccati equation, and is integrable in finite terms if 


CREO CL NG es, 2 enn he ip ara Se (2) 
where g is zero or a positive integer (Forsyth 1929). 

Thus, an analytical solution to the paraxial electron trajectory can be obtained 
for fields that are defined by eqn. (2). Of these the 1/z? field is of particular 
interest since it has the strongest divergence. This lens field may be 
approximated to by a lens having the geometry illustrated by fig. 1, with the 
object mounted in the divergent part of the field, this arrangement allowing the 


» object to be situated in a region of high field intensity. The electron imaging 


properties of this particular lens field are examined in the following section. 


$3. THE DIOPTRICS OF THE 1/27 LENS FIELD 
The trajectory of a paraxial electron in the 1/z? lens has the solution 
r=Dz2z((1+K) cos (B/z)+i(1—K) sin (B/z)); is se (3) 
D(1+4K) and D(1—K) are constants of integration determined by the initial 
conditions that specify a particular trajectory. The location of the gaussian 
image plane is given by the point of intersection with the optic axis of the SGes 


paraxial ray 7,(z). Equation (3) gives for this trajectory 7,(3)= —(z/q) sin %, 
where w=B/X, b=B(z1—X- ). It follows that paraxial trajectories intersect 


the optic axis at points given by the relations 


RRO eels Bema ney ONT. FF wads (4) 


This value of z and the positive values of % represent the infinite number of 
virtual intersections of the trajectory with the axis in the region preceding the 


object point. Hence these values do not influence practical considerations 


and will be neglected. ‘The N negative values of ys correspond to N successive 
image points. Consequently the condition for single image formation may be 
expressed 27 >w2>7. Accordingly, the position of the gaussian image plane 


under conditions of single image formation is given by the expression 


2,=B/(w—7z). The corresponding magnification of the gaussian image is given 
by the relation m=w/(w—7) where w is again to satisfy the condition expressed 


for single image formation. 
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It will be observed that when w=7 the image is formed at infinite distance - 
with infinite magnification, the corresponding object point X;=B/7 being the 
effective focal point of the lens. Accordingly, the effective focal length of the 
lens is given by the relation f=lim,_, .7(2) where 74(z) is the paraxial ray that 
originates in the focal point with unit inclination to the axis. Accordingly, — 
eqn. (3) gives for the effective focal length f= B/7, giving a simple measure of the 
refractive power of the lens. 

The general gaussian properties of the lens are now determined explicitly 
in terms of the field parameters. It remains to evaluate the aberrations. _ 


24 
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Fig. 1. The diagrammatic cross section Fig. 2. -The variation of the least resolvable 
of a magnetic lens which is suitable distance Ar with the lens power parameter B 
for reproducing a part of the of the 1/z? magnetic ‘lens field at specific 
hypothetical field distribution electron energies. The power of the lens is 
H(0, 2)=y/z?. maintained at the high magnification condition 


throughout, and the loci of the field strengths 
that are appropriate to this condition are given 
in kilogauss. 


§4. THE RESOLVING POWER OF THE 1/z? MAGNETIC LENS FIELD 


For this particular lens field the spherical aberration, defined above in third 
order approximation, becomes 


aC Z| (w— m+ 5 43]. Sto (5) 


In addition to the spherical aberration, a chromatic aberration produced by 
fluctuations eA¢ in the energy of the electron beam is important in considerations 
of resolving power. Fora purely magnetic lens this aberration may be expressed, 
in third order approximation (Glaser 1940 b), 


Ad 44 
A= — Try’ i Pred ee ee rs Ss Ais 
3 I 


where 7,() is again the trajectory that originates in the object point with unit 
inclination to the axis, with the added condition that it vanishes in both object 
and image points. For the 1/z* lens field this aberration may be expressed 


A,/rx' = CAg/} =(7/2B)X*Ad/d. 
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A comparison of the quality of different electron lenses is difficult because 
of the many independent parameters involved, and it is customary to take the 
dimensionless constants C,/f and C,/f as characteristic figures representing the 
aberrations produced by lenses of equal focal length. For the 1/s? lens field 
with the object maintained at the high magnification position, these characteristic 
figures take the values (C,/f)y,=0-18, (C,/f)x,=0:5. The corresponding high 
magnification characteristics of the 1/{1+(z/a)?} lens field are dependent upon 
the parameter k® (Glaser 1941), but for a lens power maintained at the optimum 
value for a specific voltage (Cosslett 1946) the aberration characteristics take the 

values (C/f,)x,=0-96, (C,/f)x,-=0-76. In terms of these characteristics alone 
the 1/z? lens field gives a marked improvement. However, little indication is 
given of any improvement in resolving power that may be obtained in practice, 
bearing in mind the limitation set by magnetic saturation of the pole pieces. 
A more satisfactory comparison is given by the limit of resolution obtained under 
high magnification conditions with the lens aperture adjusted to equalize the 
discs of confusion produced by spherical aberration and by diffraction 
(Zworykin et al. 1945). This limit A, is plotted for the 1/z? lens in fig. 2, in 
dependence upon the power parameter B for specific values of the accelerating 
voltage ¢. ‘The contribution of the chromatic aberration is calculated for a 
voltage ripple of 2-5 x 10-5, and a spread of one electron volt in the energy of the 
imaging beam. The loci of the appropriate field strengths at the object, H(0,X;), 
that are necessary to maintain the high magnification condition are also shown. 

A comparison of the resolving power of this lens field with that given by the 
1/{1+(z/a)*} field (Cosslett 1946, 1947, Gianola 1950) shows that it is little 
different. It is not possible to ascertain the absolute ultimate resolving powers 
of the two lens fields without some knowledge of the practical field intensities 
that are feasible, but it is expected that the peak field strength of the 1/z? lens field 
may approach saturation field strength, since the object may be placed close to 
the pole pieces. Consequently some small improvement in resolving power 
may be possible. 


§5. FURTHER INVESTIGATION OF THE GENERAL 1/2 LENS FIELD 


It follows from eqn. (2) that the fields that can be analysed in the manner 
described all have a field divergence equal to or less than that of the 1/2? lens. 
The spherical aberration of the latter field was not expected to differ appreciably 
from that of the 1/{1 + (z/a)*} lens field with an infinite asymmetry index, g= © 
(Dosse 1941 a). For this reason the performance of a rather more divergent 
field has been examined. 

Referring again to the paraxial ray equation, eqn. (1), and making the 
substitution uw=1/z, p =1/2(n—1), the ray equation is reduced to the form of the 
general Bessel equation. Accordingly, for non-integer values of p the trajectory 
may be expressed in terms of Bessel functions of the first kind, namely 

raw (A,T()+ AyT (x), nea (7) 
where A, and A, are integration constants and «=2Bpu'?”. The position of 
the gaussian image plane can be determined as previously by finding the axial 
intersections of the specific trajectory 7,(z); this gives the general relation 
Ff (x;)I_(%_) =I _(*,)J_,(%,) where x, and x, represent the values of x at the 
object and image points respectively. 
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Further, the requirement that the trajectory should approach zero slope at 
infinite distance gives for the focal point the relation J ,(xp)=0. The effective 
focal length of the general 1/z” field may now be obtained as previously, and 
can be expressed in the form 


es Cena 
x.~~ 3) saphampa 70 


where (7) is the gauss pi-function. | 

It remains finally to obtain an explicit relation for the spherical aberration | 
that is produced by the general 1/z” lens field. Substituting for the axial field 
(0, z)=y/z” and for the paraxial trajectory, eqn. (7), into the expression for the 
aberration constant 


f-sey aap) | |"* temas 
Hy (I g(8) —Tgl)T (0) #OP devas (8 


It is now possible to determine the main characteristics of a specific example : 
of the general 1/z” lens field, given tabulated values of the ap PSOpE=s Bessel ~ 
function. For the particular instance of the 1/z? lens field, p=, this general 
analysis yields expressions which correspond to those obtained by direct 
integration of the ray equation. 7 

In particular the 1/z? lens field has been investigated, it being anticipated — 
that distributions of greater divergence would have too small a refractive power 
to be of importance. For this distribution, p=}, an inspection of tables of 
J ..1j4(*) (National Bureau of Standards 1949), gives for the condition of single 
image formation 5-91 >B/2.X7>2-79 and for the position of the first effective 
focal point x-=2-79. In addition, a numerical integration of eqn. (8) between — 
the limits x,,=2-79 and x, =0, gives for the aberration constant of the 1/2 lens — 
field, under alinane ue high magnification, the value C,/X,;=0-10. 


§6. COMPARISON OF EXAMPLES OF THE 1/2” LENS FIELDS 
The quantity C,H(0,X;) is best taken as a characteristic figure in the © 
comparison of 1/z” lens fields, as it gives a measure of the aberrations produced 
by fields having objects situated at points having comparable field strengths. 
The table gives this aberration characteristic, together with the effective focal 
point and focal length characteristics, for the 1/s*°, 1/s?, and 1/z° lenses, to show 
the variation over a practical range. 


n 4/3 2 3 

xe 4-49 3-14 2:79 
(e/8mdx)!? XH(0, Xp 1-50 3-14 5-58 : 
(e/8mdx)"? fH, Xp) 27-6. 3-14 2-52 : 
(e/8my)! C,H(0, Xp) 2-6 0-57 0-56 


It is apparent that little practical improvement in the quality of the lens is 
obtained by increasing the divergence of the field beyond that of the 1/z? lens. 
For this reason the resolving power of fields of high divergence will not be | 
considered further. 

This analysis gives some indication of the lower limit that the spherical | 
aberration may be reduced to in practice by a straightforward modification of 
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the field of an electron lens. It is concluded from the results of the investigation 
_ that, in the case of the magnetic lens at least, no outstanding improvement is to 
_ be obtained in this way without introducing subsidiary corrector fields. 

The fact that the 1/2” lens fields can be subjected to rigorous analysis 
facilitates their use in design problems, in particular those concerned with the 
investigation of suitable corrector lenses, when their low aberration characteristics 
may prove advantageous. 
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ABSTRACT. Expressions are derived for the intensity of the analysed light reflected from 
absorbing crystals (a) under crossed polarizing units and (b) with the analyser slightly 
rotated from the crossed position. Four intensity maxima per complete rotation of the 
specimen are obtained when the polarizers are crossed, but under certain circumstances, as 
the analyser is rotated from the crossed position, the number of maxima may be reduced to 
two. ‘These relationships are illustrated by experimental curves on specimens of tin, 
uranium, bismuth and stibnite. 


si. INTRODUCTION 
N an earlier paper (Mott and Haines 1951 a) the effect of crystal symmetry 
on the reflection of a plane polarized beam of light at normal incidence was 
summarized, together with the modifications produced by anisotropic surface 
Ims and controlled etching of the reflecting surface. 


Eh 
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According to the theory developed by various authors (Szivessy 1929, Wright | 
1919, Orcel 1925, Capdecomme and Orcel 1941), four equal intensity maxima 
and minima should be obtained when an anisotropic specimen is rotated over _ 
a range of azimuth angles from 0°-360° between crossed polarizers. This is — | 
confirmed by observation in some cases, but many anomalies have been reported. — | 
Farnham (1931) and Short (1940) classify anisotropic minerals according to 
whether they give four or two intensity maxima when the crystal is rotated through 
360°. Eash and Upthegrove (1933) state that the anisotropic 5’ constituent in — 
ternary copper-—nickel-tin alloys gives two extinctions per complete stage 
rotation. Sampson (1923) measured the rotation angle produced by reflection 
from various faces of a manganite crystal, unequal maximum rotations in the 
positive and negative directions being obtained. In addition, he found that the 
angle between successive intensity maxima and also between successive minima 
was not exactly 90° as would be expected. Sampson (1929) subsequently 
amplified his results and showed that as the analyser is rotated from the crossed 
position intermediate maxima are reduced in intensity until finally only two 
maxima are obtained when the offset of the analyser is equal to the maximum 
angle of rotation of the plane of polarization. The maximum contrast between 
the grains however is obtained when the analyser is rotated still further. Orcel 
(1925) also points out that, when the effect of reflection is essentially a rotation 
of the plane of polarization (i.e. the ellipticity of the reflected beam is very ~ 
slight), the contrast between the grains in a multi-grained specimen is most 
marked when either the polarizer or analyser is rotated slightly from the crossed 
position. Jones (1924) found that most etched metals gave four maxima, but 
with bismuth lightly etched in nitric acid a marked difference in the alternate. 
pairs of maxima can be detected when the analyser is rotated only slightly from 
the crossed position. Jones suggests that this affords an accurate method of 
crossing the polarizer and analyser. On the other hand, Morrogh (1949), for 
example, reports that the anisotropic 7-phase of the iron-silicon system gives 
four maxima. 

Variable results were obtained in the earlier work at A.E.R.E., especially 
on electrolytically polished uranium, although some specimens always gave four 
intensity maxima per rotation. For instance, in a given field of a uranium : 
specimen, some grains gave two intensity maxima, others gave four, and only | 
in a few cases were the four maxima of equal intensity, two being more marked 
than the other pair. ‘The maximum grain contrast for uranium and many other 
anisotropic metals as polished was obtained when the analyser was rotated 
through 1°-2° from the crossed position. When a Foster prism (1938) was 
used as the vertical illuminator, the results were generally similar to those with 
polaroid units and a glass slip vertical reflector, although when four maxima 
were obtained they appeared to be more equal in intensity. If a quarter-wave 
plate was inserted between the Foster prism and the objective, it was often 
possible to rotate the plate to such a position that four nearly equal maxima 
were obtained for most grains on a particular specimen. The setting of the 
quarter-wave plate varied with the specimen under examination and the 
objective lens used. 

In view of these anomalies it was decided to investigate the phenomenon: 
both in theory and in practice. The intensity measurements were made with. 
the photometer unit previously described placed in the monocular eyepiece tube 
of a Beck Universal Metallograph No. 50 (Mott and Haines 1951 a). 
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§2. EXPRESSIONS OF PRACTICAL INTEREST 
2.1. Low Symmetry Crystals 


2.11. General case for low symmetry crystals (i.e. triclinic, monoclinic and 
orthorhombic systems). Whilst transparent crystals of these three systems 
are optically biaxial, absorbing crystals with the same symmetry do not possess 
true optic axes. Fora given direction of the wave normal, such crystals in general 
transmit two elliptically polarized waves. The ellipses representing the 
polarization of the two waves have the same ellipticity, and are described in 
the same sense; the major axes are mutually perpendicular. 

A normally incident wave reaching the surface of an absorbing crystal enters 
it without any change in the direction of the wave normal. ‘Two principal 
elliptically polarized waves can therefore be associated with each case of normal 
incidence, and will vary with the orientation of the surface and crystal symmetry. 
These two principal refracted waves are transmitted with different velocities, 
and are absorbed to different extents; asSociated with them are two principal 
reflected waves with different reflectivities and phase shifts. 


Fig. 1. Fig. 2. 


Only a principal incident beam can be transmitted and reflected without any 
change in polarization. In dealing with a plane polarized incident beam it is 
first necessary to resolve it into two principal elliptically polarized beams. 


2.12. Resolution of plane-polarized normally-incident beam into two principal 
elliptically polarized beams. ‘The axes of coordinates for the surface are most 
conveniently taken to coincide with the axes of the principal ellipses of 
polarization for a wave normal perpendicular to the surface. 

Let the azimuth of the plane of polarization of the normally incident beam 
be ¢ and let the amplitude of vibration be c (fig. 1). The vibration representing 
the plane polarized beam is 


W=ccosocos?: y=csindcos@ = -  kasuien (1) 
where 6=2zft, f being the light frequency. 


The two principal elliptically polarized waves that can be transmitted without 
change of polarization correspond to vibrations with loci 
x2 Fi x2 Ve 
et aes seeder SN Bee (2) 


The vibrations corresponding to the two elliptically polarized beams are 
therefore 
x, =acos(6+«) ; Nyon (Oa aa vees (3) 
and Hy =bcos(O-+ 8);  yo=asin(@+B), “ ...... (4) 
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where « and f are phase angles relative to the incident plane polarized beam 
given by eqn. (1). 
Let vibration (1) =A x vibration (3) +p x vibration (4) so that 


~ 


ccos ¢ cos 0 = Aacos (0+) +b cos (6+ f), 
csin¢d cos 0 =Absin (9+) +pasin (6+). 


i 


By equating coefficients of cos@ and sin@ we find for the unknowns «, f, A, » 


Aa cos «+ ub cos 8B =c cos ¢, Aasina«+pbsin B =0, 
Ab cos «+ pacos 6B =0, Absinw +pasin B =csin ¢. 


Eliminating » and f, 


\(a?—b?) cosa=accos¢, A(a?—b?)sinn=—besing, ...... (5) 
whence tan « = —(b/a) tan¢, 
(a? cos? 6+ 6? sin? ¢) 
r2 — C2 “(gaye ee o 0 © ES Pe Pei Medcin ste.e (6) 
By symmetry tan B = —(a/b) tand 
b? cos? 6+ a” sin” 
pe = Ce EEE a Pear ry (7) 3 


2.13. Intensity of reflected beam. Let the reflectivities based on amplitudes 
be r and s for (3) and (4) respectively; let the phase shifts on reflection be y and 8 — 
respectively. 

The principal reflectivities and phase shifts are given by the Fresnel equations 
r= {(n—1)?+R73/{(n +1)? +h7}; tany = —2k/(k?+n?—1) and similar equations 
forsandd5; mand k denote the principal refractive indices and principal absorption 
coefficients for the two transmitted principal elliptically polarized waves. 

The reflected wave is given by 

Xp=rdAa cos (8+a+y)+sub cos (+8 +8), 
Yp=rarb sin (6+a+y)+sua sin (0+6+5). 
The intensity of the reflected (and unanalysed beam) is therefore — 
Ip = {rAa cos («+ +) + sub cos (68 +8)}*+ {rAasin (a+) + sub sin (B +8)? 
+ {rAb sin (a + y) + sua sin (6 +8)}* + {rAb cos (a +) + sua cos (B +8)? 
= (a2 + b®)(77?2 + 52u2) + 4rshuabcos(a—-B+y—8). ss aaaaee (8) 
\? and y” are given by eqns. (6) and (7), Acos« and Asin« are given by eqn. (5); 
by symmetry we also have 
u(a?—b?) cosB = —becos¢, p(a?—b*)sinB=acsing. ...... (5 a) 
With these substitutions eqn. (8) becomes i 
(@— oy 5 = (g2 1b2\p2ateost 252 sin? 252 cos? Scant 
et ae n= (a? +6?)(7?a? cos’ 6 + 1°b? sin? ¢ + s*b? cos? f + s2a? sin? ) 
+ 4rsab{ — ab cos (y — 8) + (a — 8?) sin ¢ cos ¢ sin (y — 8)} 
= 3(7° +3°)(a? + b?)? — 4rsa®b? cos (y — 8) + 4(r? — s”)(a* — b+) cos 26 
+ 2rsab(a® — b*) sin (y —8) sin 24. 
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Let the ellipticity b/a of the principal vibrations corresponding to the section 
Sbee. Then 


Vis) oe 
aes 


= = Oa — 495 008 (y- 5) + 5 cos 2h 


+2rs 


Gs 8 

itesin (y= SSID deg = ee ae (9) 
As the plane of polarization is rotated through 360°, the intensity I, will pass 

through two maxima and two minima. These points are given by 


4rs 4Ke ; 
tan 24) = a3 <= sin (y—4) = (KP 11+) sin(y—5), ...... (10) 
where K is the ratio of principal reflectivities r/s. 
Hence the directions of maximum and minimum reflectivity do not in general 
lie along the principal directions of the section. These directions coincide if 
the ellipticity e is zero as for uniaxial crystals, or if the phase difference y—8 


!} happens to be zero. 


2.14. Intensity of analysed beam (general case). Let ss be the azimuth of the 


% analyser (fig. 2). The transmitted wave is 


X, =Xpcosp+ yr sin ys 
= {rAa cos (0+%+y) +spb cos (6+ 8 +8)} cos x 
+ {rrb sin (0 +a+y)+spa sin (0 +8+5)}sin& 
=[{rAa cos («+ y) +sub cos (8 +5)} cos + {rAb sin (« +) 
+ suasin (6 +6)}sin %] cos @+[— {rdasin(«+y) 
+ sub sin (8 +5)} cos ¢% + {rAb cos («+ y) 
+spacos(8+5)}sin x] sin 0. 
_ The intensity of the transmitted beam is therefore 


I=[{rdacos (a+) +spb cos (6 +8)}? + {rAasin («+ y) + sub sin (8 +5)}?] cos? ys 
+[{rAb sin («+ y) +suasin (B +6)} + {rAb cos (~ +) 
+spuacos(B +6)}*] sin? % + 2[ {rAa cos (« +7) 

+ sub cos (6 +8) }{rAb sin («+ y) + suasin (6 +4)} 
— frdasin («+ y) + sub sin (B +8)}{rAb cos (a + y) 
_+spacos (8 + 8)}] sin x cos 
= (r?2a? + 52.267) cos? ys + (7?A2b? + 52a) sin? ys 
— 2rsAp(a? — b) sin (« — B +y—8) sin ys cos fs + 2rsAwab cos (% —B +y—8). 
Substituting from eqns. (6), (7), (5) and (5a), we obtain 


(a? — b?)?I/c? = {(r?a* + s*b*) cos® ¢ + (r? + s*)a*b? sin? h} cos? yf + {(7? + s*)a*b? cos? ¢ 
+ (r?b4 + s?a*) sin? f} sin” ys + 2rs(a? — b?){(a? — b”) cos ¢ sin ¢ cos (y — 9) 
+ ab sin (y —5)} sin ¢ cos bs + 2rsab{ — ab cos (y —8) 
+ (a? — 6) cos ¢ sin ¢ sin (y —4)}. 
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Polarizer and crystal fixed; analyser rotated. In this case ¢ is constant and 
ib varies; expressing the intensity as a sum of harmonics in ¥, we obtain 
(a? — b?)?I/c? = 4 (a? + b?){(r2a? + s?b®) cos? b + (7b? + 82a”) sin? d} 
+ 2rsab{(a? — b?) cos ¢ sin ¢ sin (y — 8) — ab cos (y —8)} 
+ 4(a® — b®)[{(r2a? — sb?) cos? 6 + (72b? — s?a”) sin? 6}” 
+ 47?s?{(a2 — b?) cos d sin d cos (y — 8) + ab sin (y — 8)}*}1 cos 2( — fo), 


Pee tats (11) 
\ 
where the phase angle ¢ is given by 
Pope 2rs{(a2 — b?) cos ¢ sin ¢ cos (y—5) + ab sin (y —8)} 
: (72a — s*b) cos? ¢ + (77b? — s?a”) sin? 6 ; 
tpe = tae Oe 
es ea) quem UL a nee (112) 


(K2—e?) cos? ¢ + (Ke? —1) sin? d 
Maximum illumination occurs for analyser setting % and %)+7; minimum 
illumination occurs at %y+7/2 and %)+37/2. ‘Thus two maxima occur per 
revolution of the analyser. 

When ¢=0, i.e. when polarizer is set parallel to principal direction 
tan 2) =2Ke/(K?—e?) sin(y—5). Hence the maximum illumination does not, 
in the general case, occur when the polarizer and analyser are parallel toa principal 
direction of the section. 

Polarizer and analyser fixed; crystal rotated. If the polarizer and analyser 
are set at angle c, 4=¢+«. In this case both ¢ and ¢ vary but « remains constant. 
From eqn. (11) 

(a® — b*)?I/c? = 4(a? + b?)?(7? + s*) — 2rsa*b? cos (y — 8) 
+ 4(a* — b*)(r? — s?)(cos 26 + cos 25) + rsab(a? — b?) sin (y — 8)(sin 2¢ + sin 2:5) - 
+4(a? — b?)?(7? + s”) cos 2¢ cos 2% + $rs(a? — b?)? cos (y — 8) sin 2¢ sin 2b 
= 4(a? + b?)?(r? + s”) — 2rsa*b? cos (y — 8) + (a? — b?)? {72 + 8? + 2rs cos (y — 8)} cos Ze 
+ 4(a* — b*)(r? — s?) cos € cos (26 + €) + 27sab(a? — b?) sin (y —8) cose sin (24 +€) 
+ 4(a* — b*)*{7? + 5? — 2rs cos ana d)} cos 2(2¢ + €). 


aL 2 

5 =}(r? +5?) a —2rs or cos (y —8) + }{r? +s? +2rs cos (y —8)} cos 2e 
+Hh(r 2 ay! — 5 cose cos (2¢ +¢) +2rs i 5 Sin (y — 8) cose sin (24 +) 
+h? A es cos (y= 0)} cos 2(20-+ 6). = ee (12) 


The intensity of illumination contains second and fourth harmonics; four 
maxima are therefore possible per revolution of the stage. The fourth harmonic 
I, is a function of the reflectivities and phase difference, the second harmonic J, 
is a function, in addition, of the ellipticity, and the relative setting of polarizer 
and analyser. 


I, 9° A (1 + e)2 . e2 psi 1/2 
2 = 4cos € {( = Ss) (l - 2p + 16r2s2 aie sin? (y nee 5)} 
: = ${7? +s? —2rscos(y—8)}, 

whence 

f, A apeee K? —1)°(1 + e?)?/(1 — e?)? + 16.K%e?/(1 — e?)? sin? (-y — 8) }¥/2 


io aI [ose or a ey 


Reflection of Polarized Light by Absorbing Materials 609 


-| 2.15. Analyser and polarizer in crossed position. When «=7/2 the second 


harmonic in the expression for the intensity vanishes and 


5 =ir+) ‘ ene -1 _ bs {1 + TH cos (y—8) 
— Hr? +s" 2rs cos (y—8)} cos 44. 


This expression is a minimum for ¢=0, 7/2, 7, and 37/2 and a maximum for 
$=7/4, 37/4, 57/4 and 77/4. The minimum intensity is given by 


TO even 88 : {r? +s" —2rs cos (y—5)}. 


lee? 
‘Therefore when a given multi-grained specimen is rotated between crossed 
polarizers every grain will show four equal maxima and minima intensity per 
complete rotation. ; 


2.16. Analyser rotated from crossed position. As the analyser is rotated from 
the crossed position the second harmonic term appears in the intensity and 
increases with the deviation from the crossed position. One pair of peaks is 
increased and the other pair decreased until, under certain conditions, the 
second pair of peaks disappears. ‘These conditions are worked out below for the 
special case when the ellipticity of the principal waves is zero. 


2.2. Uniaxial Crystals 

2.21. Uniaxial crystals. For the general case of uniaxial crystals, those 
belonging to the tetragonal, hexagonal and trigonal systems, the ellipticity e is 
zero, i.e. the principal waves reduce to plane polarized waves in the two principal 
directions, corresponding to the given case of normal incidence. 

The same is also true for special cases of non-uniaxial crystals; the normal 
' to the surface must then lie in a certain conical surface, as explained in a previous 
publication (Mott and Haines 1951a). To obtain the appropriate expressions 
for these cases e is put equal to zero in all the above results. 

Thus the expression (9) for the intensity of the analysed reflected beam 
becomes J,/c? =4(7?+5s7)+4(r? =k] cos2¢, whence Ip(0)/Ip(m/2) =1?/s? = K?. 


® Equation (10) shows that the maximum and minimum values of J, for uniaxial 


crystals occur along the principal directions of the section. 
The variation of intensity with rotation of the analyser is obtained by putting 
5=0 in eqn. (11) and e=0 in (11a): 
I/c? =}3(r? cos? d + s* sin? ¢) + 4{(r? cos? 6 —s* sin? f)? 


+4r?s? cos? (y — 6) cos? d sin? fd}? cos 2(b— iy) wee ee (14) 
where the phase angle ys is given by 
__2rs cos ¢ sin d 
ee cae Sok 
2K t 
that is tan 2p) = a = omy Ody toy Be Pla (14a) 


The analyser gives the maximum i ape when f= and P=Y%y+7. 
Minimum illumination occurs for analyser settings p=y)+7/2 and p=%)+37/2. 
Tf the polarizer is set parallel to a principal axis 6=0, and, therefore, %=0. 
For uniaxial crystals, therefore, maximum illumination occurs when the polarizers 
are parallel to a principal direction of the section. 
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The variation of intensity with the rotation of the stage is obtained by putting | 


e=( in eqns. (12) and (13): 
I/c2 =4(r? +8?) + 4{7? + s? + 2rs cos (y —8)} cos 2e + 4(7? — 5) cos € cos (26 + €) 
+ ${r? + s? —2rs cos (y —8)} cos 2(2¢ + €). 


The intensity of illumination, as in the general case, contains both second and 
fourth harmonics; the ratio of the two harmonics is obtained by putting e=0 
in (13): 


Ip 4 s virete RES Se 
. lL (ReDim Soy 


In the crossed position «=7/2; the intensity is then given by 
I/c? = }{r* + s* — 2rs cos (y — 8) }(1 — cos 4¢). 


The minimum intensity for uniaxial crystals with crossed polarizers is therefore 


Zero. 


2.22. Analyser rotated from the crossed ERE Putting 24+¢=.4, the 
equation for the intensity becomes 


I=1, +1, cosx+I, cos 2x, 
2 
= = —(J,+ 4], cos x) sin x, Ss = —I,cosx—4I, (2 cos? x—1). 

The turning points in the intensity occur when either (a) sinx,,=0 or 
(b) cosx,= —I,/4,. The first pair of turning points (a) always exist, the second 
pair (b) exist if |I,/47,| <1; if |J,/4J,| =1, the two pairs coincide. The first 
pair of turning points occurs at x,,=0or7; cosx,= +1; cos2x,,=1. Therefore 
(471 /dx"),= F1,—41,=41,{ + |1,/41,| — 1}. 

I,=4c?{r? +s —2rs cos(y—8)}>0. Therefore (d?I/dx*), is negative for both 
signs if |J,/4I,] <1, and is positive for one sign if |J,/4J,| >1. At the second 
pair of turning points 


DIN tela” Ue = I, 
(Fa), - SE, Ads (5 - 1) =41, (1 - isa) >oit| a |, 


4], 


that is, if the second pair exists. Hence if |1,/4J,| <1 there are two maxima in | 


the range 0<x<27; if |J,/4J,| >1 there is only one maximum in the same 


range. As x=2h+e, the number of maxima in the range 0<¢d<2z are four 


and two respectively. 
With eqn. (15), the condition for the existence of four maxima becomes 


K?-1 
(K—1)?+4K sin? $(y—53) 
This condition will be satisfied for all values of cif K?— 1<(K—1)?+4Ksin? 3(y—9), 


[cos « lx 


Zi aacunl geo (16) 


ie. if 2K {1 —2 sin® }(y —8)} <2 if Kcos(y—8)<1. If, therefore, K, the ratio of 


the principal reflectivities, is less than sec(y—8), four maxima will occur in the 


intensity—stage curve for all settings of the polarizers. If K>sec(y—8), two of 
the maxima existing for crossed polarizers will decrease as |cose| increases, 
and will disappear by coalescing with neighbouring minima at a setting obtained | 


by treating inequality (16) as an equation. At this setting the intensity curve 
shows a very flat minimum, due to the vanishing of the first three derivatives. 
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An upper limit to the differential phase shift y—8 can be obtained from 
eqn. (16). If the analyser setting 47—e’ just suppresses one pair of maxima, 
sine’ ={K?+1—2K cos (y—8)}/(K?—1) and 

cos (y—8) =4(1 —sine’)K+(1+sine’)/2K. 
Considered as a function of K, cos (y—5) is a minimum for 
K?=(1+sine’)/(1—sine’). 
The minimum value is given by cos (y—56) =(1 + sine’)/K =(1—sin?«’)"? =cose’. 
Hence the differential phase shift y—8 cannot exceed the critical offset ¢’. 


2.23. Angle through which the plane of polarization is rotated. If the differential 
phase shift y—6 in the case of a uniaxial crystal is small, the reflected beam is 
approximately plane polarized; this plane of polarization can be taken as 
coinciding with the analyser setting for maximum intensity as given by eqn. (14a). 

Assuming that y—6 is small, this becomes tan 2%) =2K tan d/(K?—tan? 4), 
i.e. tanyy=(tan¢)/K. The rotation w of the plane of polarization is d—y, so 
that tanw =(K—1)tan¢/(K +tan?¢), which agrees with Capdecomme and Orcel 
(1941). The maximum value of this rotation occurs when tan?¢=K, and the 
maximum rotation of the plane of polarization is given by tan w,,,, =(K —1)/2+/K. 

If y—6 is small, more precisely if sind(y—5)<(K—1)/2./K, the analyser 
setting at which one pair of intensity maxima disappear is given by 

| cose | =(K—1)?/(K?-1)=(K—1)/(K+1). 
If <’ is the departure from the crossed position «=47—e’ then 
sine’ =(K—1)/(K +1). 
Co RA LY Kee ah es 
Grok Ke 1 Ka, las 
Thus, provided that y—5<K-—1<1, then, as the analyser is rotated from the 
crossed position, the intermediate maxima vanish for an offset equal to the 


maximum angle of rotation of the plane of polarization, a result found by 
Sampson (1929). 


Further if K-—1<1 


§3. EXPERIMENTAL OBSERVATIONS 
3.1. Procedure 


Specimens of tin, uranium, bismuth and stibnite were chosen for the tests 
and were prepared as follows: 


Tin (tetragonal). Deformed by compression, mechanically polished and slowly 
cooled after recrystallization at 100°c to reduce the incidence of twinning. The 
surface was finally prepared by electrolytic polishing in 60% perchloric acid as 
recommended by Puttick (1949). 


Uranium (orthorhombic). Cold-worked and recrystallized in the upper limit 
of the alpha-range. The surface was given a rough mechanical polish before 
electrolytic polishing in a solution containing 2 volumes sulphuric acid, 1 volume 
orthophosphoric acid and 2 volumes distilled water as described by Mott and 


- _ Haines (1951 b). 


Bismuth (rhombohedral). Deformed and _ recrystallized. | Mechanically 
polished prior to 4 minutes electrolytic polishing in the solution described for 
uranium. 


25-2 
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Stibnite (orthorhombic). Polished mechanically on a lead lap and fir 
on soft felt using a suspension of rouge in oil. 

The vertical illuminator on the microscope was a zinc sulphide coated glass 
slip and both the polarizer and analyser units were of polaroid sheet, the setting 
of the analyser being controlled and measured to an accuracy of about +0-05°. 
A Wratten C filter (green) was used in conjunction with a mercury vapour lamp 
which was not stabilized against mains variations but showed only slight changes 
in intensity throughout a given experiment. The intensity of the reflected beam 
was measured by a photomultiplier tube of a microphotometer unit (AC Type 
1175A constructed inthe Electronics Division, A.E.R.E.) insertedin the monocular ~ 
eyepiece tube of the microscope. 

The microscope was set up by centring all objectives with respect to the 
rotating stage, which had been centred by the makers with respect to the optical 
axis. Using a highly polished, stainless steel specimen with the analyser, eyepiece 
and objective lens removed, the polarizer was rotated to give maximum intensity. 
This ensured that the plane of polarization was in the plane of incidence at the 
vertical reflector. The analyser was then rotated to give a new minimum with 
the lamp iris and field iris closed to minimum aperture, so that the specimen _ 
was illuminated by a narrow cone at near-normal incidence. The objective lens 
was then inserted and rotated in a special mount (which ensured that true centring 
was maintained for all settings) to a position in which minimum intensity was _ 
obtained. In this position the birefringence due to the strain in the objective 
was compensated as completely as possible by the stray polarization effects in 
the rest of the optical system, i.e. those due to the mounting plates used for the 
polaroids and the reflection at the glass slip illuminator. Under these conditions 
the four maximum intensities obtained on rotation of a single grain of an aniso- 
tropic metal under crossed polarizers should be equal. This was checked for 
each specimen and slight adjustments in the setting of the objective made, if 
necessary, to achieve this result. 

A grain about 0-1—0-2 mm diameter and free from twinning was chosen in 
each specimen in turn. The iris diaphragm controlling the field aperture was 
closed down so that the image of the iris was totally enclosed within the chosen 
grain. In the first set of experiments intensity measurements were made as 
the stage was rotated through intervals of 10° except in the region of the maxima 
and minima intensities when some readings were taken every 5°. This was 
carried out for the crossed position and with the analyser rotated through various 
small angles up toa maximum of 4°. This gave experimental curves for the change 
in the intensity with azimuth angle ¢ for various values of e’ in the region of the 
crossed position. In addition, for some specimens, the position of the analyser 
was found at which the intensity was a minimum for various stage settings. 
From this the relationship between w, the angle of rotation of the plane of 
polarization, and the azimuth angle ¢ was obtained. 


3.2. Results 


Tin. The curves shown in figs. 3 (a) and 3(5) were obtained for two grains 
in the tin specimen differing markedly in orientation. From fig. 3 (a) it will be 
observed that for the crossed position the intensity—azimuth angle curve shows 
four peaks occurring at intensity values of 3-3, 3-4, 3-3 and 3-4 respectively. 
Alternate pairs of peaks are slightly lower than the others, which suggests that 
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_ the system as a whole was not exactly equivalent to the crossed polarizer condition ; 
_ this may have been due to a slight inaccuracy in setting the analyser or to incomplete 
compensation of the stray polarization effects in the microscope. ‘The difference 
in alternative peaks in fig. 3(a) is more marked in the curve obtained after the 

“4 


KEY 


Analyser Rotation 
from Crossed Position 


20 --K---- AL? 


Intensity 


8 xr : 
x : x 
a) x 
Sn . , \ 
bik 7 
4 Lg a Me ae ! a 
oN, ¢ 
7 2 . ee Le x = f 
ao, J \ vo el) LoS 
No IN: he 


0 40 80 120 160 200 240 280 320 360 
Stage Rotation (degrees) 


Fig. 3 (a). Electrolytically polished tin. 
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Fig. 3 (6). Tin electrolytically polished in perchloric acid solution. 
Note. Upper scale should read -—0-8 to 0-8 degrees, not —8 to 8. 


analyser has been rotated by 0-5° from the crossed position, the ratio between 
the pairs being about 2-5: 1. After an analyser rotation of 1°, the intermediate 
peaks are reduced still further and are only just discernibie. With further 
rotation of the analyser to 2° and 3° respectively from the crossed position, only 
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two maxima are obtained and the minima have become more sharp, as predicted — 
in §2. 

Since the alternate maxima are almost eliminated at an analyser rotation — 
of 1°, the probable maximum angle of rotation of the plane of polarization for 
this particular grain is of the order of 1-0°, which agrees with the value obtained 
from the (w, ¢) curve (not reproduced). 

It will be observed that, although the minimum intensities for the crossed 
position are equal, they have a positive value instead of being zero as predicted 
by theory. This is due to low intensity scattered light in the laboratory reflected 
by the specimen, and also to light scattered by imperfections in the polished 
surface. Further, the increase in the minima obtained by increasing the rotation 
of the analyser is evidence of the phase difference between the two principal — 
reflected beams. 
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Fig. 4. Uranium polished in standard electrolyte. 


Both the (w, ¢) and the (intensity, 4) curves are reproduced in fig. 3 (0) for 
the second grain. The (intensity, 6) curves show similar characteristics to those 
in fig. 3(a), and four peaks are obtained after an analyser rotation of 0-5° and 
two after 1°. [he maximum rotation angle w obtained from the (@, 4) curve is 
from 0-7-0-8°. 

Uranium. ‘The results given in fig. 4 show that, for the grain chosen, four 
peaks are just apparent after an analyser rotation of 1°, but only two are obtained 
after 15°. ‘The maximum angle of rotation w as obtained from the (w, d) curve 
is about 0:5-1-0°. 

Bismuth. 'Vhe results are reproduced in fig. 5 and indicate a larger rotation — 
angle for the chosen bismuth grain than that given by the grains of tin and 
uranium. Four peaks were obtained after an analyser rotation of 3°, but two 
flat minima were evident after 4°. This agrees with the measured maximum 
angle of rotation w of 3-6°. 

Stibmte. 'Uhe rotation angle for stibnite is usually greater than that for metals, 
and this is illustrated by the curves shown in fig. 6. Four maxima and minima 
were obtained after an analyser rotation of 4° from the crossed position but only 
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two after a 6° rotation. The (w,¢) curve is not reproduced but gave a maximum 
rotation angle w for the grain investigated of about 5-2°. It is interesting to note 
that, as in the case of the metals, the minimum intensities increase slightly on 
rotation of the analyser, indicating different phase changes for the two principal 
_components on reflection. 
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Fig. 5. Bismuth. : Fig. 6. Stibnite. 


§4. DISCUSSION 

It has been shown both mathematically and experimentally that on reflection 
of polarized light at normal incidence from anisotropic surfaces, when the analyser 
is exactly crossed with the polarizer and there are no stray polarization effects 
due to the instrument, four intensity maxima and minima are obtained for a 
change in azimuth angle of 27. This condition appertains to all materials, 
irrespective of the crystal symmetry and of the absorption characteristics. ‘The 
effect of rotation of the analyser from the crossed position, however, does depend 
on the nature of the reflecting surface. For those materials for which the 
differential phase shift is small, alternate maxima are gradually reduced until 
eventually a setting of the analyser is reached in which only two maxima per 
complete stage rotation are observed. If there is no differential phase change 
between the two principal components on reflection at a uniaxial crystal, the 
minimum intensities should be zero for all settings of the analyser up to that 
which just eliminates the alternate maxima; if the differential phase change is 
small, the minimum intensities will increase slightly with offset. Ifthe differential 
phase shift is sufficiently large, then four maxima and minima will be obtained 
for all settings of the analyser. 

The above results apply only to optical systems which are free from stray 
polarization effects, and, in general, ideal conditions are not obtained. ‘The 
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main sources of uneven polarization over the field of view are strain in the objective 
and polarizer mounts and refraction and reflection at oblique incidence to the 
vertical reflector unit. In the present work the setting of the objective was 
such that the elliptical effects due to strain partially corrected for the non-uniform 
polarization resulting from the slip illuminator. The precise effect on the intensity _ 
analysis of the combination of errors due to defects in the optical system will be — 
discussed in a later paper. It would appear, however, that many of the anomalies 
reported in the literature can be adequately explained on the basis of the stray 
polarization effects in the optical systems used. In the case of etched metals 
the ellipticity of the reflected beam is considerably more marked than for polished 
surfaces for reasons discussed previously (Mott and Haines 1951a), so that 
the results are only slightly affected by stray polarization. 
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The Magneto-Resistance of Ferromagnetic Al-Si-Fe Alloys 
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ABSTRACT. By extending the treatment in an earlier paper an equation is derived for the 
saturation magneto-resistance of some mixed ferromagnetic alloys as a function of composi- 
tion and temperature. Results calculated from this equation are in good agreement with 
experiment for the Al-Si-Fe alloy system. 


§1. INTRODUCTION 
N a recent publication (Parker 1951 a, referred to as I) the author has suggested 
| that by making a number of plausible assumptions an equation can be 
derived for the saturation magneto-resistance coefficient (Ap/p),4 of certain 
ferromagnetic binary alloys; (Ap/p),,, may be defined as the relative change in 
electrical resistivity when an initially demagnetized material is magnetized to 
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technical saturation. The equation derived describes (Ap/p) sat as a function 
of temperature and composition of the alloy. There is good quantitative 
agreement between the calculated value of (Ap/p),;,4. which is the value of 
(Ap/p)sa, When the magnetic field is applied along the direction of current flow 
in a specimen, and the experimental values for alloys of iron with silicon 
(Shirakawa 1939) and iron with aluminium (Parker 1951 b, referred to as II). 
It is shown in § 2 below that the treatment in I may be extended to the case 
of alloys containing more than two elements. In §3 some calculated values of 
(Ap/p) sa, for dilute solid solutions of aluminium and silicon in iron will be 
compared with experimental values obtained by the author, and some interesting 
magneto-resistance features of this alloy system will be discussed. 


§2. THE SATURATION MAGNETO-RESISTANCE COEFFICIENTS OF 
SOME MIXED ALLOYS 

As already noted in I, the resistivity p of a metal containing a dilute solid 
solution of other elements in its crystal lattice may be expressed by p=pp+ps; 
py tepresents the resistivity of the pure metal and pg the contribution to p of 
the solute atoms. For many such alloys it may be stated as a good approximation 
that over a wide temperature range py increases almost linearly with the 
temperature ¢, while pg is independent of ¢ (Matthiessen’s rule). 

According to the quantum theory (cf. Mott and Jones 1936), py» at ordinary 
temperatures arises from the scattering of conduction electrons by the thermal 
vibrations of the ionic lattice of the metal; pg is due to scattering of electrons 
by the aperiodicities in the electrostatic field in the lattice resulting from the 
presence of the solute atoms. For sufficiently dilute solutions pg is usually a 
linear function of the solute concentration. 

The above statements can be summarized in the following way. If Ny and N,; 
are the numbers of atoms per cm? of solvent metal and of the 7th solute element 


| respectively in the lattice, then 


p=Nyoyt+d,N,o; and Ny>>,N,, ene (1) 


where o, and o; are measures of the efficiency of scattering of conduction 


* electrons, averaged over all crystal directions, due to thermal vibrations of the 


lattice and to solute atoms of the ith element respectively. As in I, the 
contribution to p of the thermal vibrations of the solute atoms is neglected. 
Let A again denote the changes when a specimen is brought from the 
demagnetized state to technical saturation; then, by definition, 


(Ap/p) sat =(NpAcy + U,N,Ag;)/(Npop + U,Nio;). 
This equation reduces to 


(Ap/P) sat =[Pr(Apr/pr) + &,pi(Api/pi)]/Ler + psl- 
The last equation can now be applied to the case of dilute solutions of aluminium 
and silicon in iron. Let the suffixes 1 and 2 refer to the elements aluminium 
and silicon respectively, then 


(Ap/p) sat =lPx(Apr/pr) + pi(Api/ps) + p2(Ape/p2)I/(pr+ ps): ++ +++ (2) 
For the case of the longitudinal saturation magneto-resistance all quantities 
on the right-hand side of eqn. (2) can be obtained from the results of independent 
experiments, and contain no unknown parameters. ‘Thus direct comparison 
betwen (Ap/p)) .,4 obtained from eqn. (2) and the experimental value can be made. 
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In the calculated curves of §3 below, the value of py is the mean quoted for 
iron by Cleaves and Thompson (1935), and Apy/py is obtained from the 
experimental results of Shirakawa (1939, 1940); Ap,/p, and Ap,/p. are taken from 
the results on the binary alloys reported in I and II. Values of pg were obtained 


from resistivity measurements, and it was ascertained that pg of each alloy q 
considered is approximately independent of temperature. The only novel — 


assumption here arises from the estimate of p, and p,. According to eqn. (1) 
solutions of two separate elements contribute independently to py. If this 
equation applies to the aluminium-silicon—iron system it is possible to write 


PS§=Pi + Po coeoeneu 
The results of a very extensive investigation at room temperature by Masumoto 
(1936) suggest that eqn. (3) is closely obeyed for solute concentrations not — 
exceeding 7% by weight. Resistivity measurements on the specimens discussed 
in §3 are in good agreement with Masumoto’s results. The values of p, and p, 
in eqn. (2) have therefore in each case been taken from the well-known results 
for binary alloys. 


§3. EXPERIMENTAL RESULTS AND DISCUSSION 
Measurements of longitudinal magneto-resistance in fields of up to 
1200 oersteds and in the temperature range —78°c to +300°cC were made on _ 
three alloys. ‘The chemical analysis, in per cent by weight, of the alloys is 
given in the table: 


Alloy Al Si Cc P Mn Ss Cu Ni 
A 0-45 0:46 0:01 0:02 0-12 0-02 0-09 0-06 
B 0-44 1:24 0:01 0:02 0:16 0-01 0-08 0-05 
c 0:96 0-81 0-01 0-03 0-13 0-02 0-08 0-04 


Values of (Ap/p)) a, Were estimated in the usual manner. ‘The results, 
together with the calculated curves, are shown in fig. 1. ‘The agreement obtained 
is about the same as that reported in I and II and is well within the limits of the 
experimental errors in the determination of the quantities used in (2). 

The transverse saturation magneto-resistance coefficient (Ap/p), 4, of the 
alloys for magnetization perpendicular to the direction of current flow was also 
measured between 15°c and 300°c, the results being recorded in fig. 2. It was 
pointed out in II that direct comparison with eqn. (2) is not possible owing to 
lack of data of (Ap/p), 44, for pure iron. A short calculation, however, shows that 
(analogous to the results for aluminium—iron alloys) the variation of (Ap/p), sat 
in the above temperature range should not exceed 0-3 x 10-3 for any of the above 
alloys. It can be seen from fig. 2 that this is the case. 

So far discussion has been confined to the case of magneto-resistance 

_coeficients associated with technical saturation. Unusual features in the 
longitudinal magneto-resistance coefficient (Ap/p), below technical saturation 
were observed in the cases of alloys B and C. For most ferromagnetic materials 
(Ap/p) increases or decreases steadily with H until saturation is reached if the 
usual experimental errors are eliminated (Stierstadt 1930). For the above two 
alloys, however, it was found that the [(Ap/p),, H] curves exhibit a marked initial 
minimum, which becomes less pronounced with increasing temperature. ‘Typical 
curves are shown in fig. 3. Somewhat similar characteristics have previously 
been reported by Shirakawa (1939) for 1-66% Si-Fe, and by Parker (1950) for 
1-8%, Si-Fe and for a single crystal of 2:8% Si—Fe. 


The Magneto-Resistance of Ferromagnetic Al-Si-Fe Alloys 619 


e Experimental Points 


~100 0 100 200 300 
Temperature (°c) ° 


Fig. 1. (Ap/p)y sa¢ plotted against temperature. 


(a) Alloy A: (6) Alloy B: (c) Alloy C. 
Curves are calculated from eqn. (2). 
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Fig. 2. (Ap/p), sat Plotted against temperature. 
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Fig. 3. The longitudinal magneto-resistance of 0:96% Al-Fe, 0°81% Si-Fe. 


(a) 18°c, (b) 242°c. 
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The writer suggests the following tentative explanation of the above 
phenomenon. It has already been pointed out that og and oy make respectively — 
negative and positive contributions to the value of (Ap/p)) .a-. In the early part 
of the magnetization process the temperature independent scattering makes a 
larger contribution to the magneto-resistance coefficient of the alloy, but, as 
technical saturation is approached, the contribution from the thermal vibrations 
of the iron lattice becomes the dominant term. 

This accounts quite simply for the rare occurrence of the initial minimum 
in the [(Ap/p),, H] curves in ferromagnetic metals and alloys. Conditions for 
this occurrence are favourable only in the case of alloys where Ap,/pp and 
Aps/ps are of opposite sign and where py and pg are of approximately the same — 
magnitude, i.e. in 1-3°% Al-Fe or Si-Fe alloys at room temperature. 

The observed variation of the minimum for a given alloy with temperature 
may be predicted qualitatively. As the temperature is raised, the contribution 
of py becomes relatively more important for a given state of magnetization. 
Consequently the minimum should become less pronounced and occur in lower 
applied fields at high temperatures. ‘These predictions have been verified for 
both alloys A and B. A quantitative treatment of this problem cannot be given 
without an extensive experimental study of the magneto-resistance properties 
and the magneto-crystalline energy of single crystals as a function of temperature ~ 
and composition. 
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Sparking Potentials in a Transverse Magnetic Field 


By J. M. SOMERVILLE 
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ABSTRACT. A theory of sparking between parallel plates in a transverse magnetic field, 
put forward by Valle in 1950, is discussed and shown to be in qualitative disagreement with 
observation. By making allowance for the distribution of electron free paths, and for electron 
recapture by the cathode, a modified theory is produced which is in better agreement with 
observation. Further refinement of the theory requires an approximate expression for the 
first Townsend coefficient which is valid when a magnetic field is present. The possibility 
of obtaining such an expression is discussed. 
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§1. INTRODUCTION 
ALLE (1950) has recently published a simplified theory of sparking between 
\ / infinite parallel plane electrodes in the presence of a uniform transverse 
magnetic field parallel to the electrodes. It is the object of this paper to 
outline a theoretical treatment which takes into account a number of factors 
ignored by Valle and which is in better agreement with observation. 
Valle’s theory is based on the equations 


Pesto Aa tia ye Tem Tile ee BE love reer (1) 
and peep (pie at) ae Ve) De (2) 
where « is the first Townsend coefficient of ionization by electron collision, y is the 
probability of a positive ion releasing an electron by impact at the cathode, dis the 
distance between the electrodes (in centimetres), F the electric intensity (in volts/cm) 
and A and B are empirically determined constants. 

Equation (1) is the form taken by the well-known Townsend condition for 
sparking when electron emission from the cathode by positive ion bombardment 
is the chief secondary source of electrons. Equation (2) is the approximate 
expression for «/p as a function of E/p derived by Townsend. 

To allow for the presence of a transverse magnetic field Valle replaces the 
pressure p in eqn. (2) by an ‘equivalent pressure’ p, derived as follows. 


Fig. 1. 


Assume that the electric field E between the electrodes is uniform and that 
electrons lose all their energy at collision. ‘Then the path of an electron between 
collisions will be anarcof acycloid (fig. 1). The lengthof onecomplete arch of this 
cycloid we may call a,* where 


eee Vel? ces the) 9 ilk ai ete eS (3) 
Hi being the magnetic intensity in oersteds and e and m the charge (e.m.u.) and 
the mass (g) of an electron. 


If Sis the length OP measured along the arc of such a cycloid and / the distance 
of P from the base, then 


INVER IO Sean ean ern (4) 
if P lies on the first arch or, more generally, 
h={S—(m—l)a}[1-{S—(m—l)a}/a] isn (5) 


if P lies on the mth arch. 

Valle assumes that all the electrons have the same free path \ and confines his 
attention to the first arch of the cycloid so that, if hy denotes the distance travelled 
in the direction of the electric field by an electron whose free path is A, then 


hy =A(1 —A/a). ee os ee (6) 


* Valle denotes the length of a complete arch of the cycloid by 8a. Throughout this paper 
Valle’s results are expressed in the notation most suitable to our purpose. 
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He then takes as the ‘equivalent pressure’ 


Peep Ne (7) 


= Like, ‘einer (8) 
where Lepr; Set ged» eat tora eee (9) 


and is thus equal in magnitude to the mean free path at unit pressure. The 
equivalent pressure p, is inserted in eqn. (2) instead of the actual pressure /, giving 


ahy=AL exp(—BL/Ehy). snes (10) 


If V is the potential difference between the electrodes expressed in volts then 
E=V/d and (10) may be written 


oh AL exp (= BLG/V hg)... bons nee (11) 
Elimination of « between (1) and (11) gives 
hy ALexp(—BLd/Vhy) 


te log (1+ 1/y) 
For a given gas and cathode material, i.e. for given A, B, L and y, eqns. (3), (6), (9) 
and (12) determine the sparking potential V as a function of pdand Hd. Valle has 
devised a convenient scheme for the computation of V from these equations and 
has carried out the calculations for the particular case when the gas is air and the 
cathode materialisaluminium. Hestatesthat the consequences of his theory are in 
qualitative agreement with observation. There are, however, several major 
qualitative differences between theory and observation which make it difficult to 
agree with this statement. Some of these are as follows : 

(i) The theory predicts (cf. Valle, fig. 4 (0), (c), (d)) thatina given magnetic field 
the same value of sparking potential may, in general, be achieved at four different — 
pressures. ‘There does not appear to be any observational support for this, and 
observers like Meyer (1919), Wehrli (1922) and Penning (1937) agree in finding, 
for a given magnetic field, only two pressures corresponding to each value of the 
sparking potential. 

(ii) The theory predicts that under some circumstances (cf. Valle, fig. 4(6)) a 
range of pressures exists for which no sparking at all is possible at any voltage, 
although sparking may take place at pressures above or below this range. This 
again lacks observational support. 

(ii) At low values of the product pd the theory yields curves of sparking 
potential against magnetic field strength (cf. Valle, fig. 4(c)) which move to the 
left, i.e. towards lower magnetic fields, as pd decreases, and which in general 
intersect one another or else lie entirely outside one another. Figure 2 is a sketch 
showing the type of (sparking potential, magnetic field strength) curves obtained 
experimentally by Penning (1937) in argon between coaxial cylindrical electrodes. 
It is seen that the curves corresponding to different pressures do not intersect and 
that any curve corresponding to a lower pressure lies entirely within any curve 
corresponding to a higher pressure. Moreover, the curves move to the right as 
the pressure decreases. 

(iv) he theory predicts that, no matter how low the pressure may be, the 
minimum sparking potential for a given pressure may always be reduced to the 
same value by asuitable choice of magneticfield. Figure 2 shows that the observed 
minimum attainable sparking potential increases rapidly as the pressure is lowered. 
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We now seek the physical reasons for these points of disagreement between 
theory and observation with a view to the development of a theory more in line 
with experiment. In so doing it will be convenient to make one modification at a 
time and follow out its implications rather than to make several modifications at 
once. 

§2. THE DISTRIBUTION OF FREE PATHS 

The first and second points of disagreement between Valle’s theory and 
experiment are directly attributable to neglect of the distribution of free paths, the 
assumption being made that all electrons have the same free path between collisions. 
This may be seen as follows : 

For any given value of V, eqn. (12) in general determines two values of hy. 
Each of these values of hy determines, through eqn. (6), two values of A (the arcs 
OQ and OR in fig. 1) and hence two values of /, giving, in all, four* values of p 
corresponding to each value of V. ~ 


Log Sparking Potential 
Pressure Decreasing 


Magnetic Intensity 


Fig. 2. Sketch (not to scale) showing types of (sparking potential, magnetic intensity) curves for 
different gas pressures obtained by Penning (1937) in argon between coaxial cylindrical 
electrodes in an axial magnetic field. 


Suppose now that we recognize that all electrons do not have the same free path, 
but that the probability of a free path length between x and x + dx is exp (— x/A)dx/A. 
We would then replace h, by A, the mean value of f at collision, giving, instead of (8), 


Lee by Saal dae al en acncaiae, BOR: (13) 
Equations (10) and (12) then become, respectively, 
ime eth CB DER) | Aerial ehh es (14) 
De Ba Vi) ee gee (15) 
d log (1+ 1/y) ; 


The value of h in terms of A and a has been given by Kugler, Ollendorf and 
Roggendorf (1933), and we may write it in the form 


P= Afeoth( seNaataliy saw 26 pn Sy dest (16) 
which is the equivalent equation replacing (6). 


* In using an equation equivalent to (6) Valle has arbitrarily limited consideration to the first 
cycloidal arch, but in general if an electron completes such an arch without collision it will immedi- 
ately begin another identical one, and other values of A (e.g. the are OPQRO’S) corresponding to 
the value iy become possible, leading to an infinite number of values of the pressure corresponding 
to each value of V. 
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This equation may also be conveniently written in the form 
hja=(A/a) {coth(a/2A)—2Aja} ie (16a) 
or hia= Fla) kis ... se chek ee an nee (165) 
where F(x) =x«{coth (1/2x) —2«}, giving hja as a function of A/a. Equation (16) — 


was derived on the assumption that the electrons were free to go on describing an 


indefinite number of cycloidal arches until collision. It is easy to see, however, 
that the mean value of h for the electrons which collide in the first arch is the same 
as the mean value for the whole group. Consider a group of N electrons starting 
from rest at O (fig. 1). This group may be divided into two subgroups, namely the 
N exp(—a/A) electrons which reach O’ without colliding, andthe N{1 —exp(—a/A)} 
electrons which collide in the first arch. Now the N exp(-—a/A) electrons reaching 
O’, and beginning a new arch from there, differ only in total number from the 
original N starting from O; their chance per unit distance of colliding is the same, 
and the probability of an electron starting from O’ reaching a certain distance 
before collision is the same as the probability of an electron reaching the same 
distance starting from O. Consequently the original N electrons and the subgroup 
of N exp(—a/A) have the same mean value of / at collision. ‘Therefore the re- 
maining subgroup, i.e. the electrons colliding in the first arch, have this same mean 
value. (This argument can clearly be extended to show that the electrons colliding 
in the second or any other arch have the same mean value of h at collision.) In | 
applying eqn. (16) it is therefore not necessary to distinguish between electrons 
generated in the gas, which are free to describe an unlimited number of cycloidal 
arches before they collide, and electrons generated at the cathode, which may be 
recaptured at the end of the first arch, since the value of his the same for either class 
of electrons. be 

Returning to eqn. (16), it is not difficult to show (see Appendix) that h is a 
single-valued monotonic increasing function of A which asymptotically approaches 
the limit a/6 as A tends to oo. With eqns. (15) and (16) replacing (12) and (6) 
respectively there will be only two values of the pressure corresponding to each 
value of V, because (15) will determine two values of h for each of V, and (16) 
will determine only one value of A, and hence of p, for each value of h. This 
modification therefore removes our first objection to Valle’s theory. 

It removes the second also, as we shall now see. In eqn. (15) the exponential 
cannot exceed unity, and hence nosparking is possible if h> A Ld/log ie +1/y)=hy 
say. Since h is a monotonic function of A, and hence of p, once p is so low that 
h>h,, further reduction in p will only result in further increase in h, so that, once 
sparking becomes impossible, it remains impossible at lower pressures. 

Some (V, Hd) curves for various values of pd computed from eqns. (15) and (16) 
are shown in fig. 3. _In computing these and other sparking potential curves in this 
paper the following values of the constants have been used, appropriate to air as 
the gas and copper as the cathode material, namely A =14-6(cmmm Hg)! 
B=365 vicmmm Hg, as given by Townsend (1915), L=0-026cmmmHg, ~ 
y =0-025, as given by Engel and Steenbeck (1932). It will be seen that our third 
objection has also been partially removed in that the curves no longer move to the 
left with decreasing pressure, but the fourth objection still remains, and the 
modified theory still implies that the sparking potential can be reduced to the same 
minimum value for all values of pd. Moreover the sparking potential curves now 
tend to a limiting formas the pressure is reduced, and there is little change in the 
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curves for values of pd below about 0:01cmmmHg. This independence of 
pressure at low pressures is due to the fact, already discussed, that h tends to a 
limiting value (a/6) as A increases. 

Physically this independence of pressure is due to the assumption that, in the 
magnetic field, electrons keep on describing a series of cycloidal arches until they 
collide with gas atoms. _ If the pressure is lowered, electrons describe more arches 
before making a collision, but the collision is ultimately made, however low the 
pressure may be, and is just as effective in promoting sparking as if it had been 
made earlier. The situation would be changed if a process existed whereby 
some electrons were removed from the space between the electrodes before collision 
and if the chance of such removal increased as the pressure decreased. One such 
process of removal is by electron recapture at the cathode. 
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Fig. 3. Calculated sparking potential V for Fig. 4. Calculated sparking potentials, for the 
sparking between parallel plane copper same conditions as those of fig. 3, when 
electrodes a distance d cm apart in air allowance is made both for the distribution 
at a pressure mm Hg in the presence of electronic free paths and for electron 
of a transverse magnetic field H oersted, recapture by the cathode. 
allowance being made for the distribution At points on the dotted line the condition 
of electronic free paths. H/E (=Hd/V)=1°1 is satisfied. Its signifi- 


cance is discussed in the text. 


§3. BLECTRON RECAPTURE AT THE CATHODE 
Electrons generated at the cathode which describe a complete cycloidal arch 
without making a collision will in general be recaptured by the cathode at the end of 
thearch. Thus, of N electrons leaving the cathode, N exp(—a/A) will complete an 
arch of length a without collision and be recaptured by the cathode, and, so far as 
their influence on breakdown is concerned, they might as well never have existed. 
Consequently, of every N electrons generated at the cathode only N{1 —exp(—a/A)} 
are effective in promoting breakdown, so that the coefficient of emission of electrons 
by positive ions at the cathode is effectively reduced by the factor 1 — exp(—a/A). 
We can therefore make allowance for electron recapture at the cathode by sub- 
stituting y{1 —exp(—a/A)} for y in eqn. (15), obtaining 
lew AL exp(—BLd/Vh) 
ad log{1+1/yfl—exp(—a/A)}] 
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Pees (7) 
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This equation, together with eqns. (3) (9) and (16), now gives values of the sparking | 
potential with allowance both for the distribution of electron free paths and electron _ 
recapture at the cathode. 4 | 

Figure 4 shows curves of sparking potential plotted against Hd derived from _ 
eqn. (17). It will be noted that in these curves the desired progression of the _ 
minima upward and from left to right with decreasing pressure has been achieved, _ 
and that the qualitative agreement with Penning’s curves (fig. 2) is better. 


§4. OTHER SOURCES OF ELECTRON LOSS 

Recapture by the cathode is not the only process causing loss of electrons from 
the space between the electrodes. Any actual discharge tube with plane parallel 
electrodes must be finite, and so must be the extent of the magneticfield. Electrons 
may be lost as effective agents in promoting sparking if they diffuse along the” 
magnetic field beyond the edges of the electrodes or out of the magnetic field or to_ 
the walls of the containing vessels. Such loss would increase with decreasing 
pressure. An analysis of its effect on the sparking potential presents considerable’ 
difficulties and is not attempted here. 

For our purpose it is not necessary to consider electron loss due to electron 
motion across the magnetic field because, although we are using a plane model for 
simplicity, we have in mind the qualitative comparison of theory with experimental 
measurements, like those of Penning, of sparking potentials between coaxial — 
cylindrical electrodes in an axial magnetic field. In such systems electron loss 
perpendicular to the magnetic field does not occur. 


§5. MODIFIED EXPRESSION FOR a IN A TRANSVERSE 
MAGNETIC FIELD 


One aspect of the theory as developed so far which is particularly open to 
criticism is the replacement of the pressure ‘n Townsend’s approximate expression 
for «/p (eqn. (2)) by the ‘equivalent pressure’ p, defined by eqn. (13), giving 
eqn. (14). E 

Equation (2) is derived as follows. Assume that electrons lose all energy at 

collision, that an electron with energy greater than eW’” ionizes at collision and that 
one with less does not. In one centimetre path length an electron will make onthe _ 
average 1/X’ collisions. (‘The symbol A’ is used for the mean free path in recogni- 
tion of the fact that values of A’ derived from eqn. (2) are not equal to the kinetic 
theory mean free path or the mean free path measured in other ways, and W’ is used 
for similar reasons). Of these 1/A’ collisions a fraction exp (—//A’) will be preceded 
by a free path exceeding /. Assuming the electric field to be uniform and equal to 
E, then if we put /= W’/E, we have the number of ionizing collisions per centimetre, 


namely = (1/0) expC PEN)? te (18) 
or aip=(l/)L exp (WEL 2 ee eee (19) 
writing L’ for pA’. This leads.to eqn. (2) if we write 

ABIL B= WY Lo ee (20) 


If we try to follow the same line of argument when a transverse magnetic field is 
present, we find that we cannot simply replace the free path / by the free distance 
travelled in the direction of the electric field. This is so because values of J and h 
are not distributed in the same way, and it is not true to say that the fraction of free 
distances in the direction of the electric field is equal to exp (—h/h). 
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Asa preliminary to obtaining an equivalent equation to (2) when a magnetic 

field is present, we will now derive the distribution function for h. It is clear, from 
what has already been said, that we need only calculate the distribution function for 
electrons colliding in the first arch, since the distribution function for electrons 
colliding in other arches will be the same. Referring to fig. 1, an electron starting 
from rest at O will be at a given distance h from the base line twice, at distances from 
O which we may call S, and S,, and which are the roots of eqn. (4) : 

S,=$a{l —(1—4h/a)¥?}, S,=4haf{l+(1—-4h/a)"}. (21) 
Of N electrons starting from O, the number colliding at a distance from the base 

line between h andh+dhis Nexp(—S,/d) |dS,|/A+N exp(—S,/A) |dS, |/A, where 

_ dS, and dS, correspond to dh, and, from (21), are equal to + dh/(1—4h/a)". 

The total number of electrons colliding in the first arch is N{1 — exp(—a/A)}, so 

that the probability of collision within the interval h to h + dh is 

_ P(h)dh = [ Ndh exp (—a/2A) [exp{a(1 —4h/a)"?/2A} + exp { — a(1 —4h/a)"?/2A}] 

A(1 —4h/a)2N{1 — exp (—a/A)} 


__ cosh {(a/2A) (1 —4h/a)"?} dh (23) 
Nemenialent (i sia nd a 
Remembering that / cannot exceed the height 4a of an arch, the fraction of free 
distances in the direction of the electric field exceeding h is 
sinh {(a/2A) (1 —4h/a)¥?} 
i DA east ge ORy oo wales 22 
It is not difficult to show that this expression reduces to the form exp(—ZJ/A) as 
the magnetic field H-0, i.e. as aco and hl. Consequently, if we put A” for A 
and W’”/E for h, where X” and W” are the values appropriate to this case, we obtain 
an expression for the fraction of electrons colliding with energy in excess of eW”, 
i.e. the fraction of ionizing collisions, of the form 


sinh {(a/2A”)(1—-4W"/Ea)*\/sinh(a/2\") ns a (25) 
which reduces toexp (i W" /E") (cf. eqn.(18)) when H=0. Hence we must have, 
Bom eqn, (20), HAG Bee ee P Onein ene (26) 


To obtain « we must multiply expression (25) by the number of collisions per unit 
length in the direction of the electric field, a quantity which we may denote by 
1/hk”. Doing this, and using (26), we obtain 

a =(1/h") sinh {(a/2d") (1 —4.Bpd"/Ea)*?}/sinh(a/2A") kaa (27) 
as the equation equivalent to (18). 

We now encounter a difficulty, in that we lack criteria which will allow us to 
determine \” andh”. Just as, in eqn. (18), the identification of \’ with the kinetic 
theory mean free path, or the mean free path measured in other ways, would lead to 
values of A not in accord with observation, so, in eqn. (27), the identification of X” 
and h” with similar quantities measured or calculated in other ways may prove 


unsatisfactory. 
It might seem logical to relate h” directly to” by an equation of the form 
ON bi a (28) 
similar to eqn. (166). This leads to the equation 
a ={1/aF(A"/a)} sinh{(a/2A") (1 —4Bpd"/Ea)!?}/sinh (a/20") se (29) 


4 T=2 
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which reduces, when H->0, to «=(1/A")exp(—Bp/E), so that, by comparison 


with (18) and (20), we must put A” =2’. 
Some calculations have been made of sparking potentials based on use of 


eqn. (29), with X” =’, and including the correction for electron recapture by the _ 


cathode. Compared with Penning’s observed values, the calculated minimum 


sparking potentials, and the magnetic fields at which they occur, are too small by an | 


order of magnitude, and calculations based on eqn. (18) are, in fact, in much better 
agreement with observation than calculations based on eqn. (29). It does not 
seem, therefore, that eqn. (29) is a satisfactory expression fora. _ 

We may approach the question of the value to be given to h” in another way. 
To bring eqn. (18) (for the non-magnetic case) into line with observation it is 
necessary to let the mean free distance X’ travelled in the direction of the electric 
field differ by a certain factor from A, the mean free path defined in other ways; 
X'/A=L'|/L =1/AL, of the order of two or three. When a magnetic field is present, 
we may make a similar correction, by taking h” to be equal to 1/AL times h, the 
mean free distance in the direction ne the electric field, calculated fromA. Thusif, 
in (27), we put h” =A(A)/AL =aF(A/a)/AL and X" =A, we get 


a ={AL|aF(A/a)} sinh{(a/2A) (1 —4BL/Ea)!}/sinh(a/2X) —... (30) 


A comparison of the values of « given by eqns. (18) and (30) shows that, if 
4BL/Ea<0-9, they never differ by more than about 50%, so that under these 
circumstances, using (30) instead of (18) will not make any major difference to the 
sparking potential curves. In fig. 4 the line corresponding to BL/Ea=1 is 
shown dotted. (From (3), 4BL/Ea=BLeH?/2E?m x 108 =0-84(H/E)?, so that 
this line also corresponds to H/E=1-1.) Since, according to eqn. (30), « is zero 
at all points on and below this line, the sparking potentials calculated using eqn. (30) 
must lie above it. 

It seems then that the use of eqn. (30) will result in sparking potential curves in 
general qualitative agreement with observation, but not greatly different from the 
curves of fig.4. It must be remembered that many arbitrary assumptions underlie 
eqn. (30), and it is impossible to say to what extent it is a valid approximation to « 
until measurements of « in a magnetic field are available with which it may be 
compared. 

Whatever its range of validity may be, it will certainly not be valid when H>E. 
In this region our basic assumptions are no longer reasonable approximations to 
the truth, for clearly, when electrons cannot attain from rest in the electric field 
sufficient energy to ionize, our assumption of completely inelastic collisions is 
seriously astray. A theory satisfactory under these conditions must take account 
of elastic collisions, and the possibility of an electron gaining sufficient energy to 
ionize as the result of two or more favourable collisions. This would involve an 
analysis of the situation in which the energy of random motion of the electrons is 
comparable with the energy acquired from the field between collisions. 


This region, in which H/F is relatively high, is just the region where the 


calculated potential curves fail to show even qualitative agreement with observa- 
tion, for Penning finds that at high magnetic fields the sparking potential 
decreases again, while all our calculations give a steady rise. 

It seems that a better understanding of the nature of the observed sparking 
potential curves must wait on experimental measurement of « as a function of 
E, Hand p. If such measurements were available it might be possible to find an 
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empirical analytical expression for « to replace T'ownsend’s relation or, failing that, 
to work graphically from the experimental curves. 
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APPENDIX 


To show that x-1 coth x — x~* is a monotonic function of x for x>0. 

Let y=! cothx—x, 
then dy/dx = —x~* coth x — x1 cosech?x + x3 

= (2 sinh?x— x cosh x sinh x — x?)/(«? sinh 2x) ; 
hence, when x>0, dy/dx has the same sign as 
cosh 2x—1— $x sinh 2x — x? = 5 eae 
n=1 (2n)! 
MED he wang he = 1 A 
- 2am? Ze 20" en - Ge 
But 4(2n—1)!<(2n)! for n>2, so that dy/dx<0 when x>0 and y is mono- 

tonic if x>0. 
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ABSTRACT. Experiments are described which investigate the metal transfer which 
takes place when a radioactive hemispherically ended slider is pressed normally into a flat 
surface of the same or another metal. 

It is found that metal transfer in the form of a number of small fragments is always 
observed, and this supports the view that metallic junctions are formed between metal 
surfaces in contact. Experiments in which the load is varied suggest that at the higher loads 
the oxide layer on the metal surfaces is broken up to a larger extent than at lower loads, anda 
more than proportional increase in. metallic interaction and transfer takes place. An 
analogous effect is observed in the presence of boundary lubricants. 

When the surfaces are impacted together, very similar results are obtained. Somewhat 
less pick-up is observed than for static loading, and the difference is probably due to the fact 
that it takes time for strong junctions to be formed. Impact experiments with surfaces 
covered by lubricants show that a lubricant layer may be trapped between the surfaces, and 
this produces a large reduction in pick-up without greatly reducing the amount of plastic — 
deformation. : 


§1. INTRODUCTION 
HE theory of Bowden, Moore and Tabor (1943) on the nature of the 
friction force postulates that, when clean metal surfaces are pressed 
together, metallic junctions are formed between them and these have to 
be sheared before sliding can occur. 

A direct test of this theory would consist of a measurement of the adhesion 
of clean metal surfaces pressed normally together. However, in general no 
such adhesion is observed, because when the force pressing the surfaces together 
is removed, a release of elastic stresses occurs and the junctions are broken one 
by one. In the case of soft metals such as indium these elastic stresses are small, 
and McFarlane and Tabor (1950) and Parker and Hatch (1950) have recently 
shown the considerable adhesion of clean indium surfaces. 

Another method of investigating these metallic junctions, applicable to hard 
as well as to soft metals, consists in making one of the metal surfaces radioactive 
and then, after pressing the surfaces together, measuring the radioactivity of 
the other. The transfer of radioactivity under these conditions gives direct 
proof of the formation of metallic junctions. 

A point of some interest is the time factor when the surfaces are pressed 
together. Crook and Hirst (1950), working with indium, have shown that the 
force of adhesion is much smaller for impacted surfaces than for surfaces pressed _ 
slowly together. 

Experiments were therefore carried out in which a metal slider with a 
hemispherical end, previously made radioactive, was pressed into a flat surface 
of the same or another metal, slowly or under impacting conditions, and the 
transfer of activity to the other surface investigated. 


* Now at the Lubrication Laboratory, Department of Mechanical Engineering, Massachusetts 
Institute of Technology. 
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_ §2. EXPERIMENTAL 
_ In preliminary experiments the hemispherical slider was fixed into a simple 
holder and pressed into the other surface by a spring. ‘The apparatus was not 
_very rigid and no precautions were taken to eliminate tangential motion as the 
two surfaces came into contact. Under these conditions it was found that some 
lateral motion—estimated at 10-4 to 10-3 cm—occurred, and this had a profound 
effect on the amount of radioactive transfer, so that later experiments were 
carried out with a much more rigid apparatus. A special holder was made which 
fitted a Vickers Hardness Tester, and the hemispherically ended slider then 
occupied the position normally taken up by the diamond indenter. Loads of 
1-100 kg could be applied, and the time of loading was normally 10 seconds. 
For the impact experiments the slider was screwed into the bob of a compound 
pendulum (fig. 1). Impact energies of 104 to 108 ergs could be used. The time 
of impact was estimated (cf. Bowden and Tabor 1941) as about 10~ sec. 


LATHE 


TOOL POST PENDULUM 


ELECTROMAGNET 


Bee 
RADIOACTIVE 
SLIDER 


FLAT SURFACE 


Fig. 1. Sketch of apparatus used for the impact experiments. 


Before the experiments the surfaces were cleaned to remove surface 
contamination, the radioactive hemispherical slider by rotation in the shaft 
of an electric motor whilst pressing fine carborundum paper against its surface 
and the flat surface by lapping on fine carborundum paper under water, followed 
by drying on clean filter paper. 

The radioactive sliders were prepared by irradiating sliders of the required 
shape in the B.E.P.O. at Harwell. ‘The amount and nature of the resultant 
activity varies widely for different metals. Sliders made of copper, cadmium, 
platinum and stainless steel (18%, chromium) have proved suitable. 

The radioactivity transferred on to the flat surface after indentation 
experiments had been carried out was measured by placing the flat surface on a 
sheet of photographic film, Ilfex x-ray film being normally used. After 
development of the film, a blackened region corresponding to the transferred 
radioactivity could be seen. The method was made quantitative by visual 
- comparison of the amount of blackening with that produced when a thin foil of 

the same metal, irradiated together with the slider and thus of the same specific 
activity, was placed on a piece of the same film for known periods of time. The 
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amount of blackening was found to depend only on the total amount of radiation 
received by the film, not on the time or intensity separately, i.e. the reciprocity _ 
law was obeyed (Rabinowicz 1950). a 

The sensitivity of the method is such that, in the case of copper, metal | 
transfer of 10~! g per impression may be fairly readily detected. The accuracy _ 
of the method depends on whether it is required to determine the amount of 
pick-up absolutely or to compare two impressions made at the same time. 
In the former case the probable error is about 30%, in the latter case 15%. — 
The reproducibility of normal loading experiments is fairly good, results 
obtained in different experiments, nominally the same, showing mean variations 
of about 40%. 


Si [oe 
0) 10x. () 100u 
(a) 


Fig. 2. Autoradiographs of copper surface after indentation by radioactive copper hemisphere 
under a load of 10 kg. 


(a) Clean. Metal transfer=6 x 10-" g. 
(6) Surface covered with cetane. Metal transfer=10~!° g. 


(Exposures have been adjusted to give equal intensities of blackening.) 


§3. RESULTS OF STATIC LOADING EXPERIMENTS 


(i) Effect of Lateral Movement 


Experiments made with the first, not very rigid, apparatus showed the great 
importance of eliminating lateral movement when bringing the two surfaces 
together and removing them. In a typical case, copper pressed into copper, 
clean, under a load of 4 kg, the amount of transfer varied between 3 x 10-8 and — 
3x10 g, the exact value depending probably on the amount of lateral 
movement whilst making the impression. When great care was taken to 
eliminate this motion consistent values of 10-!° g were obtained. 

This value was also obtained when experiments were carried out under the 
same conditions in the Vickers Hardness Tester, in which lateral motion is 


extremely small, and thus may be considered with some confidence as that for 
pure normal loading. 
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(ii) Nature of the Metal Transfer 


_ Figure 2 (a) shows an autoradiograph produced by clean copper pressed into 
copper using a load of 10 kg. It is seen that metal transfer has taken place as a 
number of fragments, each of which has produced a blackened zone. 

The total amount of pick-up for this impression was 6 x 10-!g. Assuming 
that most of this pick-up is contained in about 20 fragments equal in size, we 
can deduce for each (if spherical) a diameter of 2x 10-4cm, or 250004. This 
shows that it is probably metallic copper and not just the oxide layer that has been 
transferred, since the latter for copper is at most a few hundred Angstréms 
thick (Evans 1946), and a flat flake of oxide of the required mass would have 
an area bigger than the size of the metallic junctions. At 4kg the fragment 
_ diameter is 100004, and the same considerations apply. 


(ui) Use of Different Metals and its Effects 


Impressions were made for a number of metal combinations, and the results 
are collected in table 1. The pick-up was found to be independent of the time 
of loading within the range 1-50 seconds. 


Table 1. Metal Transfer produced using Clean Surfaces of the Metals stated. 


Load =4kg 
Metal transfer Metal transfer 
Metals (10-1g) Metals (10-°g) 

Cadmium on cadmium 2 Cadmium on steel 3 
Copper on copper 1 Copper on steel 0:2 
Platinum on platinum 4 Steel on copper 1 
Steel on steel ( Copper on cadmium 10 

Platinum on silver 15 


It will be seen that the softer metals give higher pick-up than the harder ones, 
although this effect is fairly small. Unlike metal combinations give as much 
pick-up as do like ones. 


(iv) Repeated Impressions and their Effect on the Pick-up 


The final finish given to the top surface was produced by rotating it against 
600-grade carborundum paper, and thus a surface with small but sharp 
asperities was produced. It was to be expected that the amount of pick-up 
would be large until these asperities had been worn off, and this behaviour was 
found. The first impressions produced more metallic transfer than subsequent 
ones; after the fifth impression the amount of pick-up was constant, as is shown 
in table 2. Each impression was made on a different part of the lower surface. 


Table 2. Clean Copper on Copper. Load =20kg 
Impression 1st 2nd 3rd 4th 5th 6th 7th 8th 
Transfer (10-1°g) 600 250 150 100 100 50 50 50 


In carrying out experiments it is important that sufficient impressions are 
taken to eliminate this initial effect. All other values in this paper correspond 
to the steady-state value (50 x 10-1 g in table 2). 


634 R. Rabinowicz 


(v) Effect of Surface Roughness 


Experiments were carried out using finishes of different grades of carborundum | 
paper for the lower surface. The results are shown in table 3 (the roughness __ 
was measured by profilometer). 


Table 3. Clean Steel on Copper. Load=20kg 


Finish (grade) 600 220 80 
Roughness (microinch) 1 2 10 
Pick-up (1071g) 50 60 70 


It will be seen that the changes in roughness have had very little effect on — 
the amount of pick-up in the steady state; the differences actually found are so . 
small as to be only barely significant. | 


(vi) Variations in the Load and its Effects 


Experiments were carried out in which the load was varied from 1 to 100kg. 
In all cases an arrangement of a hard indenter and a softer lower surface was 
used, since otherwise the geometry of the indenter would have changed during 
the course of the experiment. ‘ 


30 
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Fig. 3. Metal transfer as a function of load for a clean copper indenter and an annealed 
copper surface. 


Figure 3, for hard copper on annealed copper, clean, shows that at low loads — 
the pick-up is proportional to the load, but at high loads deviations occur. 
Similar data were obtained with the other metal combinations, and the effect 
appears to be a general one. All the experimental results obtained could be 


fitted to curves of the type Pick-up=c,L+c.L2, where c, and c, are constants . 
and L is the load. 
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(vii) Use of Lubricants 


It was found that use of a lubricant reduced the amount of metal transfer, as 
may be seen in table 4. 


Table 4. Copper on Copper : Amount of metal transfer in 10° g 


Load (kg) 4 10 50 
No lubricant 1 6 150 
Cetane : 0-06 1 50 
Solution of stearic acid in cetane 0-01 0-25 15 


The solution of stearic acid, known to be a better boundary lubricant than 
cetane itself in the case of sliding surfaces, produces less pick-up also in the 
case of stationary surfaces. Both lubricants lose much of their effectiveness at 
the higher loads. 

Figure 2(5) shows that the use of a lubricant has not greatly reduced the 
number of metal fragments transferred. It is the average size of the fragments 
that has been decreased. 


§4. EXPERIMENTS WITH IMPACTED SURFACES 

Experiments were carried out to check whether the behaviour of the impacted 
surfaces was the same as that of the statically loaded ones as described in the 
preceding sections. Very similar results were obtained. For example: for a 
number of impressions at a constant impact energy the metal transfer decreased 
for the first four or five impressions and then became constant, in agreement 
with the results shown in table 2; when the impact energy was varied it was 
found that at low energies the pick-up was proportional to the energy of impact, 
but at high energies deviations occurred and the pick-up was greater than would 
be given by a proportionality relationship (cf. fig. 3). 

The important experiment, however, was a quantitative comparison of the 
amount of metal transfer with that obtained under static loading. 


(1) Comparison of the Metallic Transfer for Statically and Kinetically Loaded 
Surfaces ig 


To test how the amounts of pick-up compared in the two cases, loads and 
impact energies were selected such that indentations of equal size were produced. 
‘Table 5 shows the results obtained. 


Table 5. Hardened Copper on Annealed Copper. Clean 


PDiameter of impression (mm) 0-4 0-8 f2 
Mass of pick-up { impacted surfaces 3 20 100 
(10-%g) —’\_static loading 5 30 100 


The amount of metal transferred is approximately the same in the two cases, 
Static loading produces somewhat more pick-up, but the difference is small. 


(ii) Impact Experiments with Lubricated Surfaces 


Experiments were carried out with surfaces covered by various paraffin oils 
and a range of silicone oils. It was found that metal transfer still took place, 
but on a reduced scale compared with that for clean surfaces. ‘The reduction 
was in two parts; firstly, there was an overall decrease in pick-up, similar to that 
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observed during static loading in the presence of a lubricant ; secondly, there was 
another effect as a result of which no metal transfer took place at the centre of the 
indentation. Figure 4 shows the results of an experiment in which radioactive 
annealed copper hemispheres impacted hard copper surfaces, both clean and 
covered with a viscous silicone oil. The profile of the indenter as measured by a 
Talysurf profilometer and an autoradiograph of the flat surface after impact are 
shown in each case. As a result of the use of the lubricant, metallic contact at 
the centre of the indentation has been prevented and additional deformation 
produced. 

Further experiments were carried out in which hardened steel balls of 2-5 cm 
diameter were impacted into copper surfaces. In this case the metallic 


feore 
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N 


Slider of annealed copper after impacting The same, but surface covered with oil of 
a hard copper surface. Clean. viscosity 1000 cst. 


| eee 
0 500% 0 300m 
Pick-up on to the hard copper. The same, lubricated. 
Fig. 4. 


interaction was studied by measurement of the electric conductance between the 
steel and the copper during impact. With clean surfaces the electrical 
measurements showed that, as expected, metallic contact occurred at impact. 
When an oil of viscosity 10 centistokes was spread on the surface, the collision 
process showed the same characteristics. Deformation of the surface took 
place and metallic contact occured. When the surface was covered with an oil 
of viscosity 500 centistokes a marked change was observed. No metallic contact 
or break-through occurred during the impact, but the copper surface still showed 
very marked deformation. ‘I'ypical indentations are shown in fig. 5. It will be 
seen that the deformation produced in the presence of the oil is somewhat 
different in shape from that obtained with clean surfaces. It is conical rather 
than spherical, with a marked ‘ pip’ or dimple at the centre. 
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§5. DISCUSSION 


These experiments on the pressing of a hemispherical slider into a flat metal 
surface show that, under very varied conditions, metal fragments are transferred 
from one surface to the other, and this shows in a direct way that strong ee 
are formed between metal surfaces in contact. 

The amount of metal transferred is small. In a typical case, copper on steel, 
clean, load 4kg, the pick-up is 2 x 10-" g, which corresponds to a uniform layer 
of thickness 0-1 A over the area of the indentation. This value may be compared 
with that for sliding surfaces. For copper on steel, clean, load 4kg, Rabinowicz 
and Tabor (1951) have found metal transfer equivalent to a uniform layer of 
thickness 500 A over the sliding track, i.e. there is a difference by a factor of 5 000 
between the two cases. It follows that even small amounts of lateral motion 
while the static indentations are being produced may lead to comparatively 


(a) (b) 
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Fig. 5. Micrographs and profiles of copper surface after impacting by a steel hemisphere. 
(a) and (c), clean; (b) and (d), surface covered with a silicone oil of viscosity 500 cst. 


large amounts of metallic transfer, far greater than for purely static loading. For 
example, a lateral movement of only 10-?cm may lead to a hundred-fold 
increase of pick-up over that for pure static loading. 

We must now consider the fundamental differences between the processes 
of sliding and of static loading. When the metal surfaces are brought into 
contact, junctions (usually referred to as ‘metallic junctions’) are formed 
between the prominent asperities (Bowden, Moore and ‘Tabor 1943). Each 
junction will consist partly of a region where oxide layers still cover each metal 
and constitute the interface, partly of a region where the oxide layer has been 
broken and contact is established between the bare metals. Only in this latter 
region is there any likelihood that the strength of the junction is as great as that 
of the bulk metals and, since this is a prerequisite for the transfer of metal 
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fragments, only in this region can we expect the occurrence of pick-up. It 


follows that the fragments formed during static loading will have dimensions — 


smaller than that of the junction (fig. 6). 


~~ 


If, however, the metallic junction is broken by shear, as during sliding, ; 


conditions are different. In this case, owing to the different geometry, the 
breaking of the junction is preceded by a shearing action which will tend to 
break up the oxide layer and produce metallic contact over the whole of the 
junction. Then, if a fragment is produced, its dimensions will be those of the 
junction. In addition, the whole junction may increase in size before shearing 
takes place (McFarlane and Tabor 1950). 


(q) 


Fig. 6. Hypothetical appearance of a typical metallic junction showing the oxide surface layer. 


(a) As formed during static loading. Note partial break-up of oxide layer. 

(b) Surfaces separated normally. Note formation of small transferred fragment. 

(c) Surfaces separated tangentially. Note shearing of oxide layer, resulting in a large 
transferred fragment. . 


Considering next the variation of pick-up with load, it is found that the 
pick-up increases more than proportionately at the higher loads. This is probably 
due to the fact that the more severe deformation of the surfaces at the higher 
loads leads to a more widespread breakdown of the oxide layer, and thus. 
metal-to-metal contact is made over a larger fraction of the true area of contact, 
which in turn produces a greater amount of metal transfer. An analogous effect 
has been observed in the case of sliding surfaces by Whitehead (1950), who 
observed the increase of coefficient of friction with load within a certain range of 


loads and was able to show that this did correspond to different degrees of — 


breakdown of the oxide layer. 
The behaviour in the presence of a lubricant is consistent with our present 
views on the nature of lubrication. Boundary lubricants act by reducing the 


area over which metallic contact occurs. ‘Thus in the presence of the lubricant — 


some metallic transfer occurs, though less than in the absence of the lubricant. 
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The presence of the lubricant diminishes the size, rather than the number, of 
metallic junctions. In these experiments we find that, at higher loads, the 
lubricants are less effective in reducing the amount of pick-up, and this may be 
due to the fact that, as in the presence of the oxide layer considered previously, 
the lubricant layer is broken up at excessive deformation of the surface (McFarlane 
and ‘Tabor 1950), and as a result more metal-to-metal contact occurs. Again 
the results are analogous to those for sliding surfaces, where it is found that in 
the presence of good boundary lubricants the amount of metallic transfer 
increases more than proportionately with the load (Rabinowicz and Tabor 1951). 

From the fact that impacted surfaces—in which the time of contact is of the 
order of 10-*sec—give nearly the same value of pick-up as do surfaces in contact 
for times in the range 1 to 50 seconds, it appears that the formation of metallic 
junctions takes place very rapidly. However, the results show that the metallic 
transfer is somewhat less with impacted surfaces—by about 25%. This fits 
in well with the results of Sampson, Morgan, Reed and Muskat (1943), who 
_ found that surfaces in stationary contact for very short times gave lower values 
of the static coefficient of friction than surfaces that had been in contact for 

some seconds. 

At first sight the results are in disagreement with those of Crook and Hirst 
(1950) who, working with indium and measuring the force of adhesion, found 
that impacted surfaces gave less metallic interaction than statically loaded ones 
by an order of magnitude. But it must be borne in mind that the greater strength 
of junctions formed over a long time of contact is probably due to creep. Indium 
shows marked creep at room temperature. Hence plastic flow will continue to 
take place under the applied load, and this will continue to strengthen the metallic 
junctions. The metal used in our tests, copper, shows far less creep at room 
temperature, consequently the continued application of the load will not 
strengthen the metallic junctions to anything like the same extent. 

The final point of interest in these experiments is the great difference in 
behaviour of statically and dynamically loaded surfaces in the presence of a 
lubricant. One result of using the lubricant is an overall reduction in the pick-up, 
and this effect is present in both cases. With the impacted surfaces there is, in 
_ addition, a cushioning effect produced by the lubricant. Owing to the short time 
of the collision, the lubricant at the centre of the indentation cannot escape and is 
trapped, and this in turn leads to the absence of pick-up at the centre of the 
indentation. Under suitable conditions metallic contact may be completely 
eliminated. However, deformation of the metals still takes place, and the 
indentation is characterized by a ‘pip’ or dimple at the centre corresponding to 
the trapped lubricant. 

These results have important technological significance. ‘They show that 
in cases where the loading is impactive (e.g. in gear teeth), the presence of the 
lubricant, while reducing, or even eliminating, metallic contact, may have little 
effect on the amount of deformation produced in the metal surfaces. 
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ABSTRACT. ‘This paper describes an experimental investigation designed to evaluate 
the effect of width and length on rates of evaporation from small, saturated, plane, rectangular — 
areas placed in a turbulent boundary layer. An approach is made to the problem of esti- 
mating the distribution of temperature on the evaporating surface by using a number of 
thermojunctions situated at various points under the surface. The higher rate of evapora- 
tion near the edges of the area is well illustrated by this method, greater temperature 
differences between the surface and air stream being observed near the edges due to greater 
latent heat effects. Provided the crosswind width of the evaporating area is not less than 
about half the downwind length, Sutton’s two-dimensional solution of this limited area 
problem is found to lead to tolerable accuracy in the circumstances tested. 


§1. INTRODUCTION 

N a previous paper (Davies and Walters 1951) the authors have described 
| experiments, carried out in an open circuit wind tunnel, designed to test 

the dependence of rate of evaporation on width of area. These experiments 
have now been extended to areas of length up to 15 inches and widths varying 
from 1 to 6 inches. ‘The technique followed was again the ‘gravimetric’ 
method, which has been previously described fully (Pasquill 1943). In order 
to make as close a comparison as possible between theory and experiment for 
absolute values of rates of evaporation a good estimate is needed of the saturation 
vapour concentration at the surface. This is accomplished by soldering 
thermojunctions at various points on the lower side of a thin brass plate, on the 
upper face of which the saturated pieces of filter paper are placed. The 
temperature distribution and hence the surface vapour concentration distribution — 
is then determined by connecting the thermojunctions to a calibrated 
potentiometer. ‘Ihe temperature differences between the air stream in the 
tunnel and the evaporating surface, due to the latent heat of vaporization of the 
liquid, are found to be appreciably higher near the edges of the area than — 
elsewhere inside the area. ‘This gives further confirmation of the authors’ 
previous suggestion (Davies and Walters 1951) that rate of evaporation is 
greater near the edges of the area. 
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The greater range of length and width now tested enables an interesting 
comparison to be made with experiment of Sutton’s theoretical results on 
évaporation (Sutton 1949) and the solution obtained by the present authors 
{Davies and Walters 1951). The former treatment neglects the enhanced rate 
of evaporation near the longitudinal edges and the latter assumes the rate of 
evaporation to be independent of the downwind distance from the leading edge. 


§2. EXPERIMENTAL ARRANGEMENTS 

The general method adopted for measuring rates of evaporation followed 
that already described in previous work (Davies and Walters 1951), except that 
the new glass plate used was about 17 inches long and 8 inches wide, and a 
closer approach to the temperature of the evaporating surface was attempted. 
The width of the plate was limited by the dimensions of the arm of the available 
balance. The length of the plate was determined by the length of the flat 
portion of the faired holder, into which a recess of the necessary dimensions 
was made. 

When the required evaporation observations had been completed the glass 
plate was replaced by a thin brass plate of the same area and 0-04in. thickness. 
‘Thermojunctions were soldered to the underside of the plate at the points 
marked A, B, C, D, E, F in fig. 1, and temperatures measured by means of a 
calibrated potentiometer. 


Mean Wind Direction 
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Fig. 1. Diagram indicating posi- 
tions of thermojunctions at 
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A rectangular piece of filter paper, 13$in. long and 5} in. wide, was placed 
on the brass plate in the holder in the position shown in fig. 1, and saturated 
with aniline, with the holder resting centrally in position on the floor of the 
working section of the wind tunnel. The air speed in the tunnel was adjusted 
and maintained at the constant value, 1330 feet per minute, used during the 
evaporation tests. Readings were made of the thermojunction temperatures 
(at points A, B, C, D, E, F), the air temperature in the tunnel stream, and the 
air temperature in the neighbourhood of the potentiometer, which was sheltered 
as much as possible and set just outside the tunnel door. After five or six minutes 
the readings settled down to steady values until the filter paper began to dry, 
which took about forty minutes. ‘These temperature results together with 
the new evaporation results are given in the following section. 
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§3. EXPERIMENTAL RESULTS AND DISCUSSION 

The temperature measurements are discussed first, as they lead to a correction 
factor which must be evaluated in order to compare theory and experiment for 
absolute magnitudes of rates of evaporation. During the experiment the 
temperature of the air stream in the tunnel remained constant and was found to 
be 0-5°c cooler than room temperature near the potentiometer. After this 
difference was allowed for, the temperature differences between the tunnel air 
stream and the evaporating surface were evaluated. ‘These were 1-7°c for 
positions A, B, C, E, 1-4°c for D and 1-8°c for F. 

The temperature decreases near the edges (A, B, C, E, F) of the area are 
observed to be about 20% greater than in the centre (D) of the area: this 
indicates a considerably greater rate of evaporation near the edges than at 
points in the interior of the area, and confirms the result of a previous paper 


(Davies and Walters 1951). It is interesting to note that, in the special — 


circumstances of this previous work, the theoretical variation of rate of 
evaporation across an evaporating area is proportional to y~™?(w—vy)-™?, 
where w is the width of the area, y the distance of the point considered from 
an edge, and m is a turbulence parameter, whose value was found (Davies and 
Walters 1951) to be 0-15. This formula predicts that, at a point situated jin. 
(approximately the position of thermojunctions C and E) from an edge of an 
area 54in. wide, the rate of evaporation is about 15° higher than in the centre. 
This is in reasonable agreement with the experimental value of about 20%. 
The temperature of the outer surface of the filter paper is assumed to be the 
same as that of the underside of the brass plate, as given by the thermojunction 
readings. If At is the actual difference in temperature between the two surfaces 
concerned, gq the flux of heat, k, and k, the heat conductivities of the brass and 
filter paper—aniline layers respectively, x, and x, the thicknesses of these layers, 


then 
aima( a+ ze). nae (3.1) 


From measurements of x), x, etc., the value of Az¢ is found to be 3 x 10-®°c, 
and this is a maximum effect, as we have assumed that all the latent heat of 
evaporation is supplied through the brass plate. This result justifies the 
assumption that the temperatures measured are those of the actual evaporating 
surface. 

It is assumed in theoretical treatments of evaporation problems (Sutton 
1949) that the saturated vapour pressure of the evaporating liquid is constant 
over the area and consequently the temperature of the surface should strictly 
be kept constant over the area. However, although the temperature decreases 
in the present work varied appreciably over the surface, the actual temperatures 
could be regarded as effectively constant. In the conditions considered the 
surface temperature in the centre was 17-1°c and along the edges it was 16:8°c, 
the tunnel air temperature being 18:5°c. At this temperature a suitable mean 
reduction factor thus appears to be 1:5°c. If 17-0°c is taken as a mean surface 
temperature, the variation in actual temperature on the plate is extremely small, 
and the consequent variation of saturated vapour pressure is also small (from 
the appropriate curve for aniline), and so the assumption of constant surface 
temperature leads to little error. It is reasonable to suppose that the factor of 
1:5°c will not vary significantly over the small range of room temperatures 
considered in the evaporation tests. 


a 
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Several measurements of rates of evaporation have been made which are 
additional to the results already given in a previous paper (Davies and Walters 
1951). For a constant length / of area and a constant mean air speed, log E was 
plotted against log w, E being the rate of evaporation in lb. per minute and w the 
width in inches. ‘The results are shown in fig. 2 in the form of two sets of points; 
the conditions tested are described in the figure. Twenty-two points are plotted 
in each set and lines of regression are computed and drawn. The slopes are 
found to be 0-84 and 0-86, in good agreement with the theoretical result of 0:85, 
described previously (Davies and Walters 1951). 

Measurements to determine the dependence of E on / were also made in 
the conditions detailed in fig. 3: the points corresponding to the lengths 
3-5, 5-0, 13-0, 15-0in. (marked as full circles in fig. 3) were obtained from lines 
of regression giving the dependence of E on w (see fig. 2). From fig. 3, E is 
found to be proportional to 1°76. The result given by Sutton’s well-known 
theory (Sutton 1949) which assumes unlimited width of area, is E proportional 
to [(+m(l2m), where m is a parameter determined from the velocity profile in 
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Fig. 2. Fig. 3. 


1) the turbulent boundary layer over the area. As described previously (Davies 
‘) and Walters 1951) the evaporating area tested was positioned in a fully turbulent 
boundary layer, and a value for m of 0-15 was obtained. ‘This leads to a theoretical 
value for (1+m)/(1+2m) of 0-88. The consequences of this difference in 
theoretical and measured dependence on length is clearly seen in the table. 
The considerable range of length and width of area now tested experimentally 
makes a comparison of available theoretical treatments interesting. In the 
table some of the results of the evaporation measurements are given together 
with corresponding theoretical values computed, using Sutton’s (1°49) formula 


To ae ee) I eet 7 eta me (372) 

and the formula obtained by the present authors (Davies and Walters 1951) 
1—nfm/2 = 1—m 

Re 27x a,U,-"f™?T (1 — m/2)lw Ree horse: (3.3) 


D\(m/2)1(3— m/2)0($ —m/2) 
where K has been tabulated (Pasquill 1943), y is the saturation vapour 
concentration at the temperature of the surface, a, and m have been specified in 
previous work and are easily calculated from the observed value of m, and U, is 
the mean air speed at a convenient reference height. ‘The f/ factor in eqn. (3.3) 


has been discussed previously (Davies and Walters 1951) and enters only to 
2,U-2 
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small power, m/2. Near the surface in the earth’s atmosphere it has been 
suggested (Davies and Walters 1951) that an appropriate value, in adiabatic 
conditions, is 1-6. Some recent results (Townsend 1951) on the variation of 
eddy velocities with distance from the wall in a turbulent boundary layer appear 
to some extent to be relevant here. Assuming that the results obtained by 


Townsend apply in the wind tunnel used by the authors, a value for f of 1-0 is. 


computed at a height of about 0-4 in. over the evaporating surface and of 1-6 at a 
height of about 0-2 in., the value increasing as the surface is approached. A mean 
value for f is needed over the height of the diffusion zone which is probably about 
0-4in. thick at a distance of 4-:0in., say, downwind from the leading edge. An 
exact mean value for f clearly presents a formidable problem and in the present 
calculation a value of 1:6 is taken by virtue of the above data, the choice of f 
influencing only slightly the computed theoretical rate of evaporation from 
eqn. (3.3). 


Experimental and Theoretical Rates of Evaporation from Saturated, Plane, 
Rectangular Areas 


Air temperature =18-5° c, mean air speed = 1 330 ft/min, m=0-15, mean evaporating 
surface temperature =17-0° c 


from eqn. (3.2) from eqn. (3.3) 
() (2) (3) (4) (5) (6) (7) 
3308) 4-0 3835p 0m 3°:89><10n* 102% 12945¢ 10-5 51%, 
5:0 4-0 403 sear pA 109% Paige ee 56% 
15-0 4-0 HO? > Als Rs 129% S32 mes 76% 
5-0 2D) 363 20a, 3:34 -,, 101, 1°85 Be 56% 
15-0 DS eS O Saar. S84 Ss 120% Sos.» 3: 769% 
Dis 4-0 ilowil 6 1b ns 96% 


(1) length (in.); (2) width (in.); (3) experimental rate of evaporation (lb/min) ; 
(4)-(7) theoretical rates of evaporation: (4), (6) lb/min, (5), (7) as percentage of (3). 


* This result applies to the region between 13°0 and 15-0 in. downwind from the leading edge. 


Sutton’s formula (3.2) for rate of evaporation neglects the effect of lateral 
diffusion and consequently of the longitudinal edges, but incorporates the effect 


due to the leading edge. The formula (3.3) obtained by the present authors — 


is based on a simplified problem of diffusion from an infinitely long area, with 


finite width, so that the effect of the longitudinal edges on rate of evaporation was — 
approximated to. Both treatments make use of the concept of eddy diffusivity 


together with the notion that the transfer of mass is linked to that of momentum. 


§4. CONCLUSIONS 


The results given in the table show that provided the crosswind width of the | 


evaporating area is not less than about half the downwind length, the two- 


dimensional solution of this limited area problem due to Sutton is found to lead 


to tolerable accuracy in the circumstances tested. 


It is interesting to note that the experimental value for average rate of | 
evaporation per unit length in the range 13-15 in. from the leading edge is found | 
to exceed by only 4% the theoretical value (see table, bottom row) given by | 


eqns (33), 
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Finally, the results given in this paper provide further evidence to show that 
the approach to practical problems of turbulent diffusion and evaporation, 


initiated by Sutton (1949) using the concept of eddy diffusivity, is a powerful 


one, but in order to account theoretically for the effect on rate of evaporation 
of variations in length and width of area a full analytic, computable solution is. 
needed to cover the effects of leading, trailing and longitudinal edges. 
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ABSTRACT. A theorem, relating the torsional modulus for small torsions of a rod of 
incompressible highly elastic material, held in a state of simple extension, to the extension 
ratio and force necessary to maintain the extension, is verified experimentally for rods of 
vulcanized rubber of rectangular cross section. 


$l. INTRODUCTION 

T has been shown (Rivlin 1948) that for an incompressible highly elastic 
] material which is isotropic in its undeformed state, the stored energy 

function W must be related to the deformation through two strain invariants 
I, and I, which can be expressed in terms of the deformation. Ifa rod of sucha 
material, of circular cross section, is subjected to a simple extension of extension 
ratio A and then twisted, for sufficiently small torsions the torsional modulus 
is related to the extension ratio A and the force N required to maintain the 
extension, by the relation (Rivlin 1949) 


ucla iit son re (1) 


where a is the radius of the rod in its undeformed state. M is the torsional couple 
and % the amount of torsion, measured per unit length of the rod in its extended 
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state, so that [M/}],- _» is the torsional modulus of the rod in the limiting case 


when the torsion is infinitesimally small. It is seen that this relation is 


independent of the manner in which W depends on J, and I, for the material. 


This theorem has been verified experimentally (Rivlin and Saunders 1951) | 
using a rod of vulcanized natural rubber. Green and Shield (1951) have derived — 
a more general theorem applicable to a uniform rod of any cross sectional form — 


twisted about any axis parallel to its length. They obtain the result (see Appendix) 


xLT,- = g)=g [eo Fe hee al (2) 


where A, is the cross sectional area of the rod in its undeformed state. I9/Ap is 


the square of the radius of gyration of the unstrained cross section about the | 


position of the axis of torsion in the undeformed body, it being assumed that this 


axis moves with the body when it is deformed. Sy is the torsional rigidity about _ 


the same axis for infinitesimally small torsions of a rod having a cross section 


identical with that of the unextended rod, the rigidity modulus for the material | 


being taken as unity. 

In the present paper, the relation (2) is verified experimentally for the case 
of prisms of vulcanized rubber of rectangular cross section. For such prisms, 
in which the sides of the rectangular cross section have lengths a and b, I,/ Ap is 
given by 


LjAy=(@44- 2122 ee eee (3) 

and S,/A) is given (Love 1927, § 225) by 
Se Pe UN Mae 1 (2n+1)za 4 
ZnO) aa ‘i 


§2. EXPERIMENTAL 


The four rods used in these experiments were formed by a moulding process 
from a mix having the following composition in parts by weight: natural rubber 
(smoked sheet) 100, sulphur 2-0, zinc oxide 2-0, stearic acid 0:5, M.B.T.S. 1:0, 
Nonox D 0:5. Vulcanization was effected by heating at a temperature of 140°c 
for 45 minutes. Each of the rods was 7:46 cm in length and had steel end-pieces, 
by which it could be held in grips, bonded to it. Each of the rods had a rectangular 
cross section, the lengths of the sides a and b being given in table 1. 


Table 1 
Reference No. of Rod 1 2 B 4 
a (cm) 1-900 1:270 0-952 0-951 
b (cm) 0-477 0-478 0-478 0-950 
Ratio a/b 3-980 2:°654 1:992 1:001 


The test-pieces were subjected to various degrees of simple extension and in 
each state of simple extension the relation between the torsional couple M and 
amount of torsion % for small superposed torsions was measured, as well as the 
force N necessary to maintain the simple extension with s=0. These 
measurements were carried out in an apparatus similar to that used by Rivlin 
and Saunders (1951) in verifying eqn. (1) experimentally. The grips and 
end-pieces were so designed that in each case the torsion took place about an axis 
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passing through the centroids of the cross sections. Owing to the difficulty of 
measuring the torsion ys over the central regions of the rods, where end effects 
would not be expected to be significant, it was measured over the whole rod 
(i.e. from the relative rotation of the end-pieces) and the inaccuracy involved, 
for various degrees of extension of the rods, was determined in a subsidiary 
experiment. In all cases, it was found that the error so introduced in the value 
of % was less than 1°% and was therefore neglected. The extension ratio was 
obtained by measurements over the central region of the test-pieces where 
end-effects could be neglected, in the manner adopted by Rivlin and Saunders. 

The measurements of corresponding values of M and ¢ for various values of 
the extension ratio A were carried out over a range of values of ¢4-(0-0-3 rad/cm) 


Table 2 
Rod No. 1 Iy/Ayp=0:3198 cm?, Sp/Ayp=0-06402 cm? 
A 1257; 1-512 1-762 2-030 2-284 
N (kg) 1-906 Sou li! 3:97 4:68 Sos! 
[M/tb]y—o (kg cm?/rad) 0-600 0-730 0-744 0-760 0-765 
a 0-1965 0:2525 0-270 0-290 0-301 
B 0-191 0-246 0-273 0-289 0-298 
Rod No. 2 If Ap=0°1535 cm?, S/Ap=0-05842 cm? 
A 1-264 1-524 1-784 2-036 2-294 
N (kg) 1:34 2°16 2135 3°22 3-635 
[M/b]y—o (kg cm?/rad) 0-228 0-246 0-2485 0-246 0-246 
ot 0-108 0-125 0-134 0-137 0-142 
B 0-106 0-127 OEBi7/ 0-142 0-146 
Rod No. 3 Iy/Ay=0-09457 cm?, S,/Ay=0-05214 cm? 
A 1-264 1:518 1-769 2-032 2:283 
N (kg) 0-989 heb 7iy/ 1-960 2325 2-708 
[M/b]y)=0 (kg cm?/rad) 0-1125 0-112 0-1145 0-117 0-114 
x 0:0727 0:0758 0:0847 0-0900 0-0880 
B 0:0736 0-0825 0-0870 0-0895 0:0910 
Rod No. 4 I/Ay=0°1507 cm?, S/Ay=0-1268 cm? 
A 1-262 1-520 1-769 2-030 2:288 
N (kg) 1-95 3-13 3-965 4:68 5-295 
[Mb], =0 (kg cm?*/rad) 0-421 0-412 0-392 0:375 0-358 
ot 0-137 0-143 0-144 0-143 0-141 
B 0-139 0-144 0-1465 0-148 0-149 


such that the resulting (M, 7) curves were substantially linear. ‘The values of 
[M/],-o were obtained as the slopes of these straight lines and are given in 
table 2, together with the corresponding measured values of the force N necessary 
to maintain the extension. For each rod the values of J)/A) and S)/Ay were 
calculated from eqns. (3) and (4) respectively and the values obtained are also 
given in table 2. From the values of N, [M/i/],-) and A, given in table 2, the 
corresponding values of the quantity on the left-hand side of eqn. (2), denoted «, 
were calculated and from the values of I/Ao, S)/Ap and A, the values of the 
quantity on the right-hand side of eqn. (2), denoted 8, were calculated. The 
values of « and f so obtained are given in table 2. It is seen that they are in 
reasonably good agreement. 
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APPENDIX 
In obtaining the result (2), Green and Shield (1951) derive a corresponding 
theorem for a ce nee yen 
altl- Ty — (Zo = So)Av?/A® 5) 
y=0 “Tht eM a oe 
in which A, is the extension ratio in directions perpendicular to that of the load N. 
Equation (2) is then obtained by writing A, =1/4/A for an incompressible material. 
In eqns. (2) and (5), as derived by Green and Shield, ™ is calculated as the couple 
about the axis of torsion of the forces acting on the curved surface into which 
a plane normal to the axis of torsion is bent as a result of the torsion. ‘This is 
not strictly the torsional couple required to maintain a state of torsion in the 
stretched prism which we have measured. The latter is defined as the moment 
about the axis of torsion of the forces acting on a cross section of the prism which 
is plane and normal to the axis of torsion in the deformed state and can be 
calculated by a method similar to that adopted by Green and Shield. If the 
prism is subjected to a small torsion superposed on a simple extension, while the 
surfaces of the prism parallel to the axis of torsion are force-free, then the 
expressions for the two couples are identical. However, if we consider that in 
addition to the longitudinal force a uniform hydrostatic pressure w is acting on 
the prism, Green and Shield obtain for the couple M they have calculated 
M =p, 4{(O + VA?) [AL —AW2(Lp — Sp) ] —wUp— So)}, eee ee (6) 
while we obtain for the couple M which must be applied to the ends of the prism 
in order to maintain the state of torsion 
M=fd,4(O+ FAVE Ay Sl eee (7) 
In both eqns. (6) and (7) ® and V’ are given by 9 =20W/ol, and ‘=20W/ol,, 
where W is the stored-energy function for the material and J, and J, are two of 
the three strain invariants in terms of which it is expressed. We note that 
although M given by (6) involves w explicitly, eqn. (7) does not. Thus, if we 
interpret M in the sense in which it is used in the experiments, eqns. (2) and (5) 
are valid whether or not a hydrostatic pressure is applied to the prism. 

The apparent discrepancy between the expressions for M given in 
eqns. (6) and (7) may be further explained in the following manner. The 
presence of a hydrostatic pressure results in a uniform surface traction w acting 
normally to the curved deformed end of the prism. ‘This, when resolved in a 
plane normal to the axis of torsion, gives rise to a couple about this axis of 
magnitude —yA,*a@(J)— Sp), which is exactly balanced by the couple about the 
axis of torsion due to the action of the hydrostatic pressure w on the curved sides 
of the deformed prism between a curved end and a plane section normal to the axis. 
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LETTERS TO THE EDITOR 


Transparent Conducting Films 


Research on thin electrically conducting films, in connection with electrically heated 
windows, has led to the following observations on gold films : 

1. For very thin (e.g. 60A) gold films sputtered directly on glass, the resistivity is 
much higher than that of the metal in bulk. Also, the optical absorption is markedly 
higher than that calculated from the bulk optical constants. The transmitted colour is 
variable from blue to blue-green. On heating the film, the colour becomes first definitely 
blue and then, at a higher temperature, pink. There is evidence of grain growth during 
these changes. There are, however, theoretical grounds for expecting the transmitted 
colour of such a thin film to be pale yellow, if the metal were in the form of a uniform 
homogeneous layer. 

2. The above-mentioned properties are markedly influenced by the process used for 
cleaning the glass. Cleaning the glass with dry precipitated chalk tends to impart to the 
films a low resistance, low absorption and a colour inclining toward ‘ straw’. Chemical 
cleaning, on the other hand, tends to give the films a high resistance, high absorption, and 
a consistent blue colour. 

3. The effect of other substrates is still more marked. Oxide substrates have been 
studied in particular. On a substrate of sputtered zinc oxide the gold film is invariably 
blue, and has an exceptionally high resistance. On sputtered bismuth oxide, however, the 
film is of a pale straw colour, while the resistance and optical absorption are consistently 
low. In this case, subsequent heating for a few minutes at about 200° c results in still 
lower values of resistance and absorption, not far from those based on the bulk constants 
(with a ‘ thinness ’ correction). 

Observations 2 and 3 appear to be novel, and the practical applications of 3 have been 
patented. The following are typical results for sputtered gold films of thickness 60—70A 
(on a weight-bulk density basis). The optical results are overall figures including the 
effects of the thin glass support, and of the oxide film. The transmission factor T, and the 
absorption factor « (for light incident on the coated side of the glass support) refer to 
tungsten filament light, colour temperature about 2 700° kK, and R is the resistance between 
_ two opposite edges of a square of film. 


Substrate R a If 
Chalk-cleaned glass 60-70 24 62 
1004 bismuth oxide, on glass, with final heating 10 13 U5 


It is suggested that certain substrates, like bismuth oxide, strongly favour adhesion 
of the gold, with low mobility over the surface and a minimum of aggregation, resulting 
in good homogeneity of the films, at least up to a certain thickness. Other substrates 
favour the opposite effect. In any case, the results point to a new field for the study of 
surface forces, and of the properties of gold (and perhaps other) films in relation to their 
structure and to the bulk constants of the metal. In the present work structural and other 
observations are being extended. The sputtered films were produced by the method 
described by Preston (1950). 

The work described above has been carried out in the Light Division of the National 
Physical Laboratory on behalf of the Ministry of Supply. ‘This letter is published by 
permission of the Director of the Laboratory and with the approval of the Chief Scientist, 
Ministry of Supply. 


National Physical Laboratory, E. J. GILLHaM. 
Teddington, Middlesex. J. S. PRESTON. 
5th June 1952. 


PRESTON, J. S., 1950, Proc. Roy. Soc. A, 202, 449. 
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Area Contacts on Germanium 


The behaviour of point contacts on germanium and of p—n junctions within single | 
crystals of the material has been extensively studied by many authors, but the properties 
of metallic contacts of large area have received comparatively little attention. Some of 
these properties are described below. 'The observations refer mostly to gold contacts — 
applied by evaporation in vacuum to n-type germanium of 5 ohm-cm resistivity and a 
Hall coefficient of about 104 cm?/coulomb. They apply qualitatively also to specimens 
of higher conductivity. The germanium surfaces were etched before application of the 
electrodes. ; — 

The voltage—-current relations of contacts prepared in this way are qualitatively sumilar 
to those obtained with point contacts. 'The reverse currents usually show a partial 
saturation at small reverse voltages. The saturation currents which can be evaluated by 
a simple geometrical extrapolation are found to be temperature dependent in accordance 
with an equation of the form i,=7, exp (—¢/RT), where ¢ varies from contact to contact 
between 0:60 and 0:72 ev. No significant differences are obtained with contacts of gold, 
silver and copper and no improvement of the rectifying properties could be achieved by 
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applied to the assembly as a whole. 


forming processes of the type normally employed for point contacts. This may be due 
to the fact that equivalent current densities cannot be obtained on area contacts without 
causing burn-out. 

The slope (45°) of the lower line on fig. 1 shows that the total current for a small 
reverse voltage is directly proportional to the contact area. This holds for electrodes 
of different shapes and implies that the current density is uniform (no edge effects). The 
differential conductance in the saturation region is also found to be a linear function of 
contact area. In the forward direction the series resistance of the bulk germanium becomes 
important and the total voltage applied to the specimen as a whole is then no longer equal 
to the barrier voltage. Forward currents at a given total voltage are not directly propor- 
tional to the contact area but increase less rapidly, as shown in the figure. This is due to 
geometrical factors and to the effect of injected carriers on the resistance in series with the 
contact barrier. For small electrodes the series resistance is essentially a spreading 
resistance, and the injected carriers which reduce the resistivity of the material near the 
contact are effective in the most important region. This accounts for the diminution of 
the rectification ratio with increasing contact area, and such an effect must be expected 
in all materials in which carrier injection occurs. 
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Figure 2 shows the variation of peak back voltage and peak back power with electrode 
area. In the neighbourhood of the voltage turnover the contact temperature can no longer 
be assumed to be equal to the ambient temperature. There is, however, some evidence 
(e.g. from experiments on 2lectrodes of different shapes and different conditions of heat 
dissipation) that the turnover phenomenon is not solely due to contact heating. The peak 
back voltage and the peak back power depend on the ambient temperature in a manner 
very similar to that of point contacts on germanium and lead sulphide. 

For purposes of comparison, some results for formed and unformed tungsten point 
contacts have also been included in the figures. The forming process has little effect on 
the peak back conditions, but it sensitively affects the magnitude of reverse currents for 
small voltages. According to fig. 1 the behaviour of the (unformed) evaporated contacts 
is similar to that of the formed point contacts, a fact which must have an important bearing 
on the mechanism of the forming process and which is not yet adequately understood. 

Thanks are due to Professor R. W. Ditchburn for his interest in this work. 


Department of Physics, J. W. GRANVILLE. 
University of Reading. H. K. HEnIscuH. 
5th June 1952. 


Contact Properties of p-Type Germanium 


The interpretation of many phenomena observed on contacts between n-type 
germanium and various metals or between two specimens of n-type germanium involves 
specific assumptions concerning the existence and position of surface levels. Pfann and 
Scaff (1949, 1950) and Fan and Bray (1949) have already drawn attention to the fact that 
high-voltage rectification and transistor action can be observed on p-type as well as n-type 
germanium. However, the interpretation of these effects on p-type material demands 
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Fig. 1. Voltage-current relations for point contacts Fig. 2. Typical voltage—current relations 
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on germanium: A, n-type starting material for contacts between two specimens 
etched; B, n-type starting material after polish- of p-type germanium (oscillographic 
ing; C, p-type material produced by heat shock ~ method). 


(etched); D, p-type material as for C after polish- 
ing; E, p-type characteristic given by Torrey and 
Whitmer; F, p-type characteristic of a formed 
contact given by Pfann and Scaff. (Curves 
A-D, which refer to unformed contacts, were 
obtained by oscillographic methods. The close 
similarity between the peak back voltages of 
curves A and C is fortuitous.) 


-quite different assumptions concerning surface levels. ‘These problems still await a 
-satisfactory solution. Meanwhile, a more detailed comparison between the contact 


properties of n- and p-type germanium seemed desirable and some results are described 


below. 
Figure 1 shows typical characteristics of an unformed tungsten point contact on etched 


-p-type germanium and also on a polished surface. ‘The material was prepared from n-type 


652 Reviews of Books 


germanium by heat shock (heating to about 800° c, followed by rapid cooling). Typical 
characteristics observed on the starting material are also shown, as well as p-type curves 
given by Torrey and Whitmer (1948) and (for a formed contact) by Pfann and Scaff (1950). 
For comparison, n- and p-type rectification curves have been drawn in the same quadrants. 
The peak back voltages measured on the present specimens vary a good deal from point to 
point, more so than on the n-type starting material. For a given point, the peak back 
voltage and power decrease rapidly and continuously with increasing temperature up to 
100° c. Anomalous thermo-electric effects of the type observed by Granville and Hogarth 
(1951) on polished surfaces of n-type material were not observed on the present p-type 
specimens. As for n-type germanium, there are no significant differences between the 
average peak back voltages observed for different whisker materials. Nickel, platinum, 
tungsten and copper gave substantially similar results. 

Figure 2 gives voltage-current relations obtained for four different contacts between 
two specimens of p-type material. These relations show that a barrier of considerable 
height must exist at these junctions. 

Thanks are due to Professor R. W. Ditchburn for the support of this research. 


Department of Physics, J. W. GRANVILLE. 
University of Reading. H. K. HENIscH. 
26th May 1952. P. M. TIPPLe. 


Fan, H. Y., and Bray, R., 1949, Phys. Rev., 76, 458. 

GRANVILLE, J. W., and Hocartu, C. A., 1951, Proc. Phys. Soc. B, 64, 488. 

PFANN, W. G., and Scarr, J. H., 1949, Phys. Rev., 76, 459; 1950, Proc. Inst. Radio Engrs., N.Y., 38, 
II51. 

Torrey, H. C., and Wuirmer, C. A., 1948, Crystal Rectifiers (New York : McGraw-Hill). 


REVIEWS OF BOOKS 


Oeuvres scientifiques de Paul Langevin. Pp. 687. (Paris: Centre National de 
la Recherche Scientifique; London agents: H. K. Lewis, 1950.) In paper 
44>.; in cloth 53s. 


The appearance of most of the original papers of Paul Langevin collected together in 
one volume must be regarded as an event of great importance. This collection is not quite 
complete, as many papers which are readily accessible are not included; a list of these papers 
is given. In time this collection covers the entire first half of the twentieth century, for 
his first paper appeared in 1902 and his last in 1949 ; in substance it covers a range so wide 
as to include important papers on ionization of gases, kinetic theory, thermodynamics, 
electromagnetic theory, magnetism, relativity, radioactivity, ultrasonics, nuclear physics and 
gravitation. ‘This is indeed a contribution to physics in the grand manner. An obituary 
notice of Langevin has been published in the Proceedings of the Physical Society (1947, 59, 
1041), where tribute to his immense range of achievement was made. His earlier work, 
coming as it did at a vital period in the development of physics, played a most important 
role in effecting the transition from the classical to the modern point of view, and can be 
regarded as forming a bridge between the statistical mechanics of the late 19th century and 
the early quantum ideas of the 20th century. In particular, Langevin’s work is an excellent 
example of the power and scope of classical statistics in dealing with the fundamental. 
collisional phenomena involved in the conduction of electricity through gases. Then, 
again, it is interesting to note that he was a pioneer worker in relativity, and also that Einstein 
first learned of de Broglie’s original notions of wave mechanics from Langevin. For these 
reasons this volume of collected papers is, and probably will continue to be, of profound 
interest to the student of physics, and particularly to the young research worker, since, in 
addition to important theoretical papers, it also gives numerous papers describing his. 
experiments. 
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Langevin’s first papers dealt with the properties and motions of ions (the existence of 
free electrons in gases had not yet been demonstrated), and his interest in this subject no 
doubt originated during the time he spent at the Cavendish Laboratory under Sir J. J. 
Thomson. An outstanding problem which excited the interest of physicists at the turn of 
the century was that of the explanation of the electrical conductivity produced by x-rays or 
radioactive emanations. It was the application of classical kinetic theory to the experimental 
determinations of conductivity, and of current decay, which enabled the nature of the ion 
carriers to be found, the processes of ionization and recombination to be investigated, and 
the mechanisms of ionic transport, diffusion and recombination to be understood in terms of 
the collisions of elastic sphere atom models. 

By consideration of the polarization of neighbouring neutral atoms produced by a 
charged ion, and also the resulting attraction, effecting reduction of the free path, the 
consequent collection of molecules to form large clusters about a single ion was envisaged. 
These charged molecular groups, whether singly or in clusters, were shown to be the 
carriers of electricity. It was in this particular field that Langevin made his first 
outstanding contribution. In the series of papers in 1902-1905 (Ann. Chim. Phys., 5, 245) 
culminating in the famous paper Une formule fondamentale de théorie cinétique, he treated 
the motion of ions in gases, and rigorously derived the well-known formula for the drift 
velocity W: 

W=0-815(Ze/m)(A/C){(m+ M)/M}?, 


of an ion of charge e, free path A, r.m.s. velocity of agitation C( > W), mass m moving under a 
field Z in a gas of molecular mass M in which the ion exerts no electric forces except at 
collision (low gas pressure). The ions are considered to be nearly in statistical equilibrium 
with the gas atoms : thus restriction is now known to be unnecessary, since the formula also 
applies to the motions of electrons under a regime of elastic conditions when the electrons 
are not in statistical equilibrium with the gas atoms, provided, however, that their energy 
distribution is Maxwellian. This latter consideration, again, only affects the numerical 
coefficient, which only varies from about 2/3 to unity over a wide range of electron energy 
distribution functions. Replacement of the force Ze by that due to the concentration 
gradient leads to Langevin’s formula for the diffusion coefficient of ions. The application of 
his formulae to the measurements of the mobilities of clusters, single ions, or of electrons led 
to determinations of the mass factor and of the collisional cross sections of the ions: and, 
in fact, in the hands of Townsend and Tizard, such application suggested in 1913 the 
possibility of a variation of atomic cross section in collisions with slow electrons. That 
Langevin had peculiar insight in this subject is seen by the fact that his basic formula held 
the field for about 25 years, in fact until the discovery of exchange phenomena and the 
calculations of Mott and Massey based on quantum scattering, and was used in the inter- 
pretation of the experimental work of Tyndall and Powell. 

The other great contribution with which Langevin’s name is always associated is in the 
theory of magnetism, and this was made in 1905. The theory of electrons in the hands of 
Lorentz was making great strides, but doubt was expressed whether the magnetic theories 
of Ampére and Weber could be expressed adequately in terms of the new electron theory. 
It was Langevin who set these doubts at rest. Curie’s work on the temperature dependence 
of paramagnetism was almost completely explained by Langevin on the basis of the classical 
kinetic theory on the assumption that paramagnetic susceptibility was due to the reserve of 
permanent magnetic dipoles in each molecule. ‘The results deduced that over a large range 
of temperature paramagnetic susceptibility is inversely proportional to the absolute 
temperatures—known a3 the Curie-Langevin law. ‘This work today still forms a most 
valuable introduction to the electron theory of magnetism. The full explanation, however, 
had to await the application of the quantum theory of spin, with its restrictions on the 
orientations of the molecular magnets in a magnetic field. 

This book of 687 pages gives a curriculum vitae of Langevin and a complete bibliography 
of his scientific works, and has altogether thirteen chapters. In spite of its size, this volume 
is quite pleasant to handle, and the print is clear and easy to read. The first chapter 
describes his experimental works on ionization of gases by x-rays, his experimental and 
theoretical work on ionic mobility and recombination, including the polarization of neutral 
atoms produced by the ion charge, and the properties of secondary x-rays ; he contributed 
to these subjects throughout his life. The second chapter is devoted to considerations of 
ionization in the atmosphere and the properties of ion clusters. 
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The next three chapters are of a theoretical character and give his papers on fundamental 
kinetic theory, electromagnetic inertia, relativity, and also his classic work in the theory 
of magnetism. The sixth chapter presents a collection of his papers on relativity, and this 
is followed by chapters on radioactivity and units. Then follow some twenty pages 
devoted to a critical review of classical, quantum and statistical mechanics. Chapters ten 
and eleven are concerned with his contribution to ultrasonics, under-water acoustics and 
divers technical problems : these represent his scientific war work in World War I. He 
originated and developed the piezoelectric under-water detector, and carried out trials 
in Toulon harbour. The two final chapters are devoted to his most recent work on nuclear 
physics and gravitation, and are of a fundamental nature. 

It is not possible within the scope of a short review to discuss the chapters in detail : it 
_ is sufficient to stress the scope and quality and importance of the work described, to admire 
the clarity of its exposition, and very strongly to recommend this collection to all interested 
in the history and development of an important aspect of modern physics. 

F. LLEWELLYN JONES. 
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ABSTRACTS FOR SECTION A 


Some Post-War Developments in Magnetism, by L. F. Bares. 


ABSTRACT. A survey is given of some of the more outstanding developments in 
magnetism during the post-war years, including the great extension of our knowledge of the 
ferromagnetic domain through the theoretical work of Néel and others and by the systematic 
use of the Bitter powder figure techniques; the main ferromagnetic, paramagnetic and 
nuclear resonance phenomena and experiments on the diffraction of neutrons by anti- 
ferromagnetic crystals are also described. 


Paramagnetic Resonance and Spin—Lattice Relaxation Times in the Copper Tutton 
Salts, by D. M. S. BaGcuLey and J. H. E.-Grirrirus. 


ABSTRACT. Paramagnetic resonance absorption in single crystals of the copper 
Tutton salts has been investigated at 3:10 cm wavelength between 90°k and room 
temperature. In salts diluted with zinc Tutton salt the hyperfine lines are clearly resolved 
at 90° kK but broaden and merge into one another as the temperature is raised. From 
these measurements the spin-lattice relaxation time 7 can be estimated and is found to be 
about 2x10~*sec. It is found that 1/7 varies linearly with temperature between 
150° Kk and 290° kK. An investigation of the crystalline field in the Tutton salts shows 
that for the concentrated salts the field has a strong rhombic component above 150° k, 
whilst in the salts diluted with zinc to less than 1% copper the crystalline field has 
tetragonal symmetry at 290° k, but the tetragonal axis lies in a plane perpendicular to that 
containing the tetragonal axis at 90° k. 


A Note on the Theory of Semiconductors, by P. 'T. LANDSBERG. 


ABSTRACT. One obtains two different theories of semiconductors, depending on whether 
the free energy due to the spin of electrons in donor (or acceptor) centres is taken into account 
(approach A) or not (approach B). It is shown (i) that the former approach can be carried 
through for a degenerate semiconductor, (ii) that it applies formally to the case when the 
interaction between impurity centres and conduction electrons is taken into account, 
provided one of the symbols is reinterpreted, and (iii) that it is to be preferred. Approach 
B has, however, been used almost universally. It is shown that it corresponds to the 


unlikely situation in which fifty per cent of the electrons from each impurity level come 


from paired states, and the remaining fifty per cent from an unpaired state having the same 
activation energy. 


Quantization of Einstein’s Gravitational Field: II—General Treatment, by 
Suraj N. Gupta. 


ABSTRACT. The quantization of the complete gravitational field is carried out by 
extending the work of an earlier paper. The main obstacles in the quantization of Einstein’s 
field are overcome by expressing the field quantities in the Riemannian space as expansions 
in the flat space, and then splitting the gravitational field into the linear and the non-linear 
parts. The linear part of the gravitational field is regarded as the free gravitational field, 
while the non-linearity is treated as a direct interaction between the gravitons. ‘This 
treatment is quite general, but it suffers from the usual limitations of the perturbation 
method. 

The gravitational self-energies of the photon and the electron are also investigated by 
using the linear approximation of the gravitational field. It is found that the photon 
self-energy vanishes unambiguously, while the electron self-energy is quadratically 
divergent. 
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A New Theory of Liquid Helium, by H. N. V. 'TEMPERLEY. 


ABSTRACT. A model of liquid helium is described, which incorporates the main 
features of the Landau (1941) and London (1939) theories, making use also of Mayer’s 
concept of clusters. The ordinary process of liquefaction is associated with the formation 
of clusters in coordinate space, while the A-transition is probably associated with the — 
formation of small clusters in momentum space. The latter process is probably prevented; 
or greatly modified, in helium 3. 

Possible reasons for the appearance of superfluid effects are given, and it is suggested 
that the unique properties of helium II may be traced to the fact that the interactions 
between two atoms is of about the strength needed to form a molecule. 


The Absorption Spectra of Gallium and Indium Monofluorides, by D. WeLt1 and 
R. F. Barrow. 

ABSTRACT. Gaseous GaF and InF are formed when the metals are heated with 
aluminium trifluoride, and sufficient pressures of the monofluorides for observation of their 
absorption spectra are obtained at about 1 000° c for GaF and 1400° c for InF. The three 
ultra-violet band-systems, "II, *II,, ?11)*<*X*+ expected by analogy with GaCl and InCl have 
been photographed on spectrographs of moderate dispersion, and vibrational analyses have 
been carried out. The 411 states converge rapidly to give rather sharp dissociation limits: 
the products are probably Ga or In, (?P 3/2), +F(?Pi;2). The bearing of the present results 
upon the interpretation of the electronic levels of TIF and on the dissociation energy of AIF 
is briefly considered. 


On the Utilization of Multiple Scattering Measurements, by D. S. SCHONLAND. 


ABSTRACT. ‘The determination of the product pv for particles leaving tracks in photo- 
graphic emulsions is discussed on the basis of measurements of the scattering in the track 
made by the ‘coordinate’ method (Fowler 1950), allowing for errors of measurement due 
to grain size. The use of Fowler’s method with the comparatively large interval lengths 
used in practice is shown to be justifiable in spite of its neglect of much information 
contained in the track. (This neglect increases the root mean square error in the estimation 
of the scattering over the minimum attainable error by a factor less than approximately 
2:6.) A possible method of estimating pv, which approximates to the ‘best possible’, 
based on the use of mean squares, is considered, and the root mean square relative error in 
the estimation found, the necessity for cutting off the scattering distribution at large angles 
being taken into account. A possible method of estimation based on the use of mean moduli 
is also discussed. 


Decay Schemes of *'Sm and Eu, by H. W. WILSON and G. M. Lewis. 


ABSTRACT. Examination of the relatively long-lived isotopes produced by pile irradiation 
of samarium using the proportional tube spectrometer leads to the elucidation of the decay 
schemes of 'Sm and Eu. ‘The conclusions are supported by experiments involving 
fb—e~ and f-y coincidences, utilizing Geiger—Miiller and scintillation counters. 


The Structure of Lithium 6, by A. J. M. Hircucock. 


ABSTRACT. ‘The energies of the lowest states of ®Li have been calculated from the 
standpoint of the shell model, using recent data for the proton—neutron forces. It is con- 
cluded that it is necessary to take into account the tensor force in discussing the lowest states 
of this and other odd—odd nuclei. 


The Photo—Ionization Cross Section of Methane, by A. DALGARNO. 


ABSTRACT. 'The photo-ionization cross section of methane is calculated using 
self-consistent wave functions. ‘The value at the spectral head (105 850 cm~?) is found 
to be 9:4 107" cm?. An approximate general formula is employed to obtain results in~ 
a wave number range up to 200 000 cm”! and the associated oscillator strength is computed. 
It is pointed out that exchange effects probably reduce the cross section. With the aid 
of the general formula an estimate is also made of the photo-ionization cross section of argon. 
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Fig. 1. Laue photograph. ‘Tin crystal Fig. 2. Tin crystal grown in 
grown in graphite. sindanyo. 


Figs. 3 and 4. Crystals grown in glass. 
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